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Synopsis 


This thesis concerns design and analysis of parallel algorithms on PRAMs and 
the parallel comparison model. The following problems are investigated: 

1. Comparison-based merge sorting 

2. Multiway merging 

3. Optimal sublogarithmic time 3- vertex-colouring of rooted forests 

4. 3- vertex-colouring of 4-clique free 3-degree graphs 

5. Vertex-colouring of interval graphs 

6. Edge- colouring of graphs 


The multiway merging problem, where k non-empty sorted arrays of a total 
size n have to be merged into a single array, is studied as a generalisation of sorting 
and merging; when k = 2 we have merging and when k = n we have sorting. 

For comparison-based multiway merging, the existing CREW PRAM 
upper bound of 0(log n) [97] is improved in two ways: 

• an O(logfc 4- log log n) time CREW PRAM optimal algorithm is obtained, and 


• an O(logn) time EREW PRAM optimal algorithm is obtained. 


While the former improves the time bound, the latter weakens the model from CREW 
to EREW. For both of these results, matching lower bounds are also obtained. For the 
CRCW PRAM model, an algorithm that runs in 




time, with nr processors, 2 < r < n, is presented. Also, it is shown that the lower 
bound on time for comparison-based k - way merging, with nr processors on a CRCW 
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PRAM, is 

°(-& + g + E&)- 

Combining all these results, the parallel time complexity of comparison-based k - way 
merging is 

0(MERGE(n, nr) + Sort (k, kr )) 

on all of EREW, CREW and CRCW PRAM models, with nr processors; here MERGE(n,p) 
(resp. SoRT(n,p)) denotes the parallel time complexity of 2- way merging (resp. sorting) 
n items with p processors. 

Multiway merging of integers is also studied. For k - way merging n inte- 
gers drawn from the range [0 ... m — If, an ARBITRARY CRCW PRAM algorithm is 
given; this algorithm, with p processors, runs in 

0 j~fc + loglog(min{n,logm}) + min jlog k + — — ; log log mn + ra ^°S^°S TO | 

time. It is also shown that when the integers are from the range [0. . . j — 1], k - way 
merging can be done in 

o(- + aWk) + , Jssin + log JffiU) 

\p log log (pkjn) log (pk/n)J 

time, with p processors. 

When the input items are single bit integers, the lower bound on time required 
for k - way merging, with p processors, on a CRCW PRAM, is shown to be 
this lower bound is tight. However, when the integers are from the range [0 ... f — 1], 

%o 8 lo° g f P W is a lower bound - 

With the improved parallel multiway merging routine as the basic building 
block, a new comparison-based parallel merge sorting algorithm is obtained. This 
algorithm has small constant factors, and is very different from Cole’s merge sort [19] 
in that it does not use a pipelining scheme. On the parallel comparison model, with 
nr processors, 2 < r < n, the algorithm sorts n items in 

71 _ 2^(105’ (log n/ log r)) 

logr 

time— that is, faster than Cole’s merge sort except for very small values of r. On this 
model, fl(logn/logr) is a lower bound [9] on time for sorting with nr processors; the 
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only upper bound (Alon, 1987) that matches this lower bound, being an adaptation 
of the AKS sorting network, has large constant factors. On a CREW PRAM, with n 
processors, the algorithm runs in logn • 2 0 ( lo s’ n ) time. On a CRCW PRAM, with nr 
processors, the algorithm runs in ^ t s o ^ r • 2°( lo s*( lo s n / lo s r ')) time. 

For the problem of 3-colouring a rooted forest, an O((log log n) log* (log* n)) 
time, optimal parallel algorithm is presented;, TOLERANT CRCW PRAM is the 
model used. For this problem, an O(log(log* n)) time, n processors, suboptimal, CREW 
PRAM algorithm has been known for long [48]. But a sublogarithmic time, optimal 
parallel algorithm has not hitherto been known, even on a CRCW PRAM. 

Furthermore, it is shown that if /(n) is the running time of the best known 
algorithm for 3-colouring a rooted forest on a COMMON or TOLERANT CRCW 
PRAM, a fractional independent set of the rooted forest can be found in 0(/(n)) time, 
with the same number of processors, on the same model. 

Using these results, it is shown that decomposable top-down algebraic compu- 
tation, and hence depth computation (ranking), 2-colouring and prefix summation on 
rooted forests can be done in 0(log n ) optimal time on a TOLERANT CRCW PRAM. 

These algorithms have been obtained by proving a result of independent inter- 
est, one concerning the self-simulation property of TOLERANT: an JV-processor 
TOLERANT CRCW PRAM that uses an address space of size 0(N ) only, can be 
simulated on an n-processor TOLERANT PRAM in 0(%) time, with no asymptotic 
increase in space or cost, when n = 0(1V/ log log N). 

By Brooks’ theorem, for A > 3, any (A + l)-clique-free A-degree graph 
(called a Brooks’ graph) is A-vertex-colourable. An algorithm for A-colouring a A- 
degree Brooks’ graph in 0(A 2 log A logn) time, with n/logn processors, on an EREW 
PRAM, is obtained. In particular, a 3-degree Brooks’ graph is 3-coloured in O(logn) 
time, with n/logn processors, on an EREW PRAM. The basic idea of this algorithm, 
in the context of a CREW PRAM, was treated in the author’s M. Tech, thesis./ 

Besides, the following combinatorial result is obtained: for any two vertices 
u and v that are more than a constant distance apart in a 3-degree Brooks’ graph G, 
there exist two distinct 3-colourings of G of which one has u and v coloured the same, 
and the other has u and v coloured differently. Thus, no vertex in a Brooks’ graph 
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can force the colour of vertices arbitrarily far away; in contrast, while 2-colouring a 
linked list, fixing the colour of one vertex fixes the colour of all others. This highly 
local nature of the problem can be seen as suggesting that on a CRCW PRAM, an 
fl(log nj log log nj time lower bound may not hold. 

The complexity of minimally colouring an interval graph that has a 
known interval representation (let us denote this problem by IGC) on a bounded 
fan-in circuit is studied. The following results are obtained: 

• The problem of 3-colouring a linked list is Nonreducible to IGC. 

• When the chromatic number of the input interval graph is restricted to at most 
a constant, IGC is in NC 1 . 

Using these two results, it is shown that a linked list that has at most a constant 
number of stretches, can be 3-coloured in NC 1 ; here a linked list is visualised as being 
constituted of alternating stretches of forward and backward pointers in an array. This 
complements the observation that 3-colouring a linked list, in general, is unlikely to be 
in NC 1 [72]. Note that in view of this observation, the NC 1 -reduction from 3-colouring 
a list to IGC implies that IGC, in general, is unlikely to be in NC 1 . 

The complexity of IGC on a PRAM is also studied. 

For interval graphs of o(log n) chromatic number, a o(log n) time, polynomial 
processors, CRCW PRAM algorithm is obtained. In particular, when the chromatic 
number is 0((logn) 1_c ), 0 < e < 1, the algorithm runs in 0 (log n /log logn) time. An 
0(log n) time, 0(n) cost, EREW PRAM algorithm is found for general interval graphs. 

Complementing these algorithms, the following lower bound result is obtained: 
even when the left and right endpoints of the intervals are separately sorted, IGC needs 
f2(logn/ log logn) time, on a CRCW PRAM, with a polynomial number of processors. 
This result assumes significance, because, the ft(logn/ log logn) time lower bound [16] 
on finding the chromatic number x °f an interval graph is not valid when either the 
left and right endpoints are separately sorted, or x is also taken as a parameter; in the 
former case x can be found in 0(1) time, and in the latter in 0(logx/loglogra) time, 
with a polynomial number of processors. 


By Vizing’s theorem, any A-degree graph is (A -f l)-edge-colourable. But 
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(A 4- l)-edge-colouring of arbitrary graphs in parallel has proved to be difficult; NC 
algorithms are known only for graphs with A at most a polylogarithmic in n [60]. 
To obtain efficient parallel edge-colouring algorithms, thus, more colours may have to 
be allowed. Probing in this vein, the following EREW PRAM algorithms for edge- 
colouring a A-degree graph are obtained: 

• an algorithm that finds a (A + d)-edge-colouring, 1 < d < A, in 0((logd + 
(A/d) 4 )log 2 n) time, using n + m processors 

• an algorithm that finds a A 1+e -edge-colouring, 0 < e < 1, in 0(log Alog(log“ n)) 
time, using nA 1+£ processors 
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Chapter 1 

Introduction 


With the speed of the present day processors approaching the maxim um phys- 
ically possible, we can no longer expect technological innovations to provide significant 
improvements in the speed of single processor systems [51]. Meanwhile, with the in- 
creasing use of computers in complex computational tasks, the demand for higher 
performance in computing is also on the increase. To meet this demand, thus, it is 
necessary to concentrate on parallel processing. If we were to employ p processors, 
p > 1, concurrently in solving computational tasks, even though inherent sequentiali- 
ties of the tasks may prevent an ideal speed up of 0(p), many times a non-trivial w(l) 
speed up would still be achievable. Hence, it is to be expected that in near future, 
general purpose computers may well be parallel processors. Hence the research interest 
in parallel processors and algorithms. 

Over the recent years, a large body of knowledge on design and analysis of 
parallel algorithms has developed. This thesis is an effort to contribute to the same. 

1.1 Models of Parallel Computation 

The concerns of parallel algorithm design are often very different from those 
of sequential algorithm design. The performance of a parallel algorithm depends on a 
variety of factors like processor allocation, synchronisation and resource sharing, which 
are not present in sequential computation. Because of these, a consensus has not 
yet been reached on the ideal model of computation to be used in designing parallel 
algorithms. The Parallel Random Access Machine (PRAM) is probably the most widely 
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used model of parallel computation.; in this thesis, we shall be mostly concerned with 
designing algorithms on this model. 

A PRAM [61, 55] has p synchronous processors, all having access to a common 
memory. Each processor is similar to a random access machine, the standard model of 
sequential algorithm design [2], and is assigned a unique index from the range [1 . . .p]. 
PRAM models can be classified according to the constraints they impose on the global 
memory access. An Exclusive Read Exclusive Write (EREW) PRAM does not allow 
simultaneous access by more than one processor to the same memory location, for 
read or write purposes. A Concurrent Read Exclusive Write (CREW) PRAM allows 
simultaneous access for reads, but not for writes. A Concurrent Read Concurrent Write 
(CRCW) PRAM allows simultaneous access for both reads and writes. Based on the 
write conflict resolution schemes used, CRCW PRAMs can be further sub-classified. 
PRIORITY CRCW PRAM: in any concurrent write, the processor with the lowest 
index succeeds. 

ARBITRARY CRCW PRAM: in any concurrent write, an arbitrary one of the proces- 
sors succeeds; the programs should work irrespective of the identity of the successful 
processor. 

COMMON CRCW PRAM: all processors involved in a concurrent write should write 
the same value. 

COLLISION CRCW PRAM: when multiple processors write in the same memory lo- 
cation a special collision symbol appears in that location. 

TOLERANT CRCW PRAM: when multiple processors write in the same memory 
location, the content of that location remains unchanged. 

A PRAM model is said to be self-simulating , if for all N > n > 1, a PRAM of 
that model of size n can simulate a single step of another PRAM of the same model of 
size N in O(^) time. All CRCW PRAMs that are at least as powerful as COLLISION 
or COMMON are self-simulating. TOLERANT is not known to be self-simulating [39] 
(see also [4]). 

In Chapter 4, we show that a TOLERANT PRAM of size N with a linear 
address space, can be slowed down by any factor A = fi(loglogiV), with no asymptntir 
increase in space or cost. 

The cost or work of an algorithm that runs in t time using p processors i 
time-processor prc ' 
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If Seq(n) is the worst-case running time of the fastest known sequential algo- 
rithm for a problem of size n, an optimal parallel algorithm for the same problem runs 
in 0(Seq(n)/p) time using p processors. 

Solutions for some algorithmic problems like sorting and merging are typically 
dominated by comparisons between items drawn from a linearly ordered set. For these 
problems, it is useful and instructive to design algorithms on the parallel comparison 
model [96]. In this model, there are p processors synchronised by a common dock. 
In each step of the dock, a processor can compare two items and find their relative 
ordering. Thus, in each step of an algorithm, comparisons are done between at most 
p pairs of (not necessarily distinct) items. The algorithm incurs no cost in using the 
results of prior comparisons to decide which comparisons to make next. As soon as 
enough ordering rdations have been found to formulate the output as required, the 
algorithm ter mi nates. The number of steps executed till then is the running time of 
the algorit hm . This model is clearly more powerful than all standard PRAM modds. 

Another model helpful in studying comparison intensive problems on items 
from a linearly ordered set is the comparator network. A comparator is a two-input 
module with an ordered pair of outputs: if a and b are the inputs, then the first output 
is min{a, &}, and the second output is max{a,6}. A comparator network is made only 
of comparators; its size is the number of comparators in it, and its depth is the length 
of the longest path in it from an input to an output. This model is clearly weaker than 
all standard PRAM models. 

A boolean circuit [61] is a directed acyclic graph in which each node is an 
input, or an output, or a 1, or a 0, or an AND gate, or an OR gate, or a NOT gate. For 
a boolean circuit, the size is the number of edges in it, and the depth is the length of the 
longest path in the circuit from an input to an output. When the in-degree (fan-in) of 
a node is restricted to a constant the circuit is said to have a bounded fan-in, otherwise 
an unbounded fan-in. NC*, k > 0 is the class of all problems solvable in O(log*?i) 
depth on a bounded fan-in circuit of a polynomial size. NC is the class U*> 0 NC*. AC*, 
k > 0 is the class of all problems solvable in O(log* n) depth on an unbounded fan-in 
circuit of a polynomial size. AC is the class Ufc>oAC*. It is known that NC = AC. 
Moreover, 

NC* C EREW* C CREW* C CRCW* C AC* C NC* +1 

FL [24] is the class of problems solvable in deterministic log space; here a 
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problem is a multiple valued function and solving a problem involves finding any one 
of its possible values for a given input. It is known that NC 1 C FL, but whether the 
inclusion is proper or not is still not known. 


1.2 Sorting, Merging and Graph Colouring 

This thesis is mainly concerned with sorting, merging and graph colouring. 

Sorting is a fundamental tool for algorithm designers. A large percentage of 
computer time is typically spent in sorting [9, 64]. So, this problem has been, and still 
is, intensely researched. While dealing with sorted lists, time and again, it becomes 
necessary to merge several into one. Conversely, merging has been used as a basic tool 
in designing sorting algorithms: merge sort has been one of the first O(nlogn) time 
comparison-based sequential sorting algorithm to be discovered, whereas Cole’s merge 
sort is the only optimal and yet simple comparison-based parallel sorting algorithm. 

Graph colouring assumes special significance in the parallel world because of 
its applications in symmetry breaking and scheduling. For example, a context, where 
a number of entities vie for a number of shared resources, can be captured in a conflict 
graph, and a -colouring of this graph can give a schedule for resource allocation. 

1.2.1 Sorting and Merging 

Sorting is the problem of arranging a set of elements drawn from a linearly 
ordered set in nondecreasing or non-increasing order. Merging, on the other hand, is 
combining two sorted lists into a single one. In the multiway merging problem, we are 
given k sorted arrays, and we have to merge them to get a single sorted array. When 
the value of k is understood, we shall also call this problem k - way merging. When 
k = 2, k - way merging degenerates into ordinary merging. The special case of k = n, 
where n is the total number of elements, is sorting. 

When the input elements are considered atomic entities, comparison being the 
only operation performable on them, sorting, 2- way merging and multiway merging are 
said to be comparison-based. 

Sorting n elements, using only comparisons, requires fl(ralogn) time sequen- 
tially, and there are several algorithms that match this lower bound [2, 64]. In the 
parallel setting, Batcher gave 0(log 2 n) depth, 0(n log 2 n) cost comparator networks 
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for sorting [64]. Finding a logarithmic depth network for sorting remained an open 
question for long, until Ajtai, Komlos and Szemeredi [3, 81, 82] found the O(logn) 
depth, O(nlogn) size, AKS sorting network; though this network achieves an O(logn) 
depth, the constant factor hidden by the O-notation here is quite large. Cook, Dwork 
and Reischuk [23] showed that on a CREW PRAM, finding the OR of n bits requires 
fl(log n) time, with any number of processors; the same lower bound holds for sorting as 
well. So, a simulation of the AKS sorting network is optimal on all models that are not 
more powerful than the CREW PRAM model. Cole gave a simple, logarithmic time, 
linear processors, EREW PRAM sorting algorithm that has small constant factors [19]. 
For the parallel comparison model, Azar and Vishkin [9] have shown that, sorting n 
elements with p > n processors requires at least comparison steps. The only 

upper bound (due to Alon, [9]) that matches this lower bound, being an adaptation of 
the AKS sorting network, has large constant factors. Cole’s merge sort runs on the par- 
allel comparison model in CKk ^f/n) ' -loglog(p/n)) steps. For the CRCW PRAM model, 
Beame and Hastad [10] proved that computing the parity of n bits requires 
time with p > n processors. Combining this result with the parallel comparison model 
lower bound, sorting n elements on a CRCW PRAM with p > n processors requires 
^( io^ p/nj + log*logp ) time, because, parity i/ reducible to sorting [10]. No algorithm 
matching this lower bound is known today. The fastest known algorithm takes 


O 


( , lQ , S ? s (loglog(p/ft)) 5 2 Q ( log * n-io S *(p/»)+D + 
\iog(p/n) 


logn 
log log p 


time with p > n processors [45]. Cole’s merge sort runs on a CRCW PRAM in 

0 ( log log(p/ n) ) Re- 
sorting algorithms have been designed on other models of parallel compu- 
tation too. In particular, on the parallel comparison model, Columnsort of Leighton 
[68] runs in 0((logn/logr) 3 - 419 ) steps, and Cubesort of Cypher and Sanz [25] runs in 
O(2 lo8 *( n / r )(logn/logr) 2 ) steps, with 2 < r < n processors. - 

Comparison based 2-way merging, in contrast, takes only linear time sequen- 
tially. In the parallel setting, Batcher gave merging networks of O(logra) depth and 
O(nlogn) size [64]. These networks are asymptotically the best possible, because, 
any merging network must have fl(nlogu) comparator modules [64]. However, on 
an EREW PRAM we can do better: Hagerup and Rub [46] describe how to obtain 
an O(logn) time optimal algorithm using Batcher’s networks. Also, for the EREW 
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PRAM model, a reduction from searching to 2- way merging, along with the logarith- 
mic time lower bound for searching [93], implies that 2- way merging requires fi(logn) 
time with any number of processors. Borodin and Hopcroft [14] showed that on the 
parallel comparison model, with p > n processors, Q(log time is required to 

merge two sorted arrays of a total size n. This matches the best known upper bound 
on the CREW PRAM model [14, 96]. From this suboptimal CREW PRAM algorithm, 
Kruskal has obtained an optimal 0 (log log n) time CREW PRAM algorithm [65]. 

A decision tree argument shows that fc-way merging by comparisons requires 
Q(nlogfc) time sequentially [97]. Also, using a heap for finding the minimum of k 
values, k - way merging can be done in O(nlogk) sequential time [64]. Lee and Batcher 
[67] have generalised the odd-even merge principle to find an O(log n log k ) depth, k - way 
merging network. Wen [97] has obtained an O(logn) time optimal parallel algorithm 
on a CREW PRAM. 

If the input elements are assumed to be integers, the lower bound of fl(nlogn) 
is no longer valid for sorting. The best known upper bound is 0(n loglog n) [7], whereas 
no stricter lower bound than the trivial f 2(n) is known. Parallel integer sorting on a 
CRCW PRAM must obey the time lower bound of Beame and Hastad, be- 

cause of the reduction from parity. Andersson et al. [7] showed that on an ARBITRARY 
CRCW PRAM, n integers can be sorted in O(logn) time, with O(nloglogn) opera- 
tions. Bhatt et al. [12] gave an + log log m) time, 0(n loglog to) operations 

algorithm, for sorting integers in the range [0 . . .m — 1], on an ARBITRARY CRCW 
PRAM. 

For 2- way integer merging, Berkman and Vishkin [78] gave an 0 (log log log m) 
time optimal CREW PRAM algorithm for the case where the integers are from the 
range [0 ... m- 1]. For the case of m = n, they gave an 0(a(n)) time optimal algorithm 
for the CRCW PRAM model, where a is the inverse Ackermann function. Hagerup and 
Kutylowski [43] described an O(log logn + logmin{n,m}) time, optimal cost EREW 
PRAM algorithm for merging integers from the range [0 ... m — 1]. 

In Chapters 2 and 3 of this thesis, the following results on sorting and multiway 
merging are presented. 

For comparison-based multiway merging, the existing CREW PRAM upper 
bound of O(logn) [97] is improved in two ways: 
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• an <3(logA; + loglogn) time CREW PRAM optimal algorithm is obtained, and 


• an O(logn) time EREW PRAM optimal algorithm is obtained. 


While the former improves the time bound, the latter weakens the model from CREW 
to EREW. For both of these results, matching lower bounds are also obtained. For the 
CRCW PRAM model, an algorithm that runs in 


0 ( log §7 + ' (l08l ° 8 r ) s ' 20< ‘° 8 ' r+1 > + 


log k 


loglog(fcr) 


time, with nr processors, 2 < r < to, is presented. Also, it is shown that the lower 
bound on time for comparison-based k- way merging, with nr processors on a CRCW 
PRAM, is 

V 5 log r logr log log (kr) J 

Combining all these results, the parallel time complexity of comparison-based fc-way 
merging, with nr processors, is 


0(Merge(to, nr) + Sort (k,kr)) 


on all of EREW, CREW and CRCW PRAM models; here Merge(to,p) (resp. Sort(to,p)) 
denotes the parallel time complexity of 2-way merging (resp. sorting) to items with p 
processors. 

Multiway merging of integers is also studied. For k- way merging n integers 
drawn from the range [0 ... m — 1], an ARBITRARY CRCW PRAM algorithm is given; 
this algorithm, with p processors, runs in 


0 ( ■ —tt-t + loglog(min{TO,logrro}) + min jlogfc + - ^ — - , log log ttoto iSilSI — jA 

\ log log A; l p V iJ 

time. It is also shown that when the integers are from the range [0 ... j — 1], A;- way 

merging can be done in 


0 



log A; . _logn_\ 

log log(pk/n) ° S log(pkJn)J 


time, with p processors. 

When the input items are single bit integers, the lower bound on time required 
for A:-way merging, with p processors, on a CRCW PRAM, is shown to be 
this lower bound is tight. However, when the integers are from the range [0 . . . £ - 1], 

fl ( cs§fe sj) is a lower bound - 
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Using the improved parallel multiway merging algorithm as the basic building 
block, a new comparison-based parallel merge sorting algorithm is obtained. This is 
presented in Chapter 3. This algorithm has small constant factors, and is very different 
from Cole’s merge sort [19] in that it does not use a pipelining scheme. On the parallel 
comparison model, with nr processors, 2 < r < n, the algorithm sorts n items in 

logn . 2°( lo g*qogVlogr)) 
logr 

time — that is, faster than Cole’s merge sort except for very small values of r. On a 
CREW PRAM, with n processors, the algorithm runs in logn • 2°fi°8* n ) time. On a 
CRCW PRAM, with nr processors, the algorithm runs in • 2 °( lo s*flogn/lo g r)) 

time. 

1.2.2 Graph Theoretic Terms and Notations 

We follow the graph theoretic notations of Harary [50] and Golumbic [36]. 

A graph G = (V, E) consists of a finite nonempty set V of vertices and a set 
E C V X V of unordered pairs of distinct vertices [50]. Each e € E is an edge of G. 
Throughout this thesis, we shall denote the number of edges | E | by m, the number 
of vertices | V| by n. Two vertices u, v € V are said to be adjacent iff {«, u} € E; in 
this case, we also say that u and v are neighbours of each other. We denote the set 
of neighbours of v by N(v). When e = {u, v}, we say that e is incident with u and v. 
The number of neighbours of a vertex is its degree. The maximum vertex degree of a 
graph G is denoted by A(G); whenever there is only one graph G under consideration 
we shall use A for A(G). A graph in which every vertex is of degree d, is a d-regular 
graph. A 3-regular graph is also called a cubic graph. A graph that has a maximum 
vertex degree of 3 is called a subcubic graph. 

In a directed graph every edge is an ordered pair; if (u, v) 6 E we say that v 
is an out-neighbour of u and u is an in-neighbour of v. The number of out-neighbours 
of a vertex is its out-degree, and the number of in-neighbours is its in-degree. 

G' = (V', E') is said to be a subgraph of G if V' C V and E' C E. If 
E' = {{u, v} | {u,v} € E and u,v € V '} then, G' = G[V'] is the subgraph induced by 
V‘. We use G - U to denote G[V - U] for U C V, G - v to denote G - {n} for v € V, 
and G — e to denote G' = (V, E - {e}) for e € E. 
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A graph is called complete if for every pair (u, v ) of distinct vertices, u and 
v are adjacent. A complete graph on n vertices is called an re-clique. A maximum 
clique of G is a largest cardinality complete subgraph of G; its cardinality is the clique 
number v(G) of G. 

An independent set / of G is a subset of V such that for every u, v 6 I, {re, re} 
is not in E. In a maximal independent set (MIS for short) I , in addition the following 
condition holds: for each v 6 V either v € I or v has a neighbour in I. 

In the vertex colouring problem on graphs, we seek to assign colours to the 
vertices so that no two adjacent vertices get the same colour. The minimum number 
of colours needed to colour G is the chromatic number of G and is denoted by x{G). 
Clearly, x(G) is not smaller than ui(G), the clique number. If, for every V' C V, 
x(<J f [V r/ ]) = w(G[V'']), then G = (V, E) is called a perfect graph. 

A graph G = (V, E) is called an interval graph, if for some set of intervals 
of a linearly ordered set, there is a bijection / : V — *■ 5 so that two vertices u and re are 
adjacent in G iff f(u) and /(re) overlap. An interval graph is known to be perfect [36]. 

A bipartite graph G is a graph whose vertex set V can be partitioned into 
two subsets Vi and V? so that every edge in G has one endvertex in Vj and the other 
endvertex in V<i. A bipartite graph is 2 vertex-colourable. Any graph without cycles is 
a forest. Each connected component of a forest is a tree. A rooted forest F = (V, E) is 
a directed forest in which the out-degree of each vertex is at most one. Roots in F are 
the vertices with out-degree exactly zero. Each edge in E is from a child to its parent; 
that is, if ( u , v) € E, v is the parent of u, and u is a child of v. Forests are bipartite 
graphs. 

In the edge-colouring problem on graphs, we seek to assign colours to the 
edges so that no two edges incident with the same vertex get the same colour. For a 
graph G = (V,E), the minimum number of colours required to edge-colour G is the 
chromatic index x'{G) °f G. Clearly, x' ^ A. 

1.2.3 Vertex Colouring 

Finding the chromatic number of a 4-degree planar graph, let alone an arbi- 
trary graph, is NP complete [29]. Minimally vertex colouring an arbitrary graph is, 
thus, computationally difficult. 

Hence, to colour graphs efficiently we should either restrict our attention to 
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certain types of graphs that are easy to colour minimally, or settle for a possibly non- 
minimal colouring. 

In the sequential setting, both these approaches have given efficient algo- 
rithms. A simple 0(m) time search can 2-colour a bipartite graph. The proof of 
4-colour theorem gives a polynomial time algorithm for 4-colouring a planar graph; 
this algorithm has large constant factors. But planar graphs can be 5-coloured in 0(n ) 
time: see [42] for a survey of such algorithms. Several special classes of perfect graphs 
have efficient algorithms for minimal colouring, though no polynomial time algorithm 
for recognising general perfect graphs is known. For example an interval graph can be 
recognised and minimally coloured in 0(n) time. A graph can be easily (A-f-l)-coloured 
in 0(m) time. Brooks showed that any connected graph is A-colourable if it is neither 
a complete graph on A + 1 vertices nor a circuit of odd length [99, 73]. A A-colouring 
has often been called a Brooks’ colouring and a A-colourable graph a Brooks’ graph. A 
Brooks’ graph can be A-coloured in linear time (see [73], problems 9.12 and 9.13). 

In the parallel setting, a bipartite graph can be 2-coloured in O(logra) time. 
In particular, a linked list can be 2-coloured in O(logn) time optimally on an EREW 
PRAM [6]. A rooted forest can be 3-coloured in O(log(log* ra)) time with n processors 
on an EREW PRAM [48] Planar graphs can be 5-coloured in 0(logralog* n) time and 
O(n) work on an EREW PRAM [42]. Given an interval graph in adjacency list form, 
an interval representation for it can be found in O(log 2 n) time using n + m processors 
on an EREW PRAM [63]. Given an interval representation, the interval graph can be 
minimally coloured in O(logn) time with n processors on an EREW PRAM [89]. 

Closely associated to (A + l)-colouring is the MIS problem; this is because, 
there is an NC reduction from (A+l)-colouring to the MIS problem [74]; thus an NC- 
algorithm for the MIS problem implies an NC-algorithm for (A-f l)-colouring as well. 
Karp and Wigderson [62] describe an 0((logn) 4 ) time algorithm for the MIS problem, 
with 0((j^) 3 ) processors on an EREW PRAM, thus proving the problem to be in NC. 
Alon, Babai and Itai [5] and Luby [74] have obtained fast and simple randomised algo- 
rithms for the MIS problem; these algorithms take O((log n) 2 ) expected time with 0(m) 
processors on an EREW PRAM. Alon et.al. [5] and Luby [74] also present derandomi- 
sation techniques to obtain fast deterministic versions of their respective algorithms. 
The first linear processors NC algorithm for the MIS problem is due to Goldberg and 
Spencer [34], and takes 0((logn) 4 ) time on an EREW PRAM; they later improved 
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the time bound to O((logn) 3 ) [33]. Luby [75] describes a derandomisation technique 
without processor penalty; this results in a linear processors O((log ti) 3 log log re) time 
algorithm for (A+l)-colouring on a CREW PRAM. 

For the special case of A being bounded, [32] gives an 0(n) processors 0(log* n) 
time algorithm for the EREW PRAM model whereas an optimal algorithm that takes 
O(log (fc) n) time with - processors, again on an EREW PRAM, is given in [87]; 
for general graphs these algorithms take, respectively, 0(Alog A(log* n + A)) time and 
0(A 2 log A log^ n) time. 

The problem of Brooks’ colouring, is known to be in NC [57, 59, 47, 80]. 
Panconesi and Srinivasan [80] prove a 0(log A re) bound on the search radius required 
to fix the colour of a lone uncoloured node of a Brooks’ graph. Using this result, they 
obtain an time reduction from Brooks’ colouring to to (A-f-l)-colouring [80]. 

On bounded degree graphs their algorithm can be implemented with 0(n ) processors 
in 0(log 2 nlog* n) time on am EREW PRAM. 

For Brooks’ colouring of cubic graphs (i.e., A = 3), Karloff [58] give an algo- 
rithm that can be implemented in 0(log 2 n) time with processors on an EREW 
PRAM. Furthermore, Karloff [58] shows that Brooks’ colouring an o- degree graph G 
can be reduced in NC, to Brooks’ colouring an (a-l)-degree subgraph of G, for a > 3. 
This reduction requires a spanning tree in an 0(n)- vertex graph to be found, and so, 
even for bounded degree graphs, can at best be implemented in 0 (log re) time with 
(n+ To ) g a l m ’ n) processors on a CRCW PRAM [21, 54] and in 0 (log nlog log n) time with 
n processors on an EREW PRAM [17]. So, on bounded degree graphs, the conse- 
quent Brooks’ colouring algorithm, in [58] takes 0(log 2 n) time with n processors on an 
EREW PRAM and 0(log 2 n ) time with processors on a CRCW PRAM. 

In Chapters 4, 5 and 6 of this thesis, the following results related to vertex 
colouring are obtained. 

For the problem of 3-colouring a rooted forest, an 0((loglogn)log*(log* n)) 
time, optimal parallel algorit hm is presented, in Chapter 4; TOLERANT CRCW 
PRAM is the model used. For this problem, a sublogarithmic time, optimal paral- 
lel algorithm has not hitherto been known, even on a CRCW PRAM. Furthermore, it 
is shown that if f{u) is the running time of the best known algorithm for 3-colouring a 
rooted forest on a COMMON or TOLERANT CRCW PRAM, a fractional independent 
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set of the rooted forest can be found in 0(f(n)) time, with the same number of proces- 
sors, on the same model. Using these results, it is shown that decomposable top-down 
algebraic computation, and hence depth computation (ranking), 2-colouring and prefix 
summation on rooted forests can be done in O(logn) optimal time on a TOLERANT 
CRCW PRAM. 

In Chapter 5, we study Brooks’ colouring. An algorithm for A-colouring a A- 
degree Brooks’ graph in 0( A 2 log A log n) time, with n / log n processors, on an EREW 
PRAM, is obtained. In particular, this algorithm 3-colours a 3-degree Brooks’ graph 
in O(logn) time, with n/logn processors, on an EREW PRAM. 

Besides, the following combinatorial result is obtained: for any two vertices 
u and v that are more than a constant distance apart in a 3-degree Brooks’ graph G, 
there exist two distinct 3-colourings of G of which one has u and v coloured the same, 
and the other has u and v coloured differently. Thus, no vertex in a Brooks’ graph 
can force the colour of vertices arbitrarily far away; in contrast, while 2-colouring a 
linked list, fixing the colour of one vertex fixes the colour of all others. This highly 
local nature of the problem can be seen as suggesting that on a CRCW PRAM, an 
fl(logn/log log n ) time lower bound may not hold. 

In Chapter 6, the complexity of minimally colouring an interval graph that 
has a known interval representation (let us denote this problem by IGC) on a bounded 
fan-in circuit is studied. The following results are obtained: 

• The problem of 3-colouring a linked list is NC 1 -reducible to IGC. 

• When the chromatic number of the input interval graph is restricted to at most 
a constant, IGC is in NC 1 . 

Using these two results, it is shown that a linked list that has at most a constant 
number of stretches, can be 3-coloured in NC 1 ; here a linked list is visualised as being 
constituted of alternating stretches of forward and backward pointers in an array. This 
complements the observation that 3-colouring a linked list, in general, is unlikely to be 
in NC 1 [72]. Note that in view of this observation, the NC 1 -reduction from 3-colouring 
a list to IGC implies that IGC, in general, is unlikely to be in NC 1 . 

The complexity of IGC on a PRAM is also studied. 

For interval graphs of o(logn) chromatic number, a o(logn) time, polynomial 
processors, CRCW PRAM algorithm is obtained. In particular, when the chromatic 
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number is (^((logn) 1 e ), 0 < e < 1, the algorithm runs in O(logn/loglogn) time. An 
O(log n) time, 0(n) cost, EREW PRAM algorithm is found for general interval graphs. 

Complementing these algorithms, the following lower bound result is obtained: 
even when the left and right endpoints of the intervals are separately sorted, IGC needs 
ft(logn/ log log n) time, on a CRCW PRAM, with a polynomial number of processors. 
This result assumes significance, because, the ft(logn/loglogn) time lower bound [16] 
on finding the chromatic number x of an interval graph is not valid when either the 
left and right endpoints are separately sorted, or x is also taken as a parameter; in the 
former case x can be found in 0(1) time, and in the latter in 0 (log x/ log log n) time, 
with a polynomial number of processors. 

1.2.4 Edge Colouring 

The chromatic index x' of a graph G is at least A (G). Vizing [26, 13, 76, 56] 
and Gupta [56] showed that x' < A + 1; this result is popularly known as Vizing’s 
theorem. Thus, the chromatic index of a graph is either A or A + 1; accordingly, 
graphs are said to be in either Class 1 or Class 2. Deciding the class of a graph is known 
as the classification problem. Even when restricted to cubic graphs, the classification 
problem is NP-complete [52]. Optimally edge colouring an arbitrary graph is, thus, 
computationally difficult. 

To get efficient edge-colouring algorithms, hence, we should either restrict our 
attention to certain types of graphs that are easy to edge-colour optimally, or settle for 
a possibly nonoptimal (A+l)-edge-colouring of general graphs. 

On the sequential front, both approaches have given efficient algorithms. Pla- 
nar graphs with A > 8 and bipartite graphs are in Class 1 [26]. A planar graph with 
A > 8 can be A-edge-coloured in 0(n 2 ) time (the result is quoted in [18]), whereas, a 
bipartite graph can be A-edge-coloured in 0 (min {m log n, re 2 log n}) time [20, 28]. An 
arbitrary graph can be (A-l-l)-edge-coloured in 0(min{mrc,mAlog n, my/nTog n }) (the 
result is quoted in [18]). 

In the parallel setting, NC algorithms have been obtained for A-edge-colouring 
planar graphs with large A, as well as bipartite graphs. Planar graphs with A > 19, 
can be A-edge-coloured in 0(log 2 n) time on an EREW PRAM with n processors [18], 
while planar graphs with A > 8, can be A-edge-coloured in 0(log 3 n) time on an 
EREW PRAM with n 2 processors [18]. A bipartite graph can be A-edge-coloured in 
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O(log 2 ralog A) time on an EREW PRAM with n + m processors [69], and 2 -edge- 
coloured in 0(lognlog A) time on an EREW PRAM with n + m processors [69]. 

In contrast, (A-f-l)-edge-colouring of arbitrary graphs is not yet known to be 
in NC; Karloff and Shmoys [60] give an algorithm that is in NC only when A is at most 
polylogarithmic in n. This algorithm runs in 

0( A 5 log n( A + log n + min{log 3 n, A 2 log A(log* n + A 2 )})) 

time on an EREW PRAM with n + m processors. Liang, Shen and Hu [71, 70] give 
an 0(A 4,5 log 3 Alogn + A 4 log 4 n) time, n 3 logA + nA 3 processors CRCW PRAM 
algorithm for (A-f l)-edge-colouring a general graph. 

To edge-colour an arbitrary graph in NC, thus, we may have to use more than 
A + 1 colours. Fiirer and Raghavachari’s results [27] are in this vein. They give an 
O(lognlog A) time, n + m processors, CREW PRAM algorithm for cA-edge-colouring 
a general graph, c > 1, and an O(log* n) time, n + m processors, CREW PRAM 
algorithm for A 2 -edge-colouring a general graph. 

In Chapter 8 of this thesis, the following EREW PRAM algorithms for edge- 
colouring a general A- degree graph are presented: 

• an algorithm that finds a (A + d)-edge-colouring, 1 < d < A, in 0((logd + 
(A/d) 4 )log 2 n) time, using n + m processors 

• an algorithm that finds a A 1+£ -edge-colouring, 0 < e < 1, in 0(log Alog(log* n)) 
time, using nA 1+£ processors 
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Chapter 2 

Multiway Merging in Parallel 

Multiway merging is the problem of merging k sorted arrays into a single 
sorted array. When the value of k is understood, we shall also call this problem k - way 
merging. This problem can be seen as a generalisation of sorting and merging: the 
special case of k = 2 is ordinary merging, and that of k = n is sorting. 

Parallel solutions for the multiway merging problem form the topic of this 
chapter. An algorithm schema for multiway merging is described in Section 2.1. The 
implementation of the schema on various PRAM models for the comparison based 
version of the problem is discussed in Section 2.3; matching lower bounds are also 
presented. These results improve the existing CREW PRAM upper bound of O(logn) 
[97]. Multiway merging of integers is considered in Section 2.4. 

2.1 The Algorithm Schema 

The following definitions are due to Cole [19]. For three keys x,y and z, we 
say that y is between x and z. or that x and z straddle y, if a: < y < z. For a key x and 
a list X, the rank of x in X is the number of elements in X that are not larger than x. 
A list L\ is said to be ranked in another list X 2 if every element of Xi knows its rank 
in X 2 ; Xi and Z 2 are cross ranked if they are ranked into each other. 

The algorithm schema is now described: 

Input: Arrays A; of size n t - for 1 < i < k, each sorted in increasing order, where 
5Zi=i n « = n; p processors are available. Without loss of generality, we assume that all 
elements are distinct, because otherwise, for each i, the j-th element aj- of Aj, can be 
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replaced by a triplet with the linear ordering extended lexicographically. 

Output: A\ U . . . U -Afc sorted in increasing order in an array 2?[1 . . . n]. 

Step 1: For a parameter A to be appropriately chosen, if n < kX, sort the elements 
using the fastest available sorting algorithm with p processors; exit. 

Step 2: Partition A,-, for 1 < i < k, into segments of size A each. Here, by a segment 
we mean a subarray of consecutive locations. (If n,- is not divisible by A, the last 
segment may be of a smaller size.) There are f,- = segments in A,-. So the 
total number of segments t = Yli=i <» < & + x — T* ^or eac ^ segment, let its 
first element be called its leader. Let an array A(, for 1 < i < k, contain the set 
of leaders from Ai arranged in increasing order. 


Step 3: For 1 < i,j < k, i ^ j, cross-rank the leader arrays A\ and A'-. Each leader 
x now has k ranks, r,-(x) for 1 < * < k, assigned to it, r,-(x) being the number of 
elements less than or equal to x in A\. 


Step 4: For each leader x, in parallel, add together all r,-(x)’s. Consequently, each 
leader knows the total number of smaller leaders over all arrays; that is, now, all 
leaders can be placed in an array L of size f, sorted in increasing order. 


Step 5: Each leader x knows the two consecutive leaders of A,-, for 1 < i < k, that 
straddle x; note that when x £ A,-, the smaller of the two leaders is x itself. In 
other words, x knows the segment s,(x) of A,- to which it belongs. For each x € L 
and for each array A,-, search s,-(x) in parallel to find the two consecutive elements 
in it that straddle x. Let 22, (x) be the index in A,- of the smaller one of these 
two. Thus, Ri{x) is the number of elements in A,- smaller than or equal to x. 

Step 6: For each leader x, compute J(x) = ELi R i( x )> whidl clearl y is the index of 
x in jB, the output array. 

Remark: For any two consecutive elements x and y of L, there can be at most A 
elements in each A, - , 1 < * < fc, that x and y straddle. That is, Hi(y~) — 22,'(x) < A, for 
all i. Hence, I(y ) — I(x) < fcA. 


Step 7: For every pair of consecutive elements x and y of L, let 


5,-(x) 


Ri(y) - 22, (x) when y <£ A,- 
22, (y) - 22i(x) - 1 when y € A,- 
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Note that 8{(x ) is the number of elements that are larger than x but smaller than y 
in A{. Let (Ai(x), . . ., Afc(x)) be the prefix sums computed over (£i(x),. . . ,^(x)); 
Ao(a:) is set to zero. 

Step 8: Copy A,-[iZ,(x) + 1 . . . 12;(x) + £,(x)] to 5 c [A,_x(x) + 1 . . . A,-(x)]. That is, we 
isolate all elements that are larger than x but smaller than y into an array S x . 

Step 9: Sort the elements of S x for each x € X, in parallel, using p processors overall, 
and place the output in consecutive locations of the subarray of B that begins at 
index 7(x) + 1. 

Remark: Note that N(x) =| S{x) |= I(y) - I(x) - 1 < kX. 

2.2 Preliminaries 

The time needed to solve instances of size n, respectively of, prefix summation 
of b bit numbers, merging, searching, sorting and fc-way merging, using p processors in 
each case, will be denoted by PREFix(n,p, b), MERGE(ra,p), SEARCH(n,p), S0RT(n,p) 
and fc_WAY_MERGE(n,jp). 

We use in the sequel, the following upper bounds that are well known in 
literature. 

Upper Bound 2.1 [55] When comparison is the only operation allowed on the keys, 
on a CREW PRAM, MERGE(n,p) = 0(loglogn + |). 

Upper Bound 2.2 [55] When comparison is the only operation allowed on the keys, 
with p > n, on a CREW PRAM, Merge(ti,p) = O(log 

Upper Bound 2.3 [78] When the keys are integers drawn from the range [0 . . . m — 1], 
on a CREW PRAM, Merge(tz,p) = 0(logloglogm + j). 

Upper Bound 2.4 [78] When the keys are integers drawn from the range [0 . . . n — 1], 
on a CRCW PRAM, Merger, p) = 0(a(n) + f). 

Upper Bound 2.5 [46] On an EREW PRAM, MERGE(n,p) = 0(logn + *j). 

Upper Bound 2.6 [22, 91, 41] On a CRCW PRAM, 

f f b . ) log n n\ 

PREFIx(n, P ,6) = 0 (min (logj^.log"} + + ~j ■ 
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Upper Bound 2.7 [55] On an EREW PRAM, Prefi x(n,p,£>) = 0(logn + £). 

Upper Bound 2.8 [55] On a CREW PRAM, SEARCH(n,p) = 0 (\ ||^). 

Upper Bound 2.9 [3, 81, 19] When comparison is the only operation allowed on the 
keys, on an EREW PRAM, 

SoRT(n,p) = 0 (logn + -- " g7Z j . 

Upper Bound 2.10 [45] When comparison is the only operation allowed on the keys, 
on a CRCW PRAM, 

S 0 RT(n, P ) - O ( i ^)(loglo S ( I ./n)) 5 2O^-“-W0./.)«) + -fe.) . 

Upper Bound 2.11 [7] When the keys can be treated as integers, on an ARBITRARY 
CRCW PRAM, Sort (n,p) = 0(logn + ^g - lo s n ). 

Upper Bound 2.12 [12] When the keys are integers drawn from the range [0 . . . m-1], 
on an ARBITRARY CRCW PRAM, SoRT(n,p) = 0 (j^L 2 L_ + loglogm + 

Upper Bound 2.13 [55] The minimum ofn numbers can be found in 0(1) time, with 
0(n 1+c ) processors, on a COMMON CRCW PRAM, for any fixed e > 0. 

Following [40], we say, a set of n objects drawn from a linearly ordered set are 
padded-sorted if they are arranged in an array of size 0{n) in non-decreasing order. 

Upper Bound 2.14 [45] A set ofn integers drawn from the range [0 . . .n - 1] can be 
stably padded-sorted with a constant padding factor in 0 ((logn) 1 / 2 (loglogn) 3 / 2 ) time 
and 0(n log log n) operations on an ARBITRARY CRCW PRAM. 

We shall encounter several times in the course of this chapter, the following 
allocation problem 7 Z. Given an array A = (nj,. . . ,njt), where J2i=i n i = n, fill an 
array P[l...n] so that, for 1 < * < n, P[i] = j iff Y? x = o n * < * < o n xi where 
no = 0 . 

When p processors are available, the problem 71 can be solved as follows. Find 
the prefix sums ( 771 , . . . , T]k) of A. Let B[i] = ( i , 1 ), and C\j] = (rjj, 2 ) for 1 < i < n and 
1 < y < fc. Cross rank B and C. Now, P can be filled as required in the problem: if 
the rank of P[i] in C is j then P[i] = j. Hence we have the lemma below: 
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Lemma 2.1 The problem TZ can be solved in 0(PREFlx(A:,p,logn) -f MERGE(n,p)) 
time. 


We assume that on all computational models under consideration here, the 
upper bounds for merging, prefix summing and sorting satisfy the the “regularity” 
condition that larger problem instances with the same processor-advantage ( p/n ) will 
take more time. That is, Merge (n u ni/3) < cMerge(u 2 , n 2 /i), for n x < ra 2 and a 
constant c > 0; similarly for sorting and prefix summing. 

2.3 Comparison Based Multiway Merging 

2.3.1 The schema on PRAMs with Concurrent Reads 

First, we consider some aspects common to CREW and CRCW PRAM im- 
plementations of the algorithm schema presented in Section 2.1. Let us set A = k 2 . 

In Step 1, where we are sorting all the elements together, n < k 3 . Since, 
solving a larger problem instance with the same processor advantage cannot be easier, 
the time taken by this step is 0(SoRT(n,p)) = 0(Sort(A: 3 , 2 ^-)). 

We can both partition the input arrays and form the leader arrays in 0(1) 
time, using n processors, once they are assigned one per element. But here we have p 
processors. Assign one processor for every Q(^) elements of the input, using Lemma 2.1. 
Each processor can now sequentially select leaders from the elements assigned to it. So, 
Step 2 takes + Merge(u,p) + PREFlx(fc,p,logn)) time. 

In Step 3, merges, each of size at most t, are to be performed in parallel. 

Assign 0($-) processors per merge instance. Solve each instance independently. The 
time taken is clearly 0(MERGE(t, ^-)). 

Steps 4 and 6 each requires every leader to perform an addition of size k, of 
log n bit numbers. Use Lemma 2.1 to assign O(^) processors per leader. Find the sums 
using a prefix summation algorithm in 0(Prefix(&, f,logn)) time. The total time 
taken is 0(Merge(h,p) + PREFlx(fc, |,logn)). 

Each leader has to perform k searches in Step 5, each over a segment of size & 2 ; 
that is, a total of tk searches have to be performed in parallel. At most [^J processors 
are available per search. If k processors were to be used per search, the time taken would 
be SEARCH(fc 2 , ife), which is 0(qg£) = 0(1), by Upper Bound 2.8. With processors 
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per search, hence, we can finish the step in 0(jj^ + Search (k 2 , k )) time. As tk 2 < 2n 
(see Step 2 of the schema), the time needed for Step 5 is 0(f + Search(A: 2 , k)) = 0(f). 

Step 7 requires every leader to perform an addition of size k , of log k bit 
numbers. With a processor allocation similar to that in Step 4, the sums can be found 
in 0(PREFlx(fc, f ,logfc)) time. Since the leaders are now merged into the array L , the 
allocation incurs only 0(1) time. 

In Step 8, we prepare the sorting instances. Each leader x, for 1 < i < k, if 
Si(x) ± 0, then tags element A,[J?,(x) + l] with the address of the location 5 r [A,_i(x) + 
1]; note that A,[Jf?i(x) + 1] is to be copied to 5 x [A,_i(x) + 1]. Assign 0(f) processors to 
each segment as in Step 7; and for each element find the largest tagged element smaller 
than it in its array. As the second element in the segment is guaranteed to be tagged, 
the first being a leader, this does not require any inter-segment information access; a 
prefix minima over each segment, done in parallel, will do. Thus, each element can 
compute the address of the cell it must be copied to. Finally, the copying itself can be 
done in 0(f) time. So, the total time taken is 0(f + Prefix(& 2 , f ,21og&)). 

In Step 9, we again assign 0(f) processors per segment, as in Step 7. For 
each leader x, dedicate to x, a fraction 0(- p^-) of the processors assigned to s t (x), 
the segment of A,- to which x belongs, by carrying out another prefix summation over 
each segment of A,-, for 1 < i < k. That is each leader x is now assigned Q(M£l£) 
processors. Since the prefix summation can be clubbed with that of Step 8, this step 
can be finished in 0(SoRT(fc 3 , 2f-)) time. 

On a PRAM that allows concurrent reads, SoRT(m, 5 ), Prefix(to, 5 , b) and 
Merge(to, q ) are all at most 0(j- + log to). Hence, for a constant c, 

SoRT(TO c+1 ,qm c ) = 0 (Sort(to,§)) 


PREFlx(m c+1 , qm c ,b) = 0(Prefix(to,5, b )) 


and 


MERGE(m c+1 ,qTO c ) = 0(MERGE(m,g)). 
Thus, we have the following lemma: 


Lemma 2.2 When concurrent reads are allowed, k sorted arrays of a total size n, can 
be merged in 

0 fs0RT (k, f£) + PREFIX(£, £,log n) + MERGE(n,p) + 
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time using p processors. 


On a CREW PRAM, we can sort n elements in SoRT(n,p) = 0(- ^ -—+log n ) 
time (Upper Bound 2.9), find the prefix sum of n 6-bit numbers in Prefi x(n,p, 6) = 
0(j + logn) time (Upper Bound 2.7), and merge two arrays of a total size n in 
MERGE(n,p) = 0(z + loglogn) time (Upper Bound 2.1). 

Thus, we have, when p < n, 


SORT(fc, ^) = 0 (~j^“ + ‘ogi) = 0 + logfc) 

PREFIx(k, £,logn) = 0 ^L^ + log kj = 0 (j + log/:) 


Applying Lemma 2.2, we have the following theorem: 


( 2 . 1 ) 

( 2 . 2 ) 


Theorem 2.1 On a PRAM that allows concurrent reads, using p processors, k sorted 
arrays of a total size n can be merged in time 

fc_WAY_MERGE(n,p) = 0 



Theorem 2.2 On a CREW PRAM, using p > n processors, k sorted arrays of a total 
size n can be merged in time 


fc_WAY_MERGE(n,p) = 0 



logn 

log(p/n) 


+ logfc 


Proof: Use Upper Bound 2.9 for sorting, Upper Bound 2.2 for merging, and Upper 
Bound 2.7 for prefix summation, in Lemma 2.2. a 


Corollary 2.1 On a CREW PRAM, 

&_WAY_MERGE(7i,p) = 0(SORT(fc, *£) -f MERGE(n,p)). 


Theorem 2.3 On a CRCW PRAM, k sorted arrays of a total size n can be merged in 

°( l 0 g (i^)) + S0RT(i: '" ) ) 


time, using p > n processors. 
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Proof: On a CRCW PRAM, with p > n processors, we can find the prefix sum of n 
6-bit numbers in 


PREFIx(n,j>,») = 0 (min {log log »} + 

time (Upper Bound 2.6) and merge two arrays of a total size n in Merge(u,p) = 
O(log i 0 g°p" n ) ) time (Upper Bound 2.2). If, in addition, we use the best known CRCW 
PRAM sorting algorithm, by Lemma 2.2, the algorithm schema can be implemented in 
time 


0 (log 


log n 
log (pjn) 


+ min 


{ l0g i^kr l0 «*} + 


log A: 


= 0 ^log : 


logn 


+ 


log A; 


log {p/ n ) loglog(pk/n) 

Since SoRT(A,f) = !l( EiI ^ !w ), the theorem follows. 


log Iog(p/tfc) 

+ SOKT(t,^)). 


Sort (t,£)) 


□ 


Corollary 2.2 On a CRCW PRAM, k sorted arrays of a total size n, can be merged 
in A;_WAY_MERGE(n, p) 

■ 0 (‘ OS ( w) + + w&) 


Jog (p/n)) log (p/n) 

time, using p> n processors, 


Proof: Apply Upper Bound 2.10. □ 

Corollary 2.3 On a CRCW PRAM, 

A:_WAY_MERGE(n,p) = 0 (Sort(A:, *£) + MERGE(rc,p)). 


2.3.2 The Schema on an EREW PRAM 

Theorem 2.4 On an EREW PRAM, with p processors, k sorted arrays of a total size 
n can be merged in time 

( Tl lo£ k \ 

Ar_WAY_MERGE(n,p) = 0 (^logn-l — j 

Proof: We consider an implementation of the algorithm schema of Section 2.1. Again, 
we take A = fc 2 . Each step of the schema can be implemented on an EREW PRAM 
with q = Ll^rJ processors, as follows: 
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Step 1 : Since in this step, k > n 1 / 3 , log A: = 0(logn) and hence, Upper Bound 2.9 
can be used to sort the input elements 0 ( 2 J^Si _|_ l 0 gn) = O(logn) time. 

Step 2 : We can both partition the input arrays and form the leader arrays in 0(1) 
time, using n processors once they are assigned one per element. We allocate the 
processors uniformly over the input using Lemma 2 . 1 . This would take 0 (logn + 2 .) 
time. Each processor can sequentially select leaders from the elements assigned to it. 
So, the total time taken for this step is 0(logn + 2 .) = 0 (log nj. 

Step 3: Make k copies of each leader array; that is, a total of kt < 2 ±— < ^ entries 
have to be made. With ©( ) processors per leader, the copies can be made in 0(log k) 

time. Distribute processors per leader as in Step 2 and finish copying in O(logn) 

time. The total number of processors used is at most < q. For 

merging a pair of leader arrays A\ and A'-, for 1 < i,j < k, i ^ j, use the j-th copy of A(- 
and the i-th copy of A'-. Since, on an EREW PRAM with p processors two arrays of a 
total size n can be merged in 0 (logn+ 2 ) time (see Upper Bound 2.5), here with [ 53 — J 
processors per pair, the time is O(logrc). Since, there are me rge instances, the 

total number of processors used is at most 2 f^ n < jjj2_ < q. Note that the processor 
allocation incurs only 0(1) time in the case of merging. Thus, the total time required 
for this step is 0 (log n ). 

Step 4: Once the k ranks for each leader have been found, add the cross-ranks 
together in O(logn) time with L{^J processors per leader (see Upper Bound 2.7). 

Step 5 : Positioning the leaders in their corresponding segments has the difficulty 
that several leaders may want to search the same segment simultaneously. To avoid 
concurrent reads, make k copies of X, the merged array of all leaders. Let these copies 
be, U for 1 < i < k. With [j^J processors per element of X copying can be finished 
in 0(log n) time. Merge X; with A;. Assign processors for the i- th merge. Note 

that this processor distribution may require O(logn) time, as in Step 2. The merges 
and hence the cross- ranks can be found in O(logn) time (see Upper Bound 2.5). The 
total number of operations used is Yli=i + ««) = 0{tk -I- n) = 0 (n). 

Step 6 , 7 : Now that for all x 6 L, all R,(x)’s have been found, computing I(x) involves 
prefix sums computation, and takes at most O(logn) time using Upper Bound 2.7. 
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Computation of A,-(x)’s for all l £ L also requires prefix sums computation, and simi- 
larly, can be done in O(logn) time. 

Step 8: We prepare the sorting instances in this step. Each leader x, for 1 < i < k, 
knows 6{(x), the number of elements to be copied from A,- as well as the base address 
5 r [A,_i(x) + 1], where the smallest of those elements must go. If we have one processor 
per element, that is, n processors overall, each element can be informed the address 
of the location to which it must be copied in O(logfc) time using a broadcast over a 
distance of at most 6i(x) < k 2 and subsequently, the copying itself can be done in 0(1) 
time. With q processors this step would hence take 0( - = O(logn) time. 

Step 9: As n elements can be sorted on an EREW PRAM in 0(logn + nl ° gn ) time 
(see Upper Bound 2.9), with processors assigned to x, we can sort S(x) for each 

x 6 L in parallel in Q( AMl£gAM ) = O(logn) time. The processor allocation here 
would require a prefix summation over L that takes O(logn) time. 

That is, with q = [ n ^° s n fc J processors, k- way merging can be done in O(logn) 
time on an EREW PRAM. Hence the theorem. □ 

Corollary 2.4 On an EREW PRAM, 

fc_WAY_MERGE(n,p) = 0(SORT(fc, ^) + MERGE(n,p)). 

From Corollaries 2.1, 2.3, 2.4, 

Corollary 2.5 On a PRAM, 

fc_WAY_MlERGE(n,p) = 0(SORT(fc, *£) + MERGE(n,p)). 

2.3.3 Lower Bounds 

In this section we prove lower bounds on time for comparison-based k- way 
merging. Our main results are Theorem 2.5 for EREW PRAM, Theorem 2.6 for CREW 
PRAM, and Theorem 2.7 for CRCW PRAM. 

Lemma 2.3 On all models at least as powerful as EREW PRAM, 


ifc_WAY_MERGE(n,p) = ft(MERGE(n,p)). 
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Proof: Suppose that A and B are two sorted arrays of size n/2 each. Divide both 
A and B into k/2 segments of a non-zero size each. Here, by a segment we mean a 
subarray of consecutive locations; each of the segments hence is a sorted array. These 
k segments together form an instance of k- way merging of size n , a solution to which 
is clearly equivalent to the merge of A and B. Hence the lemma; □ 

Lemma 2.4 On all models at least as powerful as EREW PRAM, 

MERGE(ra,p) = fl(SEARCH(n,p)). 

Proof: Given an instance of search, A[ 1 . . . n] with search key x, and p processors, 
for two arrays B and C of size f + 1 each, for i = 1, ...,§, let B[i] = A[i] and 
Gfi+l] = -A[f + *]• If x > A[j\ then let JB[| + 1] = x,C[l] = -oo, else if x < A[f] then 
let B[ ^ + 1] = +oo,C[l] = x. Here — oo and +oo are distinct keys respectively smaller 
and larger than all other keys. Both of B and C being sorted, merging them will answer 
the search. Thus, Merge(u + 2 ,p) = fi(SEARCH(n,p)), and hence the lemma. □ 

Corollary 2.6 On all models at least as powerful as EREW PRAM, 

fc_WAY_MERGE(n,p) = ft(SEARCH(n,p)). 

Theorem 2.5 On an EREW PRAM, the time required for comparison-based merging 
of k sorted arrays of a total size n, using p> n processors, is 

fc_WAY_MERGE(n,p) = fl(logn). 

Proof: Snir [93] has shown that SEARCH(n,p) = ft (log n), on an EREW PRAM for 
any value of p. Thus the theorem follows from Corollary 2.6. O 

Lemma 2.5 On a CREW PRAM, even when the keys are single bit numbers, 

fc_WAY_MERGE(n,p) = ft(logfc). 

Proof: Given an instance (xj, . . . ,Xk) OR of size k, construct sorted arrays Zi, for 
1 < i < A:, of size n/k each as follows: set Zi[^\ to x,- and all of Zi[j] for 1 < j < 
to zero. Find the merge of all Zfs. If the last number of the output is one, then the 
OR is one, else it is zero. Cook et.al. [23] have shown that fl(log k) is a lower bound 
for finding the OR of fc-bits on a CREW PRAM. 1=1 
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Theorem 2.6 On a CREW PRAM, the time required for comparison-based merging 
of k sorted arrays of a total size n, using p> n processors, is 

&_Way_Merge(ti,p) = Q flog ( t—^j ^ ■+• log k 

V Uog (p/n)J 

Proof: Merging two sorted arrays of size n/ 2 each with p > n processors requires at 
least 

“('•■aft) 

time on the parallel comparison model [55, 14], and hence on a CREW PRAM, because 
the latter is strictly weaker than the former. Hence, by Lemma 2.3, 

fc_WAY_MERGE(n,p) = 0 . flog , *°. gre . 

V log (p/n)J 

For the second term, we use Lemma 2.5. □ 

Lemma 2.6 The number of comparisons required for merging k sorted arrays of a total 
size n, in t comparison rounds on the parallel comparison model is Qft^nk 1 ^ — n)). 

Proof: Consider the problem V of independently sorting [_n/fcj arrays of k elements 
each. Let be the i-th sorting instance, for 1 < i < [n/k\. Since, an ad- 

versary can always make sure that each element of the i-th instance is smaller than 
every element of the ( i + l)-th instance, the minimum number of comparisons required 
for solving V is asymptotically the same as [n/kj times the minimum number of com- 
parisons required for sorting one instance. Azar and Vishkin [9] have shown that the 
minimum number of comparisons required to sort k elements in t comparison rounds 
is Ct(t(k 1+1 ^ — k )). 

Now, construct a 2-dimensional array S of size [n/k\xk and let S[i,j] = (i,x*-). 
Comparing its elements lexicographically, S has the property that each element of the 
i-th row is smaller than every element of the [i + l)-th row. Hence, V can be solved by 
k- way merging all columns of S. The number of comparisons between the composite 
keys of S, that such an algorithm uses, should be at least as much as the minimum 
number of comparisons between the original keys that are required to solve V. Hence 
the lemma. a 

Lemma 2.7 The number of processors required for merging k sorted arrays of a total 
size n, in t units of time on the CRCW PRAM model is 

n( i"j 2 ^/') 




Chapter 2. Multiway Merging in Parallel 


27 


Proof: Consider the problem “P, dealt in the proof of the above lemma, again. The 
minimum number of processors required for solving V is asymptotically the same as 
[n/fcj times the minimum number of processors required for sorting one instance. Beam 
and Hastad [10] show that sorting fc keys in 0(t) time on a CRCW PRAM requires 
£l(2 ckl/C ) processors, for a constant c. □ 


Theorem 2.7 On a CRCW PRAM, the time required for comparison-based merging 
of k sorted arrays of a total size n, using p> n processors, is 


fc_WAY_MERGE(n,p) = Cl 



+ 


logfc 


log(l + p/n) 


+ 


logfc \ 

loglog(pfc/n)y ' 


Proof: The first term can be proved as in Theorem 2.6. Equating fc_WAY _MERGE(n, p) 
with t, the second and third terms follow from Lemma 2.6 and Lemma 2.7 respectively. 

□ 


Corollary 2.7 On a PRAM, 

fc_WAY_MERGE(n,p) = ft(SORT(fc, &) -f MERGE(n,p)). 

Proof: On a PRAM that does not allow concurrent writes, SoRT(fc, = ©(logfc). 
Hence, from Lemma 2.4, Theorem 2.5 and Theorem 2.6, the corollary can be shown to 
hold for EREW and CREW PRAMs. From Theorem 2.7, the CRCW PRAM model 
also satisfies the result. n 

From Corollaries 2.5 and 2.7, we have the following theorem: 

Theorem 2.8 On a PRAM, 

fc_WAY_MERGE(n,p) = 0(SORT(fc, $) + MERGE(n,p)). 

2.4 Multiway Merging of Integers 

2.4.1 An Upper Bound from the Schema 

On an ARBITRARY CRCW PRAM, we can sort n integers in SoRT(n,p) = 
0(logn+ nlo S > lo S n ) time (Upper Bound 2.11). Thus, 

SoRT ( fc, ft = 0 (*™ + log :*) = 0 (=!»]=** + log *:) (2.3) 

Using this, together with, Eq. (2.2) and Upper Bound 2.1, we get: 
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Lemma 2.8 On an ARBITRARY CRCW PRAM , when the array elements are given 
to be integers, using p processors, k sorted arrays of a total size n can be merged in 
time 

fc_WAY_MERGE(n,p) = 0 (\ogk + log log n + 

However, if the integers are from the range [0 ... m - 1], we may use the algorithm of 
Upper Bound 2.12 to sort integers in SoRT(n,p) = O( lo 1 ° 1 ^ n + loglogm -f nlo 8*°g TO ) 
time. Thus, 


SORT(fc,2i) 


0 

0 


Ik loglogm log k 

V pk/n loglogk 

In loglog m logfc 

V p log log k 


+ loglogm 
+ loglogm 


5 


For prefix summation, we can use the faster algorithm of Upper Bound 2.6. Hence, 
PREP,x(*, g .log n) = O (mi. { lo8i -^_,logt} + 

- 0 {****'+ 


Combining these with Upper Bound 2.1, from Lemma 2.2, we have the following lemma: 


Lemma 2.9 On an ARBITRARY CRCW PRAM, when the array elements are inte- 
gers drawn from the range [0 ... m — 1], using p processors, k sorted arrays of a total 
size n can be merged in time 


0 


f logfc 
V log log & 


+ log log mn + 


nlog log m 

V 


Alternatively, we can use the MERGE(n,p) = 0 (log loglogm + j*) time algo- 
rithm for merging integers in the range [0 ... m — 1] (Upper Bound 2.3). Together, with 
Eq. (2.3) and Eq. (2.2), this gives, from Lemma 2.2: 


Lemma 2.10 On an ARBITRARY CRCW PRAM, when the array elements are as- 
sumed to be integers drawn from the range [0 ... m — 1], using p processors, k sorted 
integer arrays of a total size n can be merged in time 

n log log k \ 


O ^log k + log log log m + 


V 


Combining Lemmas 2.8, 2.9, and 2.10, we have: 
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Theorem 2.9 On an ARBITRARY CRCW PRAM , when the array elements are as- 
sumed to he integers drawn from the range [0 . . . m — 1], using p processors , k sorted 
integer arrays of a total size n can be merged in time fc-WAY-MERGE(n,p) = 

0 (-j?£p + log log(min{n, log m}) + min {log* + 1^2^, log log™, + 

2.4.2 Lower Bounds 

Lemma 2.11 On a CRCW PRAM, even when the keys are single bit numbers, for 
p>n,. 

fc_WAY_MERGE(n,p) = ( ~T— 1 • 

V log log py 

Proof: Given an instance (®j, . . . , xCj of Parity of size k, construct sorted arrays Z{, 
for 1 < i < k, of size n/k each as follows: set Z,[|] to x, and all of Z,jJ] for 1 < j < J, 
to zero. Find the merge of all Z? s. If the last zero of the output is at index m, then 
the number of ones in this input is n — m. Since the lower bound for finding the parity 
of k bits with p processors on a CRCW PRAM is [10], the lemma follows. 

□ 

Given an instance of k- way merging of n single bits numbers, the solution 
can be obtained by adding together the number of zeroes in each array. This takes 
PREFlx(fc,p,log7i) time, which, by Upper Bound 2.6, is G( l0 gf 0 ~ )- That is, the lower 
bound of the above lemma is tight when the integers are of single bits. 

Theorem 2.10 On a CRCW PRAM, with p>n processors, the time required to merge 
k sorted arrays of integers from the range [0 ... ^ - 1] of a total size n is ^( i og iog ( pfc/n j)' 

Proof: Similar to the problem V dealt in the proof of Lemma 2.6, consider a problem 
Q of independently sorting arrays (x\,. . . ,z].), of single bit numbers, for 1 <i< [n/ fcj . 
Sorting each row independently, with p processors overall, will require 
time (see Lemma 2.7). 

Replace each x) by y) = 2 i + a;}. Thus y) are from the range [0 ... ^ - 1]. 
Construct a two dimensional array S with [n/fcj rows and k columns and let S[i,j] = yj. 
Clearly, S has the property that each element of the t-th row is smaller than every 
element of the (i + l)-th row. Since each column of S' is a sorted array, Q can be solved 
by k- way merging all columns of S. 

Hence the theorem. D 



Chapter 2. Multiway Merging in Parallel 


30 


2.4.3 Two More Upper Bounds 

Integers in the range [0 ... m — 1] can be k- way merged as follows. 

Let LesSj [i] be the the total number of elements that are smaller than integer 
i in Aj, the j-th input array, for 0 < i < m and 1 < j < fc; Lessj[i\ can be computed 
by merging the array (0,...,m - 1) with Aj. Compute Si = Y,j=i Lessj[i], for all 
i, in parallel. Let Count j[i] be the the number of occurrences of integer i in A 3 ; if 
Firstj[i] and Lastj[i] are the indices of, respectively, the first and last occurrences 
of i in Aj, then Countj[i\ = Lastj[i] — First j [j] -f 1. For each integer i, compute 
PS.Count[i,j] = ^(=i Count\[i ]. Now, copy the h-th occurrence of i in Aj to the 
location numbered Si + PS.Count[i,j - 1] + h of the output. 

With p processors, the merging in the computation of Less can be done in 
0(a(m) + y) time, by Upper bound 2.4. For obtaining First and Last, each input key 
has only to compare itself with the adjacent keys in its array; that is, with n processors, 
Count can be computed in 0(1) time. With p processors, hence, the computation of 
Count , as well as the final copying can be done in 0( n ~ ( p mfc ) time. That is, the total 
time taken is 0( a ^- +a(m) + PR EFI x(fc, ^,logn)). Using, Upper Bound 2.6 for prefix 
summation, hence, we have the following theorem. 

Theorem 2.11 On a CRCW PRAM, with p processors, k sorted arrays of integers 
from the range [0 ... j- — 1], of a total size n, can be merged in 

fc_WAY_MERGE(n,p) = 0 (i + «(»/*) + log i^. /n) + Egj^j) 

time. 

When the output does not have to be given in an array, we can do better, as the 
following theorem indicates. 

Theorem 2.12 Given k sorted arrays of integers from the range [0 . . . n— 1], an ordered 
chain of the integers can be obtained in 0( lo gfcg - fc + a(n)) time with L a ( n )+l!)g 5 i/fogi 0 gTJ 
processors on an ARBITRARY CRCW PRAM. 

Proof: Let p = L o( n )+iI^Jfc7 fog log I J • Follow the algorithm schema of Section 2.1, with 
the following minor variations. Choose A = k 2 . 
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Step 4: Use the cross ranks computed in Step 3 to find for each leader x, its largest 
smaller x t - in A,-; that is, x x = A([rj(x)] when x £ A[ and X{ = A[[ri(x) — 1] when x £ A' x . 
Find the largest among all x t -’s. We now have an ordered chain of all leaders. 

Merging two arrays once they are cross ranked takes only 0(1) time with linear 
number of processors even on an EREW PRAM [55]. With A; 2 processors assigned to 
x, the largest among all x t -’s also can be found in 0(1) time (see Upper Bound 2.13). 
Hence, with processors per leader this step takes at most 0(£) time. 

Step 6: Do not perform the addition. 

Step 9: For each leader x, padded-sort S x into Z x . Ordered-chain Z x \ this coupled 
with the ordered chain of all leaders, obtained in Step 4, gives an ordered-chaining of 
all input elements. 

By Upper Bound 2.14, padded- sorting can be done in O(j^^p-) time with 
0(n log log A;) operations overall. Ordered chaining Z x can be done in O(log* A;) time 
[79, 83] with 0(n ) operations overall, by replacing every non-zero element by its index 
in Z x . Coupling the local ordered chains with the global one can be done in 0(1) time 
with 0(n) operations. 

Hence, the total time taken is, 

0 ( iogiogfc + PRBFix(fc, ft, log k ) + PREFix(A,p,logn) + Merger, f|) + n * os i ^ - fc -^ 


Using Upper Bound 2.4 for merging and Upper Bound 2.6 for prefix summation, we 
have the theorem. a 
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Chapter 3 

Parallel Merge Sort 

Given an array of n elements to sort, we can proceed as follows: divide the 
array into 6 segments of equal size, sort each segment recursively, and <5 -way merge the 
sorted segments. Following this framework, and using the multiway merging schema 
of Chapter 2, a parallel merge sorting algorithm is obtained. This algorithm runs in 
. 20(log*(l°gn/loga)) s t e p S on the parallel comparison model. 

Valiant, who introduced the parallel comparison model, also gave, on that 
model, a simple 0(log(|^|^)) time algorithm for merging two arrays of a total size n, 
with ncx processors, 2 < a < n [96]. Borodin and Hopcroft [14] proved a matching 
lower bound. For sorting n elements with na processors on the parallel comparison 
model, Azar and Vishkin proved a lower bound of [9]. But, unlike merging,, 

no simple sorting algorithm matching the lower bound, for all values of a, has been 
found. Alon indeed obtained an algorithm [9] that matches the lower bound, but being 
an adaptation of the AKS sorting network [3, 81, 82] of Ajtai, Komlos, and Szemeredi, 
this algorithm is not simple and has large constant factors. Subsequently, Cole [19] 
gave a simple merge sort algorithm with small constant factors; this algorithm runs in 
0(^1 • log log a) steps on the parallel comparison model. Note that except for very 
small values of a, our algorithm is closer to the lower bound than Cole’s merge sort; 
moreover it is very different from Cole’s merge sort. In our algorithm too, vis-a-vis 
Alon’s algorithm, the constant factors involved are small. 

On a CREW PRAM, Cole’s merge sort is optimal— with n processors it can 
be implemented in O(logn) time; on a CREW PRAM, with any number of processors, 
S7(logn) is a lower bound on sorting [23]. In contrast, the proposed algorithm, when 
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implemented on a CREW PRAM with n processors, runs in logn .2°( lo s* n ) time; this 
implementation may be thought simpler than Cole’s merge sort because it does not 
use an intricate pipelining scheme. On a Concurrent Read Concurrent Write (CRCW) 
PRAM This algorithm runs in . 2 °( lo s (logn/ log a)) s t e p S on the parallel comparison 
model. 


3.1 A Procedure 

In this section we describe a procedure Segregate. This procedure is the 
same as our k - way merging schema, except that the last step of the schema is not in- 
cluded here. For the implementation details not covered here, please refer to Chapter 2. 

In a procedure call “Segregate(A, m, 6)”, A is an array of m elements with 
6 sorted subarrays A,-, 1 < i < 6, of size m/6 each, given one after another contiguously. 
The output of the call is a permutation of the input elements, given in the array A 
itself, that would have m/6 (not necessarily sorted) subarrays Gj(A) of size gj(A), 
gj(A) < £ 2 , 1 < j < m/6, given one after another contiguously, so that all elements of 
a subarray are larger than every element of a preceding subarray. 

The procedure is described below. The number of processors available is m/3 , 
/3 > 0. 

Procedure SEGREGATE(A,m,6) 

For 1 < i < 6, let every 6-th element from A,- be called a leader; the leaders 
from Ai, while given in the same order as they appear in A;, form, let us say, a leader 
array A(. We call the set of elements in A,- between any two consecutive of its leaders, 
a segment. There are m/6 2 leaders in each A,-. In parallel, merge every pair of leader 
arrays by assigning mfi/6 2 processors per pair. Now, for 1 < i,j < 6, each element of 
A'i knows its rank in A'-. With /36 processors per leader, in parallel for each leader x, 
add together the ranks of x; all the leaders can now be arranged in sorted order, say, 
in an array, L. 

Each leader x now knows the segment s,(x) of A; to which it belongs. For 
each x € X and for each array A,-, search s,-(x) in parallel to find the two consecutive 
elements in it that straddle x; assign (3 processors per leader per search. Thus, each 
x € X has now found its rank Ri(x) in every A,-. For each x £ X, find the prefix 
sums over X! t -(x), 1 < i < 6; x now knows its index in the output array. Also, for 
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every consecutive x,y € L, we can now tag the first element (where one exists) from 
the interval (x,y) in A{ with its index in the output array. Now. a prefix minima 
computation over each segment is enough to inform each element of .4 of its index in 
the output array. Permute the elements in A according to the indices found. 

The following lemma is now apparent: 

Lemma 3.1 The procedure call “SEGREGATE(A,m,<5)” takes 

0(PREFix(£,<5 2 /3,logm) + MERGE(m,m/3) + Search(£, 3)) 
time with m/3 processors. 


3.2 The Algorithm on a CREW PRAM 

First we present the CREW PRAM implementation of the algorithm. We 
construct a sequence of sorting algorithms, the (log* n — 2)-th one of the sequence 
achieving the claimed resource bounds, using the following idea: 

Given an 0(c r lognlog^ r ' n) time sorting algorithm A, r > 2, for = 

0(1), an 0(c r+ i lognlog( r+1 ) ra) time sorting algorithm B can be obtained. 

A high level description of the new algorithm B would be as follows: 

Divide the input array into subarrays of size 6 2 each, for a parameter 6 to 
be chosen. Sort each subarray using algorithm A. Now it remains only to 
merge the Jj- sorted subarrays. Make groups of 6 of these subarrays and 
apply the procedure Segregate on all groups in parallel. The elements of 
each group now appear in a sequence of (not necessarily sorted) subarrays 
of size at most 6 2 each, so that all elements of a subarray are larger than 
every element of a preceding subarray; sort these subarrays in parallel using 
algorithm A, and each group would be sorted. That is, now the input 
elements are arranged in j? sorted subarrays of size S 3 each. Repeat this 
process of combining S sorted subarrays at a time until a single sorted array 
is left. 

For an appropriate choice of S, this scheme gives the desired time bound. 

We now give the construction scheme more formally. The basic building 
block of the scheme is Valiant’s algorithm that merges two arrays of a total size n 
in O(loglogn) time with n processors on a CREW PRAM [96, 55]; let a be the con- 
stant hidden by the O-notation here. Guided by a complete binary tree on n leaves, 
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Procedure CREW_SSORT h.(I,n) 
begin 

8 «- (\/Iogn) 1 / los(h+2) n ; 

for 1 < i < n/8 2 in parallel do /* If = J[(i - 1 )8 2 + 1, . . . , i£ 2 ] */ 
CREW.SSORT 
for j *— 1 to (log n/ log 8) — 2 do 
begin 

for 1 < i < n/ 6 i +2 in parallel do 
begin 

Segregate(I/,6 j+2 ,<5); 
for 1 < k < 8i +1 in parallel do 

CREW.SSORT h -i(G k (lj),g k (lj)y, 

end 

end 

end 


proceeding bottom up, and using Valiant’s algorithm at each level, the standard merge 
sort thus runs in a log n log log n time on an n processor CREW PRAM [55]. This is the 
first algorithm in our sequence. We denote this algorithm by CREW_SSORTo(/, n). 

For h > 0, CREW_SSORT/i(I,n) is given above. Let If be I[(i - 1)P +2 + 
1, . . . , i£ J ' +2 ], the i-th subarray of size 6 J+2 in J; 0 < j < (log n/ log 6) — 2, and 1 < i < 
n/8 j+2 . Note that I(f2i)s+i’ • • • > 1*6 1 constitute If. 

On a CREW PRAM, PREFix(n,n,&) is O(logn), Merger, n) is O(loglogn) 
[55], and Search(<5, 1) is O(log<5). Thus, by Lemma 3.1, with one processor per input 
element, the procedure call “Segregate^/, 8 j+2 , 8)” takes 0 (log 8 + log log n) time 
on a CREW PRAM. If S h (n ) is the time taken by CREW_SSORT fc (J,n) with n 
processors, then 

Sh(n) < 0(1) + S h -i(S 2 ) + - 2^ (O(logtf + a log log n) + S h -i{8 2 )) 
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Since on a CREW PRAM, with any number of processors, 0(log6) is a lower bound 
on sorting <5 2 elements [23], log 6 = 0(Sh- i(<5 2 )). Thus, for an appropriate constant b, 

Sh ^ n) - ab ' • 1 °g 1 °S 71 + b ■ • Sh-^6 2 ) 

We claim that, 

„ , w (2b) h (4b - 1) - 2 b , , , A+2 > 

Sh(n) < 2b -1 a log n log* A+2 l n 

The proof is through induction. When h = 0, CREW_SSORTo(/, n) runs in 
a log n log log 7i time; this is the basis. Suppose the claim is true for h = H — 1. Then, 


S H (n) 


< 

< 


log 6 
logTl 

log< 


•log log n + b ■ 
• log log n + b ■ 


logn 
log <5 
logn 
log <5 


Sh-i(8 2 ) 

(26)^-^ - 1) - g i W|<2 

26—1 


Since log 6 = 


21og("+ 2 >n’ 


S*W < 2ab ■ ( 1 + ( 2t ) g "( 46 1L - Jt ] log„lo g («+2) „ 


26 — 1 


- wm^t^. alognlog mv n 


Thus we have the following theorem: 


Theorem 3.1 The procedure CREW_SSORTi og * n _ 2 (f, n) sorts an array I of size n, 
with n processors, on a CREW PRAM, in logn • 2°l log "1 time. 

On a CREW PRAM, PREFlx(n, ^,6) = O(logn), Merge^k^) = 
O(loglogn) [55], and Search(< 5, 1) = 0(log<5). The standard optimal merge sort 
runs in O (logn log logn) time on an i og fc, gn processor CREW PRAM [55]. Thus, by 
Lemma 3.1, the procedure call "Segregate^/, <5 J+2 , 6)” takes 0 (log <5 + log log n+ £) 
time with p < n processors on a CREW PRAM. Hence, we have the following theorem: 

Theorem 3.2 The procedure CREW_SSORT r (I,n) sorts an array I of size n, with 
— ^ 2 y n processors, on a CREW PRAM, in 0(logn -logl r+2 l n) time, for a constant 
integer r > 0. 

Once again, assume that n processors are available; then the entire sequence 
of the above algorithms can be captured in the following recursive structure. 
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Procedure CREWJRSORT(/,n,<5) 
begin 

if n < 6 2 or 6 = 0 then 
6 <— ylogn; 
if n = 2 then 

sort I with 2 processors in 0(1) time, exit; 

Divide the input I into subarrays I\, .... Is of size n/6 each; 
for 1 < i < 6 in parallel do /* recursively sort the subarrays */ 
CREWJRSORT(I,-, f,6); 

Segregate(J,ti,<5); 

for 1 < k < n/6 in parallel do 

CREW_RSORT(G fc (/),5jk(J),0); 

end 


The time taken by a call “CREW_RSORT(/,n,£)” is 

S(n,6 ) = 0(1) + S + ^(logtf log log n) 4* S(6^, 0) 

Thus, for an appropriate constant b, 

S(n,S ) < S +6(loglogn + 5(<5 2 ,0)) 

Unrolling the first term, 

S(n,6) < sQ^)+2&(loglogn + S(£ 2 ,0)) 

= S (j^,6^j +rb(\og\ogn + S(6 r ,0)), for r > 2 

= S{6 2 ,6) + — 2^ f>(loglogn + *S'(^ 2 ,0)), with r = (j2|| - 2 ) 

From the first statement of the procedure, 5(<5 2 , <5) = 5(^ 2 , 0). Thus, 

S(n,6) < 5'(^ 2 ,0) + 6^^y-2^(loglogn + ^ 2 ,0)) 
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To sort an array 7 of size n, call “CREW_RSORT(7,n,0)”. The time required, with 
n processors, is 


5(n,0) = S(n , \Aogra) < 26 log n + 26 • — ■ - - ■ 5(logn,0) 

log log n 

Unrolling this recurrence relation, we get, for r > 1, 


£(n,0) 



logn + 


(26) r log n 

log( r + 1 ) n 


S(\og^ n, 0) 


2°( r > logn 

logl r+1 ) n 


■ 5(log^ r ) n, 0) 


Thus, J(n, 0) = 2°( lo s* n ) logn; the procedure call “CREW_RSORT(7, n, 0)” is equiv- 
alent to the procedure call “CREW_SSORTi og * n _ 2 (7,n)”. 

Moreover, since log^ n elements can be sorted in 0(log^ r+1 ^ n ) time with 
(log (r) n) 2 processors on a CREW PRAM [77], we have the following theorem: 


Theorem 3.3 The procedure call “CREW_RSORT(7, n, 0)”, with the recursion ter- 
minated after r = 0(1) levels, sorts an array I of size n in O(logn) time with n log^ 7 "^ n 
processors on a CREW PRAM. 


3.3 The Algorithm on the Parallel Comparison Model 


A generalisation of the CREW PRAM algorithm, given in the previous sec- 
tion, gives the recursive procedure CMP_RSORT(7, n,6) for the parallel comparison 
model; the generalisation is only in the choice of 6. To sort an array 7 of size n, call 
“CMPJRSORT(7, n,0)”. We shall show, later in this section, that the time taken by 
this call, with na processors, is • 2°( lo s*( lo s n / lo S“)) 

On the parallel comparison model, since MERGE(n, na) = 0(log(|2|i)), and 
SEARCH(n,p) = 0(|3|^),by Lemma 3.1, the procedure call “Segregate( 7, n,6)” takes 
°( lo s(fe) + §1) time with na Processors. 

Following the analysis of CREW_RSORT(7,n, 6) in the previous section, 
since is the lower bound on sorting 6 elements with 6a processors, the time 

taken by a procedure call “CMP_RSORT(7,n,6)” is, with na processors, for an ap- 
propriate constant 6, 


S(n,6) 


< S 

< b • 


(?'‘) +61og GS) +6 ' M 
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Procedure CMP_RSORT(J, n,6) 
begin 

if n < 6 2 or 6 = 0 then 

6 4- \/ Q lo s( lo s n / logo); 

if n = a then 

sort I with a 2 processors in one step, exit; 

Divide the input I into subarrays ij, . . . , Ig of size n/6 each; 
for 1 < i < 6 in parallel do /* recursively sort the subarrays */ 
CMP_RSORT(/j,n/£, 6); 

Segregate^, n,^); 

for 1 < k < n/S in parallel do 

CMP_RSORT(G*(/), <?*(/), 0); 

end 


To sort an array I of size n, call “CMP_RSORT(I,n,0)”. The time required, with n 
processors, is 


S(n,0) = S(n,\/a lo 8( lo s' l /logo')) 

< 2 b ■ l — + 2b • logn/logp . 5( a log(logn/loga) ) Q) 

_ logo log(logn/loga) 


Unrolling this recurrence relation, we get, for r > 1, 


S(n, 0) 



logn (26) r logn/loga 
logo log( r ^(logn/loga) 


5( a l°6 (r) ( lo gn/log^) )0 ) 


With r = log*(logn/loga), 5(n, 0) = 2°( lo s*( lo s n / lo s a )) - ^ • 5(a ? 0). On the parallel 
comparison model, ot elements can be sorted in unit time with o? processors: com- 
pare every element with every other element. That is, S^OjO) = 1* Thus, *?(n,0) = 

9O (log* (log n/ log a)) . logn 
log a* 


Theorem 3.4 The procedure call “GMP_RSORT(/,n,0)" sorts an array I of size n 
in 2 0 ( lo s’( lo s n / lo s“)) • time with not processors on the parallel comparison model. 
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3.4 The Algorithm on a CRCW PRAM 


The CRCW PRAM implementation of the algorithm is similar to the parallel 
comparison model version, except in that when n = a, with ct~ processors, sorting 
would require 0(,jjS£|_) time instead of 0(1) time. On a CRCW PRAM, 


Prefix(»,p,6) = 0 + min (log j^.logn}) 

[91, 41], Merger, not) = O(log(g|£)), and Search(tz,p) = O(jgJ). Thus, by 
Lemma 3.1, the procedure call “Segregate(/, n y 6)” takes 0 ( ^logSa + jf||+log( ) ) 
time with not processors on a CRCW PRAM. Since, +■ jj^^) is the lower 

bound on sorting 8 elements with 8 a processors, it is clear that the recurrence relation 


5 („, 0 ) = 26 .!^ + 26 . r -!^^ 

log a log(log n/loga) 


5 ( a lo g( Io s n / log “) 5 o) 


holds for CRCW PRAM too, but with a larger value of b than for the comparison model 
version. Solving this, we have the following result: 


Theorem 3.5 The procedure call “CMP_RSORT(/, ra,!))” for sorting an array I of 
size n can be implemented in 2° - los (log 71 / logo)) . jJgjLR— time with not processors on a 
CRCW PRAM. 


For the CRCW PRAM model, Beame and Hastad have shown an fi( i 0 ~ °o g n ~ a ) 
lower bound on sorting with na processors [10]. Coupled with the comparison model 
lower bound, thus, sorting requires Q + logiogno ) ti me on a CRCW PRAM. But 
it is not known whether this lower bound is tight. Cole’s merge sort runs on a CRCW 
PRAM in C( i 0 g &g time. For a while this was believed to be the best possible. But 
off late, entirely different techniques [45, 35] have suggested faster algorithms, though 
still not matching the known lower bound. Our algorithm, though slower, could be 
deemed simpler than these algorithms. 
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Chapter 4 

Optimal Sublogarithmic Time 
3-Colouring of Rooted Forests 


Being a bipartite graph, a rooted forest is 2-colourable. But it would take 
flt l&Iw ) ^ me to 2-colour even a linked list, even with as many as a polynomial number 
of processors [32]; besides, there is evidence to suggest that this lower bound could be 
strengthened to fl(logn) [90]. On the other hand, a rooted forest can be 3-coloured 
in sublogarithmic time: Goldberg et.al. [32] and Han [48] describe an 0(log* n) and 
an O(log(log* n)) time algorithm respectively — both algorithms run on an n processor 
CREW PRAM; these algorithms are not optimal. Attempts to find an optimal parallel 
sublogarithmic time algorithm for 3-colouring have been successful only for rooted 
forests with a constant maximum degree [84, 87]. The algorithm of [84] takes 
time on a CRCW PRAM; the time bound was improved to <9(log^ n) on an EREW 
PRAM in [87]. In this chapter, we present an 0 ((log log rc) log* (log* n)) time optimal 
algorithm on the TOLERANT CRCW PRAM, for 3-colouring general rooted forests. 

The fractional independent set (FIS henceforth) problem is to find an inde- 
pendent set that constitutes at least a constant fraction of the vertices. We show that 
given an f{n ) time algorithm for 3-colouring a rooted forest F, an FIS of F can be 
found in 0(/(n)) time with the same number of processors, on both COMMON and 
TOLERANT CRCW PRAMs; in other words, the FIS problem is no harder than 3- 
colouring. Furthermore, with randomisation, an FIS of a rooted forest can be found 
in 0(1) time on an n processor COMMON or TOLERANT CRCW PRAM. We use 
this result to obtain an optimal algorithm for 3-colouring a rooted forest with high 
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probability, in 0 ( (log* n) log* (log* n)) time. 

On rooted forests, top-down algebraic computation problems like 2-colouring, 
depth computation for all vertices, and prefix summation can be solved using pointer 
jumping in O(logn) time with n processors on a CREW PRAM. If the input forest is 
either known to be a linked list [6], or given in adjacency list representation [1], then 
all these problems can be solved in O(logn) time optimally on an EREW PRAM. On 
the other hand, when each vertex knows only its parent, the known logarithmic time 
optimal algorithm [49] can at best be implemented on a COLLISION CRCW PRAM. 
The 3-colouring algorithm proposed here suggests optimal O(log n) time algorithms 
for these problems on a much weaker model — TOLERANT. On a COMMON CRCW 
PRAM, we show, these problems can be solved in 0(log n) time with a time-processor 
product of 0(n log(log* n)). 

To obtain the above algorithms on TOLERANT, we show that a TOLERANT 
PRAM of size N with a linear address space, can be slowed down by any factor A = 
0(log log iV), with no asymptotic increase in space or cost. This is a non-trivial result 
since TOLERANT is not known to be self-simulating in general [39]. Also, this result 
complements a corresponding observation in [38] that concerns factors A = 0(1). 

4.1 Preliminaries 

A PRAM model is said to be self-simulating , if for all N > n > 1, a PRAM of 
that model of size n can simulate a single step of another PRAM of the same model of 
size N in O(^) time. All CRCW PRAMs that are at least as powerful as COLLISION 
or COMMON are self-simulating. TOLERANT is not known to be self-simulating [39] 
(see also [4]). 

The approximate prefix summation (APS) problem over a sequence (xj, . . . ,x n ) 
of non- negative integers is to find a a sequence (0 = yo, yi,. . . ,y n ) so that for every 
1 < i < n, yt = OQCj-i xj) and y,- - y,-_i > s,-. The linear approximate compaction 
(LAC) problem over an array of n cells, at most m of which contain an item, the others 
being empty, is to arrange all the items into an array of size 0 (m). In ordered LAC, 
the output should have the same relative ordering of items as the input. Clearly, APS 
can be used to solve ordered LAC. Goldberg and Zwick [35] show that APS (and hence 
ordered LAC) can be solved in 0 (log log n) time using t — j*— processors on a COM- 
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MON CRCW PRAM; a matching lower bound is also known [15]. Using the 0(1) time 
n 2 processor simulation of COMMON on TOLERANT [66], it is easy to see that their 
algorithm can be extended to TOLERANT. Thus, with p < processors, on a 

TOLERANT PRAM, APS can be solved in 0(j + log log n) time: divide the instance 
into p groups of equal size, find the prefix sums in each sequentially, and finally combine 
the sums using a global APS. 

We define “an inverse” H of the logstar function as follows: H(0) = 1 and 
H(i) = for i > 1; if H[x) = n then log* n = x. 


4.2 Designing Optimal Algorithms 

We say, a problem V is e-contractable, 0 < e < 1, if for every instance X of 
V of size n , there exists an instance Y of V of size en , such that Y can be constructed 
from X , and a solution to X can be found from a solution to Y. 

Theorem 4.1 (cf. [40, 30, 31]) Suppose linear approximate compaction (LAC) of size 
n has an 4- C(n)) time algorithm with p processors. Let V be a problem that can 
be solved in T(n) time with n processors, and can be contracted in 0 ( 2 - 4 - f2(n)) time 
with p processors. If the model used is self-simulating, then V can also be solved in 
0(2 + -R(n)log T(n) + 0(n)log* T(n) + T(n)) time with p processors. 

Proof: To prove the lemma, it is enough to show that V can be j^Ly-contracted 
in + jR(n)logT(n) + C(n)log* T(n)) time with p processors, since, we can then 
solve the contracted instance in 0(T(n)) time with processors using the given 
algorithm. Starting with the given instance of size n, proceed in phases. In the i- 
th phase we perform at least an ^lf(^(0)/5'^ +1 ^(0)j -contraction followed by an LAC 
of the array containing the contracted instance. Thus, the size of the array at the 
beginning of phase i is at most A,- = n/H^\ 0). 

The contraction in the i-th phase proceeds in H^\ 0) stages, the j-th stage of 
which performs H (j — 1) consecutive ^-contractions followed by an LAC of the array 
containing the contracted instance. That is, the j-th stage performs a ^yyy- contraction. 
Clearly, the size of the array at the beginning of stage j is at most 2?,,y = A,/ n^Io R (*0- 
Hence, the size of the array after .H'M(O) stages is is at most A,/n^o^ H(k) < 
Ai/H(' +1 \ 0) < A, + 1 . The time taken by stage j is 0(H(j - 1)(^ J - + R(n)) 4- C(n)). 
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Summing the time taken, over all stages, the i-th phase takes 

0 + R(n)H{H®( 0) - 1) + C{n)H®( 0)) 

= 0 (^ + + C(n)\og' 

time. Note that H(H^( 0) — 1) = log.ff( t+1 )(0). After the last phase, say the r-th one, 
the array size is at most Hence, H^ T+1 \0) = T(n). So, summing the time taken 
over all phases, we get the claimed time bound. Hence the lemma. □ 

TOLERANT, unlike other standard PRAM models [39], is not known to be 
self-simulating. Thus, Theorem 4.1 cannot be applied on TOLERANT directly. But, 
we show that when the memory used is linear in the number of processors, TOLERANT 
is self-simulating. That is, when both the given linear processor algorithm, and the |- 
contraction routine can be executed on TOLERANT with linear space, Theorem 4.1 can 
be applied on TOLERANT as well. The following theorem states the self-simulating 
property more formally: 

Theorem 4.2 One step of a TOLERANT CROW PRAM of size N that has an address 
space of size O(N) can be simulated on another TOLERANT CRCW PRAM of size n 
in 0 ( ) time where n < log ^ gJ v > without any asymptotic increase in memory. 

For a discussion on earlier results, as well as a proof of this theorem see 
Section 4.5. 

4.3 Optimal 3-colouring of rooted-forests in sublogarith- 
mic time 

Our algorithm is based on the observation that in every rooted forest F = 
(V, E) at least half the vertices have at most one child; let Z be the set of these vertices. 
On a TOLERANT CRCW PRAM, the vertices in Z can be identified in 0(1) time, 
with 2n processors: If a node v and all its children simultaneously write in v and the 
write succeeds, then v is a leaf. On the other hand, if all children of v simultaneously 
write in v and the write succeeds, then v has exactly one child. Remove the vertices of 
Z from the forest jF[V]; the resulting forest is F[V — Z], 

Since F[Z], the forest induced by the removed vertices, is a collection of linked 
lists, we can find a 3-colouring <r : Z —* {1,2,3} of F[Z] in 0(^ + log'^ n) time with p 
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processors [48]. Using the colouring <r, the given 3-colouring <f> : (V - Z) -*■ {1,2,3} of 
the remaining forest F[V - Z] can be extended to the original forest F[V]: 

for each x 6 {1,2,3}, in turn, 

for each removed vertex z with o(z ) = x, 

set (j>(z) to the minimum colour not present in 
z' s neighbourhood in the original forest F[V]. 

Since each removed vertex z has at most two neighbours in the original forest, such a 
colour can always be found from {1,2,3}. 

Thus the problem of 3-colouring a rooted forest can be ^-contracted in 0(~ + 
log^ n ) time using p processors on a TOLERANT CRCW PRAM. Note that the 
address space used is of size 0(n ) only. Hence, from Theorem 4.1 and Theorem 4.2, we 
have the following result: 

Theorem 4.3 A rooted forest can be optimally Z-coloured in 0((log log n) log* (log* n)) 
time on a TOLERANT CRCW PRAM. 

Corollary 4.1 A fractional independent set of a rooted forest can be found in 0( l ^°f Q r ^ n ) 
time optimally on a TOLERANT CRCW PRAM. 

Proof: Given a 3-colouring of a rooted forest F, a three-fold counting of the vertices 
can give the colour that is used the largest number of times; this colour is used by at 
least one third of the vertices. Use the time optimal prefix sums algorithm 

[22] for counting. a 

With the help of randomness, both the FIS and 3-colouring problems on rooted 
forests can be solved much faster. 

Lemma 4.1 A fractional independent set of a rooted forest can be found with high 
probability, in 0(1) time on a COMMON or TOLERANT CRCW PRAM, with n 
processors. 

Proof: Let a set I be initialised to the set of leaves in F = (V,E), the given forest. 
Each vertex in T informs its parent that it is in I using a concurrent write. For each 
vertex v of F, “toss” a fair coin; let the result of the toss t(v) be 0 or 1 randomly with 
equal probability. For every vertex v which is not a root, if v does not have a child in 
I, and t(v) = 1, and t(p(v)) = 0, then add v to I. 
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Let Vj be the set of vertices with exactly one child in F and say, | Vi |= n x . 
That is, V — Vj consists of the leaves as well as the vertices with two or more children. 
So, at least a constant fraction of V — Vi are leaves, and hence, if n\ — o(n ), X is an 
FIS. On the other hand when n x = 0(n), at least a constant fraction of the vertices in 
Vi have been added to I, with high probability [55]. Hence the lemma. □ 

Since, every vertex in Z has at most one child in F[V], a 3-colouring of F[V - 
Z ] can be easily extended to F[V], As there is an optimal randomized Las Vegas 
algorithm for LAC that runs in 0(log* n) time with high probability [31, 37], applying 
Theorem 4.1, we have: 

Theorem 4.4 There is an optimal randomized 0 ( (log* n) log* (log* n )) time algorithm 
for 3-colonring a rooted forest on a CRCW PRAM. 

In the rest of this section, we show how to use the sublogarithmic time optimal 
3-colouring algorithm to get <3 (log n) time optimal algorithms for decomposable top- 
down algebraic computation problems (henceforth referred to as DTAC) on rooted 
forests, like 2-colouring, depth computation, and prefix summation. 

An instance of size n of the top-down algebraic computation problem [1] can 
be represented by a triplet (S,F,F), where S is a set, F is a set of functions from S 
to S', and F = (V, E) is a forest on n vertices in which each vertex v has a unique label 
l v such that, l v G S when v is a root and /„ € F otherwise. The problem is to compute 
a function S : V — *■ S, which can be recursively defined as follows: for v G V, 6(v) = l v 
if v is a root, and S(u) = l v (6(p(v))) otherwise. 

A top-down algebraic computation is said to be decomposable [1] if (i) each 
f € IF is computable, (ii) given an algorithm N(f) that computes / 6 F, and s € S, 
f(s) can be computed in 0(1) time sequentially, (iii) for each / 6 F, each vertex of 
the input forest that is labelled by / can compute N(f) in 0(1) time, (iv) 7 is closed 
under composition, and (v) for /i,/ 2 € F, 1V(/ X o / 2 ) can be computed in 0(1) time 
sequentially. 

When S = {0, 1} and F contains only the Boolean negation operation, DTAC 
reduces to 2-colouring. In Ranking or depth computation S is the set of integers, the 
roots in the input forest are labelled by 1, and F contains only the increment function. 
When S is the set of integers and F = {/,• | fi(x) = x + i; x,i 6 S'}, DTAC reduces to 
prefix summation. 
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Given an instance (S,F , F) of DTAC of size n, where F is a rooted forest, 
find an FIS Z of F using Corollary 4.1. For each z £ Z, remove z from F, and for each 
child w of z, set p(w) = p(z) and l w — l w o l z . Let us call the resultant rooted forest 
F' . By the definition of DTAC, once Z has been found, the construction of F' can be 
done in 0(1) time with n processors, without concurrent writes. 

Assume that the solution S' : (V — Z) — ► S for (S,F,F') has been found. The 
solution for the original DTAC S : V — ► S can be obtained from S' by the following 
computation: S(v) = l v if v £ Z and v is a root in F; 6(v) = l v (6(p(v))) if v € Z and v 
is not a root in F ; and 6(v) — S'(v) otherwise. By the definition of DTAC, again, this 
computation would take 0(1) time with n processors, without concurrent writes. 

Thus we have that DTAC can be —contracted in 0( lo 1 ° 1 ^ > ” n ) time with 0(n ) 
space and 0(n) operations on a TOLERANT CRCW PRAM. 

Hence, from Theorem 4.1 and Theorem 4.2 we have the following result: 

Theorem 4.5 An instance of DTAC of size n can be solved on a rooted forest in 
O(logn) time with processors, on a TOLERANT CRCW PRAM. 

4.4 A deterministic algorithm for the fractional indepen- 
dent set 

In this section we show that given an f(n) time algorithm for 3-colouring a 
rooted forest F = (’ V,E ), an FIS of F can be found in 0(/(n)) time with the same 
number of processors, on both COMMON and TOLERANT CRCW PRAMs; in other 
words, the FIS problem is no harder than 3-colouring. Instead of counting the number 
of vertices of each colour to find the largest colour class in a 3-colouring of F, we 
recolour some of the vertices so as to ensure that at least one-fourth of the vertices are 
of colour 1. This is done by guaranteeing that each vertex not of colour 1 has a near 
enough descendant of colour 1. The method is described in more detail below: 

Let L be the set of leaves in F; remove all leaves from F after giving them 
colour 1 (they are never to be recoloured). Let Fl = F[V — i] be the resultant forest. 
Obtain a 3-colouring of Fl. 

Every root of colour 1 (2 or 3) assumes that it has a parent, grand parent and 
great grand parent of colours 2, 1 and 2 (1, 2 and 1) respectively. Note that we are 
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not really adding any new vertex to F L . Recolour each vertex in F L with the colour of 
its great-grand parent, real or assumed. Now, every vertex is coloured the same as its 
siblings. 

Starting at a child of a vertex of colour 3, and going towards the root, the 
possible prefixes (of size three) of colour sequences encountered are 1-3-2, 2-3-1, 1-3-1 
and 2-3-2. We now affect a recolouring so that only instances of 1-3-2 remain in Fi. 
First, sequences 1-3-1 and 2-3-2 are removed by recolouring the vertex of colour 3 (the 
middle one) with colour 2 and colour 1 respectively. Now, the only possible extensions, 
terminating at a colour 2 vertex, for a 2-3-1 sequence are 2-3- 1-3-2 and 2-3- 1-2. Recolour 
these into sequences 2-1-2- 1-2 and 2-1-3-2 respectively. Note that recolouring these 
sequences does not require a concurrent write, if each vertex remembers the very first 
colouring received by itself and a few of its immediate ancestors. Every vertex of colour 
3 now has a parent of colour 2 and children of colour 1. A vertex of colour 2 that does 
not have a child of colour 1, has all its children coloured 3, and hence surely has a 
grandchild of colour 1. So, every vertex of Fi is either coloured 1 or has a child or a 
grand-child coloured 1. 

Reinsert L into Fl~ If there are colour conflicts, then recolour the parent in 
each conflict with colour 4. Every vertex of F is now either coloured 1 or has a child 
or a grand-child or a great-grand child of colour 1. At least one fourth of the vertices 
in F are now coloured 1. 

The power of concurrent write has been used only to identify the leaves as 
well as to locate colour conflicts. Hence, we have the following theorem: 

Theorem 4.6 If a rooted forest F can be Z-coloured in f (n) time on a CRCW PRAM 
model, then a fractional independent set of F can be found in 0(f(n)) time on the same 
model, provided the model is powerful enough to compute the OR of n bits in 0(1) time 
with 0{n ) operations. 

Thus, using the algorithm of [48] for 3-colouring, DTAC can be “contracted 
in O(log(log* n)) time with n processors on a COMMON CRCW PRAM. 

Hence, using Theorem 4.1 with minor modifications, we have the following 

theorem: 

Theorem 4.7 An instance of DTAC of size n, can be solved on a rooted forest in 
O(logn) time with 0{n log(log* n)) operations on a COMMON CRCW PRAM. 
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4.5 The Self-Simulation Property of TOLERANT CRCW 
PRAM 

In this section., we present the proof of Theorem 4.2, which essentially states 
that a TOLERANT PRAM of size N and memory 0(N ) can be slowed down, by a 
factor of A = fi (log log N), without any cost or space penalties. A corresponding ran- 
domised result is known: Gil, Matias and Vishkin [31] have shown that a TOLERANT 
PRAM of size N and memory 0(s) can be slowed down, with high probability, by 
a factor of A = fl(log* N) without a cost penalty but with an 
ditive overhead in space. Note that our result is limited to TOLERANT PRAMs of 
linear space and to A = fi(loglog N), but takes only 0(N) space. On the other hand, 
Grolmusz and Ragde [38] have shown the self-simulation property of a linear space 
TOLERANT PRAM for A = 0(1). 

Let ‘Tjv” denote “a TOLERANT PRAM of N processors,” in the sequel. 

Lemma 4.2 One step of P/y that has an address space of size s can be simulated on 
T n in 0( ~) time, when s,n < N and n = 2s. 

Proof: If every processor of IV that is in a conflict were to back out from writing, 
then the remaining writes can be self- simulated as in a CREW PRAM, in O(^) time. 
We use a flag for each memory cell to signify conflicts. The flag of a cell would take 
on values 0, 1 or greater than 1, depending on whether 0, 1 or more processors of IV 
want to write in that cell. Initially, all flags are set to 0. 

Divide T n into two equal halves: the guards and the actors. Each guard stands 
by a unique memory cell. Also, divide T jv into g = [~~] groups of size j each. Activate 
the groups one by one, in g rounds, using the actors; the i-th actor deputes for the 
i-th processor activated in the current round and is associated with the cell that this 
processor wants to write into. In each round, the actors and guards together can be 
viewed as forming a rooted forest, where each guard is a root and the actors associated 
with it are its children. Whether a root (guard) has 0, 1 or more children (actors), can 
be ascertained in 0(1) time on a TOLERANT PRAM: if a node v and all its children 
simultaneously write in v and the write succeeds, then v has no child; on the other 
hand, if all its children simultaneously write in v and the write succeeds, then v has 
exactly one child. 
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For each cell, in every round the guard increments its flag by 0, 1, or 2 de- 
pending on whether it has 0, 1, or more children. It is easy to see that at the end of 
g rounds the flags would correctly reflect the conflict status. The processors can now 
read the flags and back out if necessary. 

Hence the lemma. □ 

Corollary 4.2 One step of r_y that has an address space of size s can be simulated 
on T n in O(yf-) time, where n < s,N. 

Proof: Divide the memory into [2r] segments of size n each. Proceed in rounds, 
in the i-th round simulating the writes into the i-th segment using Lemma 4.2 in O(y) 
time. □ 

The problem of padded-sorting n values is to output them in sorted order in 
an array of size at most 0(n)\ unused cells of the output array are to he filled with 
some special character [44]. 

Lemma 4.3 Given m integers in the range [1 • • • r] in an array of size m, they can 
be padded-sorted in 0(f) time with p processors, on a TOLERANT CRCW PRAM, 
where r = 0(f) and p < log £ gm - 

Proof: We employ a method similar to the one used for bucket-sorting m integers in 
the range [1 • --log 0 ^ m] in [22]. The following procedure padded-sorts m integers in 
the range [1 • • • r] with the claimed resource bounds. 

Step 1: Divide the set of m integers into p groups of size \f] each; the last group 
may be of a smaller size. Assign one processor for each group and sort it sequentially. 
For a p x r array C, let C[i,j] be the number of items with key j in group i. Let the 
first C[i,j] locations of an array contain the items with key j in group i. 

Step 2: Let D be an array of size pr. Let D[(j — l)p + *’] = C[i,j], for 1 < i < p and 
1 < j < r; that is, we copy C into D column by column. Hence, D[(j-l)p+l], . . . , D\jp] 
contain the number of items with key j in groups 1, . . . ,p respectively. 

Step 3: Find the approximate prefix sums over D. For 0 < i < pr, let <*(*) be the 
i-th sum. For 1 < i < p and 1 < j < r, copy the first C[i,j] locations of A iy j into 
consecutive cells of a base array of size a(pr), starting at address o((j — l)p-f-i — 1) + 1. 
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Time taken is clearly 0(f). The objects in the base array, where there are, now appear 
in sorted order. p 

Lemma 4.4 One step of Tjv that has an address space of size N can be simulated on 
Fn in 0(f ) time without any asymptotic increase in memory, where n < i og jV ■ 

Proof: Let Pi, ... , Pjv be the processors of Ty. 

Divide the address space into |"^] segments of size n each. For 1 < / < N , 
let S(l) be the number of the segment into which P( wants to write. Form an array of 
ordered pairs (S(l),l) and padded-sort these on the first key, using Lemma 4.3; with n 
processors, this takes 0(f) time. Thus, all processors of I\y that wish to write in the 
same segment have now come together in a base array of size 0(N). 

Consider the segments one by one. Suppose the processors attempting to write 
in the i-th segment are spanning over IV,- cells of the base array. If N, > n, applying 
Corollary 4.2, the write attempts into the i-th segment can be simulated in 0(JV,/n) 
time. When iV,- < n , with n processors the simulation clearly can be done in 0(1) time. 

The total time taken for the simulation is 0(T,f-) = 0(f). Note that the 
space used outside of the approximate prefix summation routine is linear. Hence the 
lemma. □ 

We can now prove Theorem 4.2, as follows: Let us assume that the address 
space is of size cN , for some constant c. Divide the memory into segments of size N 
each. Proceed in c phases; in the i-th phase, use Lemma 4.4 to simulate the processors 
writing in the i-th segment of the memory. Each phase takes 0(f) time; so the total 
time taken is 0( £ f) = 0(f). Hence the theorem. 

In this proof, when the deterministic algorithm for approximate prefix sum- 
mation is replaced by the 0(log*n) time randomised algorithm [31], we get a corre- 
sponding randomised. result for n < This is an improvement over the earlier 

self-simulation of linear space TOLERANT [31] in terms of memory used. 
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Chapter 5 

Brooks’ Colouring 


In this chapter, an 0(A 2 log Alogn) time algorithm for A-colouring a Brooks’ 
graph, with n/logn processors, on an EREW PRAM, is obtained. In particular, this 
algorithm 3-colours a 3-degree Brooks’ graph in O(logn) time, with n/logn processors, 
on an EREW PRAM. (The basic idea of this algorithm, in the context of a CREW 
PRAM, was treated in the author’s M. Tech, thesis [85].) See also [88] This is an 
improvement over the previous best known algorithms, on bounded degree graphs [58, 
80]. (See the discussion in Section 1.2.3.) 

We also obtain the following combinatorial result: for any two vertices u and 
v that are more than a constant distance apart in a 3-degree Brooks’ graph G, there 
exist two distinct 3-colourings of G of which one has u and v coloured the same, and 
the other has u and v coloured differently. Thus, no vertex in a Brooks’ graph can force 
the colour of vertices arbitrarily far away. This highly local nature of Brooks’ colouring 
can be seen as suggesting that on a CRCW PRAM, the problem may have a o(log n) 
time algorithm. 

For the problem of A-colouring, a colour c is said to be feasible at v if v does 
not have a neighbour coloured c. An uncoloured vertex in a partially coloured graph 
is said to be at impasse if it has no feasible colour in {1,...,A}. If v is a vertex at 
impasse, then its neighbour of colour i will be denoted by Vi(i = 1, . . . , A). 

An a-P component in £ is a component of the subgraph induced by vertices 
coloured a or (i. Note that interchanging colours <x and (3 in an ct-fi component does 
not affect the validity of the colouring. So, for a vertex v at impasse if v a and vp 
belong to different ce-/3 components interchanging colours in one of them will resolve 
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the impasse. 

We define a construct called fork as follows: given a partially coloured Brooks’ 
graph G = (V, E), a fork at v € V is a minimal simple path F of at least two vertices 
with v as an endpoint, so that the other endpoint e(F) of F either (i) has a degree of 
at most A — 1 in G, or (ii) has a degree of A in G, and has two neighbours, both of 
which are not in F, and are coloured the same. If v is at impasse, and the neighbour 
of v in F is the a-coloured neighbour v a of v in G, then we call F an a- fork. 

For each vertex w in F, w ^ e(F), let the successor of u?, be the farther from 
v of its neighbours in F. 

If v is at impasse, then using F the impasse can be resolved as follows: Re- 
colour each x € F, x ^ e(F), with the colour of its successor in F. Uncolour e(F). 
Give e(F) one of the colours missing in its neighbourhood; note that, by definition, 
such a colour should exist. 

Similarly, if v is coloured, then using F we can recolour v: Recolour each 
x € F, x e(F), with the colour of its successor in F. Uncolour e(F). Give e(F) one 
of the colours missing in its neighbourhood. 

5.1 An Algorithm for Brooks’ Colouring 

First, we present a procedure that will be subsequently used in the algorithm. 
Procedure MAXIMAL _SET(A, a,p) 

Input: A A-regular graph G = (V, F), in which every a-(3 component is a simple path 
and for every vertex v at impasse v a and vg are end points of ct-j3 chains. 

Output: A maximal set of a- (3 components in G. such that no two members in the 
set touch the same impasse vertex; a component F touches a vertex v if at least one 
vertex in T is adjacent to v. 

Step 1: Form a graph H in which each vertex corresponds to an a-p component 
of G. Place an edge between two vertices of H if and only if the a-fi components 
corresponding to them touch the same impasse vertex. Since an a-P component can 
touch at most two impasse vertices, one at each end, the maximum vertex degree of H 
is 2. Remove all isolated vertices from H. 
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Step 2: Find the connected components of H and identify each as a chain or a 
cycle. From each cycle remove a vertex. Let the resulting graph be H'. List rank each 
component of H' . For every odd ranked vertex of H' interchange colours a and (3 in 
the corresponding a-/3 component of G. 


An algorithm for 3-colouring of 3-degree graphs is now presented. 

Procedure 3J2olourJ3ubicJjraphs 

Input: A cubic graph G — (V,E) in adjacency list representation. 

Output: A 3-colouring a : V — *• {1,2,3} of G. 

Step 1: Find an MIS M of G, and for each v G M let <r(t>) = 3. The graph G - M 
has a maximum degree of 2; that is, G — M consists of disjoint chains and cycles. From 
every odd cycle of G - M elect a representative into a set I. Clearly, G - M - I is 
2-colourable; 2-colour it. Each vertex of J is at impasse, and is left uncoloured. 
Remark: For each vertex v € V, we keep two flags P(u) and Q(v). Each flag can take 
on values “up” and “down”. Initially all flags are down. Every flag that is up, will be 
associated with an impasse vertex; that is, every flag that is up will hold a pointer to 
an impasse vertex. As soon as the impasse is resolved at a vertex, all its associated 
flags will be brought down. We shall use the following procedure extensively. 

Procedure RESOLVE 
begin 

For u G J do in parallel 

if v has at least one colour missing in its neighbourhood then 
give v the minimum feasible colour, remove it from J, 
and bring down every flag associated with v 

end. 

In the sequel, a vertex with its P-flag up will be called a P- vertex, and a component 
in the subgraph induced by P- vertices a P- component; and, similarly for flag Q. 

Step 2: For each odd cycle C of G - M, associate the P-flags of all coloured vertices 
in C with the impasse vertex in C, and put up all of them. 

Remark: Every odd cycle of G — M provides a path between v\ and V 2 , if v is the 
impasse vertex contained in it. Moreover, each vertex on this path is coloured 1 or 2 
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and its neighbour outside the path is coloured 3. Thus, t'j and uj belong to the same 

1 2 component, which is a simple path with ui and v 2 as its end points. 

Step 3: Recolour all colour-1 vertices with colour 2 where feasible. Recolour all 
colour-3 vertices with colour 2 where feasible. Call RESOLVE. 

Remark: Every 1 — 3 component of G is a now simple path. Thus, each of V 3 and Vi 
is an end point of a 1 — 3 chain. That is every impasse vertex has exactly two (not 
necessarily distinct) 1—3 chains going out of it. 

Step 4: Call MAXIMAL _SET(3, 1, 3). Call RESOLVE. 

Remark: Impasse is resolved for a vertex, if colour was changed in any one of the two 
1-3 components emanating out of it. For each v € I, now we have a simple 1-3 path 
in G with Vi and V 3 as its end points. But, note that now it is not necessary for Vi and 
V 2 to be in the same 1—2 component, let alone a 1-2 path. 

Step 5: For each impasse vertex v, associate all <2 -flags in the 1-3 path from Vi to 
V 3 with v, and put up all of them. 

Step 6: Recolour all colour-3 vertices with colour 1 where feasible. Recolour all colour- 

2 vertices with colour 1 where feasible. Call RESOLVE. Call MAXIMAL_SET(3,3,2). 
Call RESOLVE. 

Remark: For each v € /, now we have a simple 2-3 path in G with and as its 
end points. But, v\ and V 2 may not be in the same 1 — 2 component. And v\ and V 3 
may not be in the same 1 — 3 component. All P-components and ^-components of G 
are simple chains, by construction. Also, no vertex can have both its flags up, unless 
it is a common endpoint of a P-chain and a Q-chain. So, every impasse vertex v has a 
P- path Lp(v ) and a Q-path Lq(v ) of its own. 

Step 7: For each impasse vertex v, if Lp(v ) has at least one vertex with a non-P- 
neighbour of colour 1 or 2, find either a 1-fork or a 2-fork at v that involves only vertices 
from Lp(v). Use this fork to resolve the impasse at v. Call RESOLVE. 

Remark: When Lp(y ) has at least one vertex with a non-P-neighbour of colour 1 or 2, 
the two endpoints of Lp(v) being coloured 1 and 2, a fork at v can be found from Lp{y) 
in at least one direction. If v is still at impasse after this step, Lp(v ) is a maximal 1—2 


component of G. 
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Step 8: For each impasse vertex v , if Q p{ v ) has at least one vertex with a non- 
neighbour of colour 1 or 3, find either a 1-fork or a 3-fork at v that involves only 
vertices from Lq(v). Use this fork to resolve the impasse at v. 

Remark: So, we are left with only the case where u, and Vj are the end points of a 
simple i-j path for 1 < i < j < 3. 

Step 9: For v € I, let = x = y and V 3 — z. If y is adjacent to both x and 2 , 

* 

then x and z are not adjacent because G has no 4-cliques, and y is the only neighbour 
of colour 2 for both x and 2 ; let cr(x) = < 7 ( 2 ) = 2, cr(y) := 1, and cr(v) := 3. If y is 
adjacent to at most one of x and 2 , then interchange colours 1 and 3 in Lq(v). Now, 
clearly, Lp(v ) can provide a fork at v in at least one direction. Use this fork to resolve 
the impasse at v. Call RESOLVE. 


Theorem 5.1 The procedure 3-Colour-Cubic-Graphs 3 -colours a f- clique free cubic 
graph in O(logra) time and O(n) operations on an EREW PRAM. 

Proof: Correctness of the procedure is obvious from the remarks following individual 
steps. 

In Step 1 of the procedure, we have to find an MIS of G. We do this by first 
finding a (A-fl)-colouring C : V — ► {1,.. . , A + 1} of G using the optimal algorithm 
[87] for (A-fl)-colouring a bounded degree graph, and then, for i := 1 to A + 1, in turn 
adding w € V to M, if C(w) = i and no neighbour of w is already in M . The time 
taken is 0(A 2 log Alog^ n) (for any fixed k > 1) with - processors on an EREW 
PRAM. That is, with processors (k = 1), the time taken is 0(A 2 log Alogn). 
Besides, there are 0(1) instances of, all vertices having to scan their respective neigh- 
bourhoods in parallel; these can be solved in 0(A) time and 0(nA) operations on a 
CREW PRAM. The rest of the procedure is dominated by a constant number of in- 
vocations to the list ranking algorithm, and hence can be solved optimally in O(logn) 
time on an EREW PRAM [6]. Hence the claim on resource requirements holds for a 
CREW PRAM. 

For two adjacent vertices u and w, let [u.w] be the entry for the edge (u,w) 
in it’s edge list. We assume that [tt, w] and [to, u] have a pointer to each other. These 
pointers, if not given as a part of the input, can be easily created in 0(1) time using 
0(n 2 ) space and n processors on an EREW PRAM— form an adjacency matrix in 
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which each “non-zero” item is a pointer to an edge list entry; no initialisation is needed 
as we will never look at a “zero” . It is easy to see that with this structure available, all 
vertices can scan their neighbourhoods in parallel, in 0(A) time, without read conflicts. 
Hence, the algorithm can be run on an EREW PRAM also, with the same resource 
bounds. 

Note that Procedure RESOLVE can be implemented on an EREW PRAM 
in O(logn) time. Since, the P- vertices and the 0 -vertices associated with an impasse 
vertex form a list each, an invocation to list ranking is sufficient to bring down the flags 
in all of them. □ 

Since for any subcubic graph on n vertices, a cubic graph on 0(n) vertices 
of which the former is a subgraph, can be constructed in constant time with 0(n ) 
processors [58]. Thus, we have the following corollary: 

Corollary 5.1 A 4-clique free 3-degree graph can be 3-coloured in O(logn) time and 
0(n ) operations on an EREW PRAM. 

Now, an algorithm for Brooks’ colouring a general graph is presented. 

Procedure A-Colour- Regular-Graphs 

Input: A (A-fl)-clique-free regular graph G = (V,E) in adjacency list representation. 
Output: A A-colouring a : V -*■ {1, . . A} of G. 

Step 0: If A < 3 use the procedure 3-Colour-Cubic-Graphs. 

Step 1: Find an MIS M of G, and for each v € M let o(v) = A. The graph G-M 
has a maximum degree of (A— 1). From every A-clique of G — M elect a representative 
into a set I. As G - M - I does not contain any A-clique, it is (A - l)-colourable. 
Recursively (A - l)-colour G - M - I. Each vertex of I is at impasse, and is left 
uncoloured. For each v e I, select the minimum (3 e {1,...,A - 1} such that v A and 
are not adjacent. Swap colours between v x and vp. 

Remark: For each v € I, (?[u,^,...,»a-x] is a A-clique of G. Also, Vi and va are 
not adjacent. Observe that, v A is not adjacent to all of {ui, . . . ,u A _i} because, v A is 
adjacent to v and G does not contain a (A+l)-clique. 

Step 2: Recolour, first all colour-1 vertices, and then all colour-A vertices, with a 
colour other than 1 and A where feasible. Call RESOLVE. 
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Remark, Now, every 1 A component of G is a chain or a cycle. Also, for each v <E /, 
v\ and ua are end points of 1-A chains. 

Step 3: Call MAXIMAL_SET(A, 1, A). Call RESOLVE. 

Step 4: For each impasse vertex v, associate all P-flags in the 1-A path from tq to 
va with v, and put up all of them. 

Step 5: Repeat steps 2 and 3, with colour 2 substituted for colour 1. 

Remark: Now vi and v& are the end points of the same 2-A path, for every v € I. 
But v\ and va Reed not even be in the same 1— A component. All P-components of G 
are simple chains, by construction. So, every impasse vertex v has an exclusive P-path 
Lp(v). 

Step 6: For each impasse vertex v, if possible, find a fork at v that involves only 
vertices from Lp(v), and use this fork to diffuse the impasse at v. Call RESOLVE. 
Remark: If Lp(v ) does not provide a fork as required, then no vertex w of Lp(v) has 
two non-P-neighbours of the same colour. Also, both P-neighbours of w, when there 
are two, must be coloured the same. Since v\ and «a are the end points of Lp(v), this 
would mean that Lp(v) is a maximal 1-A component of G. That is, v\ and ua are the 
end points of the same 1— A path Si(v) = Lp(v), and V 2 and ua are the end points of 
the same 2 — A path S 2 (v). The only vertex common to S’i(u) and S 2 (v) is v&. Since 
v\ and «a are n °t adjacent, S\(v) does not degenerate into a single edge. But, 52 (w) 
may be a single edge. 

Step 7: For each impasse vertex v , interchange colours 2 and A in 5 > 2(v). No 
neighbour of v\ is now of colour 2, because, prior to this step, tq and va were not 
adjacent. Recolour V\ with 2, and the impasse at v is resolved. 


Theorem 5.2 The procedure A-Colour-Regular-Graphs A-colours a regular Brooks’ 
graph in 0 (A 2 log A log n) time with processors on an EREW PRAM. 

Proof: From the proof of Theorem 5.1 it follows that with j j—j processors on an EREW 
PRAM, all steps of Procedure A-Colour- Regular- Graphs, except the (A+l)-colouring 
and the recursive call in Step 1, can be implemented in 0 (A log n) time, and that the 
(A+ 1)- colouring can be done in (9 (A 2 log A log n) time. Note that G — Mis properly 
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coloured and does not have a vertex of colour A + 1. That is, the restriction of the 
colouring C to V - M is a valid A-colouring oiG-M. Hence, G-M -I is already 
A- coloured when it is passed on to the next lower level of recursion and need not be 
coloured again. In other words, we need only one invocation of the colouring algorithm, 
at the beginning of the procedure; the lower levels of recursion can make use of the 
same colouring. 

Hence the theorem. □ 

Corollary 5.2 A Brooks’ graph can be A-coloured in 0(A 2 log A log n) time with 
processors on an EREW PRAM. 

5.2 Local Nature of Brooks’ Colouring: A Combinatorial 
result 

We say that a pair of vertices u and v in a fc-vertex-colourable graph G are 
k-friends iff they have the same colour in every fc-colouring of G, and k-foes iff they 
have different colours in every colouring of G. For example, for a vertex v in a chain 
L, every vertex at an even distance from v in £ is a 2-friend of v, and every vertex at 
an odd distance from v in L is a 2-foe of v: in contrast, v does not have a 3-friend, and 
its neighbours are its only 3-foes. 

In the previous section, we showed that the problem of 3-colouring a 3-degree 
Brooks’ graph can be reduced to a constant number of invocations to list ranking. 
Hence it is intermediate in complexity to 2-colouring and 3-colouring of chains. 

The problem of 3-colouring a chain has an O(log(log* n)) time linear processor 
EREW PRAM algorithm [48]. In contrast, the problem of 2-colouring a chain has a 
lower bound of fI(logn/loglogn) on time, with a polynomial number of processors, on 
a CRCW PRAM [32]. This lower bound is proved using a reduction from Parity. Given 
an instance of Parity, its l’s are linked up in the same order in which they appear. The 
resulting list is then two-coloured. Comparing the colours of the first and last vertices, 
we get the solution of the instance of Parity. This reduction relies on the fact that, in 
a chain, vertices can have 2-friends and 2-foes arbitrarily far away; fixing the colour of 
one vertex fixes the colour of all others. 

Investigating the distribution of 3-friends and 3-foes in a cubic Brooks’ graph, 
we prove the following two theorems. The proofs are presented in Section 5.2.1 



Chapter 5. Brooks’ Colouring 


60 


Theorem 5.3 For every 4-clique-free 2-degree graph G , and for any two vertices u and 
v in G apart by more than 361 edges, there exists a 2-colouring of G that has u and v 
coloured the same. 

Theorem 5.4 For every 4-clique-free 2-degree graph G . and for any two vertices u and 
v in G apart by more than 362 edges, there exists a 2-colouring of G that has u and v 
coloured differently. 

Thus, in a 3-degree Brooks’ graph, vertices can have their 3-friends and 3- 
foes at most a constant distance away. Hence, it appears that there may be no 
o(logn/log log n ) time reduction from Parity to 3-colouring of 3-degree Brooks’ graphs. 
Thus, probably, for 3-colouring of 3-degree Brooks’ graphs, we can expect to find an 
algorithm that is asymptotically faster than the O(logre) time algorithm of the previ- 
ous section, unlike the 2-colouring of chains, where the O(logra) time optimal EREW 
PRAM algorithm of [6] is believed to be the best possible. 

Panconesi and Srinivasan [80] show that for a Brooks’ graph G and a vertex 
v in it, a A-colouring of G — v can be extended to the whole of G by recolouring a 
path of length 0(log A ra). They also prove that this bound is tight in that there exists 
a graphs G and a vertex v in it so that extending some A-colourings of G — v to the 
whole of G will require recolouring a path of length 0(log A n). Our theorems can be 
seen as complementing these results. Note that our conjecture is not weakened in any 
way by these results, because a Brooks’ colouring algorithm, on the one hand, need not 
rely on a partial initial colouring that is subsequently modified to cover the whole of 
the graph, and on the other, even when it does may avoid the degenerate case of [80]. 

Thus, the 3-friends and 3-foes of a vertex in a cubic Brooks’ graph cannot be 
more than a constant distance away. 

5.2.1 Proofs of the Theorems 

Theorem 5.3 is proved first; Theorem 5.4 will then be derived as a corollary. 

Let G = (V,E) be a 3-degree Brooks’ graph on n vertices, and let u and v be 
two of its vertices apart by more than 361 edges. Consider a 3-colouring o of G such 
that cr(y) / cr(u). We show that the colouring o can be modified into another valid 
colouring of G that has both u and v coloured the same. Let x, y and z be the three 
neighbours of v. 
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We use the following proposition, which will be proved in Section 3, in our 
proof of Theorem 5.3: 

Proposition V(u, v, x, y, z) 

If x and u are connected in G — {v, j/, 2 }, then the 3-colouring <7 of G can 
be modified into a 3-colouring o' of G - v so that 

• o(x) (?'(x) 

• c(y) = <r\y), a(z) = a\z), o(u) = a'(u) 

Continuing with the proof of Theorem 5.3, we assume that u and v are con- 
nected in G, because otherwise, a permutation of colours in the component of u (or of 
v ) is enough to give the same colour to both u and u. We consider the following three 
exclusive possibilities: 

1. Among the three neighbours of u, there is a cut- vertex of G (say 2 ), the removal 
of which disconnects the other two from u. That is, 2 and y are disconnected 
from u in G - z; in other words, every path from x or y to u in G contains 2 . 

2. Among the three neighbours of v, there is a pair (say, y and z), the removal of 
which disconnects the third from u. That is, x is disconnected from u in G— {y, z}; 
in other words, every path from x to u in G contains y or z. 

3. For each of the three neighbours of v, there is a path in G that connects it to u, 
but does not contain the other two. That is, x and u are connected in G — { 2 /, 2 }, 
and y and u are connected in G — { 2 , 2 }, and z and u are connected in G — { 2 , y }. 

Without loss of generality assume that cr(u) = 1 and o(y) = 2. 

Case Cl: 2 and y are disconnected from u in G — z 

Let us denote the component that contains u in G - z by A. 

Since every path from 2 or y to u in G contains 2 , there must be a path from 
z to u that does not contain 2 , y or v . Hence, z and u are connected in G — {v, 2 , y}. 

Since v is coloured 2, 2 , y and z are coloured 1 or 3. If z is coloured 1, then use 
Proposition 'Piu , v, 2 , 2 , y) to recolour z without affecting 2 , y or u ; otherwise, continue. 
Uncolour v. 

The colouring now satisfies the following conditions: v is uncoloured; u is 
coloured 1; 2 is coloured 1 or 3; y is coloured 1 or 3; z is coloured 2 or 3. 

Since v is uncoloured, one of the following four exclusive possibilities must hold. 
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(I) . Colour 1 is free at v. Here, v can be coloured 1, the same as u . 

(II) . Colour 2 is free at v. Here, z must be coloured 3. (Either z was initially 

coloured 1 and an invocation to Proposition 'P(u,v, z,x, y) has changed its colour to 
the present 3, or z was initially coloured 3 and no colour has been changed since then.) 
Swap colours 1 and 2 in component A; u gets colour 2, but x , y and z retain their 
colours. Now v can be coloured 2, the same as u. 

(III) . Colour 3 is free at v. Here, z must be coloured 2, and both x and y coloured 

1. Swap colours 1 and 3 in component A; u gets colour 3, but x, y and z retain their 

colours. Now v can be coloured 3, the same as u . 

(IV) . No colour is free at v. Since x and y are not coloured 2, z must be coloured 

2. Without loss of generality, assume that x and y are coloured 1 and 3 respectively. 

First assume that x and z are not in the same 1-2 component. Swap colours 
1 and 2 in the 1-2 component that contains x; x gets recoloured 2, while y, z and u 
retain their colours. Now, v can be coloured 1, the same as u. 

Next assume that x and z are in the same 1-2 component. So, if at all z has 
a neighbour of colour 3, it must be in A; the other two valencies of z are taken up by 
v and a vertex coloured 1. Since y is not in A, this means that y and z are not in the 
same 2-3 component. Swap colours 2 and 3 in the 2-3 component that contains z; z 
gets recoloured 3, while x , y and u retain their colours. Now, colour 2 is free at v, and 
the situation is similar to (II) above. So, swap colours 1 and 2 in component A; u gets 
colour 2, but x, y and z retain their colours. Now v can be coloured 2, the same as u. 
Case C2: x is disconnected from u in G - {y, z} 

Let us denote the component that contains u in G - {y,z} by 5. 

Since every path from x to u in G contains y or z, and neither y nor z is a 
cut-vertex of G (otherwise, the situation falls under Case Cl), there must be a path 
from z to u that does not contain x, y or u, and a path from y to u that does not 
contain z, x or v . Hence, z and u are connected in G — {u,x, y}. Also, y and u are 
connected in G - {u,z,x}. 

Since v is coloured 2, x, y and z are coloured 1 or 3. If z is coloured 1, then use 
Proposition 'Piu, z, x, y) to recolour z without affecting x, y or % otherwise, continue. 
If y is coloured 1, then use Proposition *P(u,u,2/,z,x) to recolour y without affecting 
z, x or u ; otherwise, continue. Uncolour u. 

The colouring now satisfies the following conditions: v is uncoloured; u is 
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coloured 1; x is coloured 1 or 3; y is coloured 2 or 3; 2 is coloured 2 or 3. 

Since v is uncoloured, one of the following four exclusive possibilities must hold. 

(I) . Colour 1 is free at v. Here, v can be coloured 1 , the same as u. 

(II) . Colour 2 is free at v. Here, both y and 2 must be coloured 3. Swap colours 1 
and 2 in component B\ u gets colour 2, but 1 , y and 2 retain their colours. Now v can 
be coloured 2 , the same as u. 

(III) . Colour 3 is free at v. Here, both y and 2 must be coloured 2 . Swap colours 1 
and 3 in component B\ u gets colour 3, but x , y and 2 retain their colours. Now v can 
be coloured 3, the same as u. 

(IV) . No colour is free at v. Since neither y nor 2 is coloured 1 , x must be coloured 
1 . Without loss of generality, assume that y and 2 are coloured 2 and 3 respectively. 

First assume that x and 2 are not in the same 1-3 component. Swap colours 
1 and 3 in the 1-3 component that contains 2 ; x gets recoloured 3, while y, z and u 
retain their colours. Now, v can be coloured 1, the same as u. 

Next assume that x and y are not in the same 1-2 component. Swap colours 

1 and 2 in the 1-2 component that contains 2 ; x gets recoloured 2 , while y, 2 and u 
retain their colours. Now, v can be coloured 1, the same as u. 

Next assume that y and 2 are not in the same 2-3 component. Swap colours 

2 and 3 in the 2-3 component that contains y; y gets recoloured 3, while 2 , 2 and u 
retain their colours. Now, colour 2 is free at v, and the situation is similar to (II) above. 
So, swap colours 1 and 2 in component B\ u gets colour 2 , but 2 , y and 2 retain their 
colours. Now v can be coloured 2, the same as u. 

Now the only case remaining is where 2 and y are in the same 1-2 component 
C 12 , and 2 and 2 are in the same 1-3 component C 13 , and y and 2 are in the same 2-3 
component C 23 . Here C 12 must be a chain with x and y as its endpoints; if C 12 has 
vertex w with only one colour in the neighbourhood, then w can be given colour 3, thus 
ensuring that 2 and y are in different 1-2 components. Similarly, C 13 must be a chain 
with 2 and 2 as its endpoints, and C 23 must be a chain with y and 2 as its endpoints. 
Also, C 23 does not degenerate into a single edge, because otherwise, u and v would be 
disconnected in G. Note that C 23 — {j/t z } I s a subgraph of component B. 

Swap colours 1 and 3 in C 13 ; z and 2 are now coloured 3 and 1 respectively. 
Since C 23 — {y, 2 } is a subgraph of component B, all its vertices retain their colour. 
Now, the only neighbour coloured 2 of 2 is in C 23 - {y,z}, and hence in 5; the other 
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two valencies of 2 are taken up by v and a vertex coloured 3. Thus, y and 2 are in 
different 1-2 components. 

Swap colours 1 and 2 in C u - x, the 1-2 component that contains y; y is now 
coloured 1. Now, x has two neighbours of colour 1, one each from C u and Ci 3 . 

Let a, b and c be the three neighbours of u, in G. Since, u is coloured 1, all its 

neighbours must be coloured either 2 or 3. We consider the four exclusive possibilities 
one by one: 

(i) . <z, b and c are all coloured 3. Colour both u and v with 2. 

(ii) . a, b and c are all coloured 2. Recolour x with 2; now colour 3 is free at v. Colour 
both u and v with 3. 

(iii) . Two of a, b and c are coloured 3. Without loss of generality, assume that a, b 
and c are coloured 2, 3 and 3 respectively. Let v be coloured 2. Then we can assume 
that u and v are in the same the 1-2 component, because otherwise, swapping colours 
1 and 2 in one of the components will give the same colour to u and v. Since u’s only 
neighbour of colour 2 is a, the 1-2 component that contains u and v must also contain 
a. Thus, by Proposition V(v,u,a,b,c), a can be recoloured without also recolouring b , 
cot v. Irrespective of whether the new colour a is 1 or 3, u can be coloured 2, the same 
as v. 

(iv) . Two of a, b and c are coloured 2. Without loss of generality, assume that a, b and 
c are coloured 2, 2 and 3 respectively. Recolour x with 2; now colour 3 is free at v. Let 
v be coloured 3. Then we can assume that u and v are in the same the 1-3 component, 
because otherwise, swapping colours 1 and 3 in one of the components will give the 
same colour to u and v. Since u's only neighbour of colour 3 is c, the 1-3 component 
that contains u and v must also contain c. Thus, by Proposition V(v ,u,c,a,b), c can 
be recoloured without also recolouring a, b or v. Irrespective of whether the new colour 
c is 1 or 2, u can be coloured 3, the same as v. 

Case C3: x and u axe connected in G — {t/, 2 }, and 
y and u axe connected in G — {s, z}, and 
2 and u are connected in G — {s, y} 

If a(x) = 1, then we use Proposition V(u,v,x,y, z) to recolour x without 
affecting y , z or u. If a{y) — 1, then we use Proposition T(u, v, y, 2 , x) to recolour 
y without affecting x, 2 or u. If a ( 2 ) = 1, then we use Proposition V(u,v,z,x,y) to 
recolour y without affecting x, y or u. None of x, y and 2 are now coloured 1. Thus, v 
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can be coloured 1, the same as u. 

Thus, given that Proposition V{u,v,x,y,z) is valid, we have Theorem 5.3, 
The proof Proposition V(u,v,x,y,z) is given in Section 5.2.2 

To prove Theorem 5.4, let G = (V) E) be a 3-degree Brooks’ graph on n 
vertices, and let u and v be two of its vertices apart by more than 362 edges. Let x be a 
neighbour of v. Then u and x are apart by more than 361 edges. By Theorem 5.3, there 
is a 3-colouring of G that colours x and u the same; it and v are coloured differently in 
this colouring. Hence Theorem 5.4. 


5.2.2 Proof of the Proposition 


We now prove Proposition V{u,v,x,y,z) 

Suppose x and u are connected in G - {v,y,z}. Consider the component that 
contains x in G — {v,y,z,u}. Of all breadth first search (BPS) trees rooted at x, of 
this component, pick one that has the largest number of leaves. Note that there 
could be arbitrarily large graphs H such that for a given vertex r in ff, all BPS trees 
of H rooted at r have a constant number of leaves (see Figure 5.1). Let T = {Vt,Et) 



OOO 


o-oo-o 


Figure 5.1: All BFS trees at r have a constant number of leaves 

be the BFS- tree thus obtained. This means that if {y, z, u} is a cutset in G — u, we 
have left out from T precisely those vertices in G — v reachable from x only through y, 
z or it. Note that T contains at least one neighbour of u, and hence contains at least 
361 vertices. 

For a vertex s G Vt, 1st p(s) denote the parent of s in T. Also, for 6 Vt, let 
P(s, t) be the path in T from s to t and d(s,t) be the shortest distance in G- {v, y, z, u } 
from s to t. Since T is a BFS-tree, for all s € V T> the length of P(x t s) is d(x, s ). 

If for a vertex s 6 V T , P(®, s) is a fork in G, then this fork can be used to 
recolour x without affecting y, z or u. Hence, we can make the following assumption: 
Assumption 1: For each s € Vt, P{x, s) is not a fork in G. 
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In other words, if si, S 2 and S3 are the neighbours of s, s* = p(s), then 
cr(s 2 ) or(s 3 ). 

Let l be a node of T that is at least three edges away from x, y, z and u in 
G — v. Since the combined valency of x. y, z and u in G — v is only nine, there are at 
most nine vertices one edge away, and at most 18 vertices two edges away from x, y, z 
or u. Thus, if T has more than 28 vertices, a vertex l as required can always be found. 
Let e = p(l ), / and g be the neighbours of l in G. Note that e, f and g are present in 
T, because all are reachable from x without going through y. z or u. Since P{x, l) is 
not a fork, one (and only one) of / and g, say /, is coloured the same as e. Without 
loss of generality, assume that <r(e) = <r(/) = 1, cr(/) = 2 and a(g) = 3. 

We show that for certain appropriate choices of l, one of which will always be 
possible, we can find a fork at x that consists only of a path P(x, w) in T, for some 
w € Vti and some vertices at most two edges away from l. Since l is at least three 
edges away from y, z and u in G - v, such a fork clearly will not contain y, z or u. and 
hence can be used to recolour x independently of them. 

We make the following claims that will be subsequently used in the proof. 
Claim 1: For w\,w 2 6 Vy, if W\ € P(x,w 2 ) then the path from w\ to xv 2 in T is one of 
the shortest paths from w\ to w 2 in G - {v,y,z,u}. 

Proof ojf_ Claim, 1: T is a BFS-tree. 

Claim 2: For wi,w 2 € Vt , if Wj € P(x,w 2 ) and {uq,^} € E, then Wi = p(w 2 ). 

Claim 3: For wi,w 2 € Vp, w\ 7 ^ ^2) if none w i’ s children is in P(x,w 2 ), then wj is in 
not P(x,w 2). 

Continuing with the proof of the theorem, if possible, let l be a leaf. For 
this particular case we make the following claims: 

Claim 4: l 4- f) an< f f ^ P( x i9) 

Proof of Claim f: Since I is a leaf in T . it is not an ancestor of / or g. 

Claim 5: If {e,g} E then e £ P(x,g) 

Proof 0/ Claim 5: Suppose {e,y} E and e € P(x,g). Then we have the following 
sequence of statements (see Figure 5.2(a)): 


0) 

d(e,g) < 2 

: l is adjacent to e and g in G 

(ii) 

<f(e,y) > 2 

: {e,g}tE 

(in) 

d(e,g) = 2 

: from (i) and (ii) 

(iv) 

there exists a vertex h ^ / in T 
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(a) 



/ 



(c) 


Figure 5.2: (a) Claim 5: P(z,e) is a fork (b) Claim 6: P(x,f ) is a fork (c) Claim 7: 
(P(x,g),l) is a fork 


so that p(g) = h and p(h) = e 

(v) p(h) =p(l) = e 

(vi) a{h) = a(l) = 2 

(vii) P{x , e) is a fork; a contradiction 


e 6 P{x,g) ) from (iii) and Claim 1; 
l is a leaf 
from (iv) 

<r(e) = 1 and a(g) = 3 

from (v), (vi) and Assumption 1 


Claim 6: If {/,£?} £ E then / £ P(x,g) 

Prop/ of Claim, 6: Suppose {f,g} 0 E and / £ P(z,ff). Then we have the following 
sequence of statements (see Figure 5.2(b)): 


(i) 

d(f,g)< 2 

: l is adjacent to / and g in G 

(ii) 

d(f,9) > 2 


(iii) 

<*(/,<?) = 2 

: from (i) and (ii) 

(iv) 

there exists a vertex hj^linT 



so that p(g) = h and p(h) = / 

: f £ P(x,g), from (iii) and Claim 1 

: l is a leaf 

(v) 

h £ P(a:, /) and l g P(x, f) 

: from (iv), Claim 4 
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(vi) h and l are adjacent to / in G : by definition 

(vii) a(h) = tr(l) = 2 : a(f) = 1 and <x(g) = 3 

(ix) P( x i f) is a fork; a contradiction : from (v), (vi), (vii) and Assumption 1 

In the sequel, for r distinct vertices w\ , — ,w r none of which is in P(x,w), 
by (P(x,w), w x,...,w r ) we mean the vertex disjoint path obtained by concatenating 
P{x,w) with edges {^ 1 ,^ 2 }, ...,{tu r _i,u; r }, assuming that all the required edges exist 
in G. 

Claim 7: If none of /, e and / is in P{x,g) then (P(x,g),l) is a fork. 

Prop/ of_ Claim 7: See Figure 5.2(c). 

(i) e and / are adjacent to / in E : by definition 

(ii) e and / are not in (P(x,g),l) : by the hypothesis of the claim 

(iii) a(e) = cr(/) = 1 : from the choice of a 

(iv) ( P(x,g),l ) is a fork : from (i), (ii) and (iii) 

Continuing with the proof of the theorem, we consider the following four 
exclusive possibilities, one by one. For each case we find a fork that can be used to 
recolour x without affecting y, z or u. 

Case Al: {e,<?} E and {f,g} £ E 

By Claim 4, / 0 P(x,g). Since {e,£f} £ E, by Claim 5, e £ P{x,g). Since { f,g } E , 

by Claim 6, / ^ P(x,g). Hence, by Claim 7, (P(x, </),/) is a fork. 

Case A2± {e,g} ^ E and {/, g) € E 

By Claim 4, / £ P(x,g). Since '{e,p} £ E, by Claim 5, e i P(x,p). Hence, if 
/ ^ P(x,g), then, by Claim 7, (P(x,g),l) is a fork. 

Assume that / € P{x,g). Since {f,g} € E, by Claim 2, / is the parent of g 

in T. Let h be a neighbour of g, different from both / and /. Since g is coloured 3, h 

must be coloured either 1 or 2. (see Figure 5.3(a).) We have the following sequence of 
statements: 

(i) l<£P(x,g) : Claim 4 

(ii) h£P(x,g) : {h,g} is in G, p(g) = f, Claim 2 

(iii) h is adjacent to g in G '• hy definition 

(iv) l is adjacent to g in G : by definition 

(v) if <r(/i) = 2 then 
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Figure 5.3: (a) Case A2, / = p(g): (P(x,l),g) is a fork (b) Case A3, e = p(g): 
(P(xJ),l,g) or (P(x,h),gJ) is a fork (c) Case A4: (P(z,e),s,/) or { P(x,l),f ) is 
a fork 


P(x, g) is a fork 

(vi) a(h) = 1 

(vii) cr(e) = cr(/) = 1 

(viii) / is adjacent to g in G 

(ix) f<£P(x,l) 

(x) d(h,l) = 2 

(xi) if h G P(x, l ) then h = p(p(l)) 

(xii) if h G P(x,l) then h = p(e) 
(xiii) h<£P(x,l) 

(xiv) g<£P(x,l) 

(xv) / and h are not in (P(x,l),g) 

(xvi) (P(x,l),g) is a fork 


: a(l) = 2, (i)-(iv) 

: (v) and Assumption 1 
: the choice of a 
: the hypothesis of Case A2 
:{/,/}€ J5, p{l) = e, and Claim 2 
:{h,l} & E and (iii) and (iv) 

: (x) and Claim 1 
: (xi),e = p(l) 

: (xii); from (vi) and (vii), a(h) = <r(e) 
: {g, 1} G E, p(l) = e, and Claim 2 
: from (ix) and (xiii) 

: (iii), (vi)-(viii), (xiv) and (xv) 


Case A& {e,g} G E and {/,</} 0 E 

By Claim 4, l £ P{x,g). Since {f,g} £ E, by Claim 6, / £ P{x,g ). Hence, if 
c £ P(s, <?), then, by Claim 7, (P(x,g),l) is a fork. 

Assume that e G P(x,g). Since {e,g} G E, by Claim 2, e is the parent of g in 
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T. (See Figure 5.3(b)). Let h be a neighbour of g, different from both e and /. Since 
g is coloured 3, h must be coloured either 1 or 2. We have the following sequence of 
statements: 


0 ) 

(ii) h £ P(x,g) 

(iii) h is adjacent to g in G 

(iv) / is adjacent to g in G 

(v) if a(h) = 2 then 
P(x,g) is a fork 

(vi) <r{h) = 1 

(vii) cr(e) = cr(f) = 1 
(viii) / $ P(x,f) 


: Claim 4 

: { h,g } is in G, p{g) = e, Claim 2 
: by definition 
: by definition 

: a (l) = 2, (i)-(iv) 

: (v) and Assumption 1 
: by the choice of a 
: Claim 4 


Now there are two cases: (a) h £ P(x,f) and (b) h € P(x, f). 

Consider the case of h £ P(x,f ); then the following statements are true. 


(ix) g £ P( x, f) : h is the only (if there is any) child of g, 

■ h i P( x ,f), Claim 3 

(x) e £ P(x , /) : (viii), (ix) and Claim 3 

(xi) e is adjacent to g in G : the hypothesis of Case A3 

Therefore, from (iii),(vi)-(xi), we have that if h £ P(x,f ) then (P(x,f),l,g) is a fork. 
On the other hand, when h 6 P(x,f), the following statements are true. 


(xii) / i p(x,h) 

(xiii) e £ P(x,h) 

(xiv) e and / are adjacent to fin G 

(xv) g and l are not in P(x, h ) 

Therefore, from (vii) and (xii)-(xv) we 
fork. 


: heP(x,f) 

: otherwise, e is an ancestor of h, which 
: in turn is an ancestor of /; since 
: d(e,f) = d(e,h) = 2, this cannot be 
: by definition 
: (xiii), p(g) = p(l) = e 

have that if h € F(a:,/) then (P(x,h), g,l) is a 


Thus, we can find a fork irrespective of whether h is in P(x,f) or not. 
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Case A£_ {e,g} e E and {f,g} 6 E 

Let k be a neighbour of /, different from both g and /. Since f is coloured 1, h must 
be coloured either 2 or 3. (See Figure 5.3(c)). We have the following sequence of 
statements: 


(i) e, l and / are adjacent to g 

(ii) g, h and l are adjacent to / 

(iii) either p{g) = / or p(g) = e 

(iv) either p(f) = h or p(f) = g 


: hypothesis 
: hypothesis 
: (i)i P(d) ^ l 
: (ii); p(f) jL l 


Now there are two cases: (a) h € P(x,e) and (b) h £ P(x,e). 

Consider the case of h £ P(x,e); then the following statements are true. 


(v) e?P(x,h) 

(vi) heP(xJ) 

(vii) d(h, l) = 2 

(viii) h = p(p(l)) = p(e) 

(ix) if p(f) = g then p(p(p(f))) = h 
else p(f) = h 

(x) h<=P(x,f ) 

(xi) f P(x,h) 

(xii) P(x, h) is a fork 


h 6 P(x,e) 
h 6 P{x,e ); e = p(l) 
f is adjacent to both h and l 
(vi), (vii), Claim 1 

(iii), (viii), (iv) 

(ix) 

(x) 

(v), (viii), (xi), <r(e) = a(f) = 1 
he following statements are true. 


Hence, by Assumption 1, h P(x,e)‘, and the 


(xiii) h £ P(x,e) 

(xiv) l?P(x,e) 

(xv) if p(g) = / then 

p(p(g )) = h, and so g P(x,e) 
else p(g) = e 

(xvi) g<£P(x,e) 

(xvii) f#P(x,e) 

(xviii) if a(h) — 2 then 

( P(x,e),g,f ) is a fork 

(xix) f,h,g $ P(x,l) 

(xx) if c(h) = 3 then 


e = p(l) 

(iv), (xin) ^ 

(iii) 

(xv) 

(ii), (xiii), (xiv) and (xvi) 

(ii); (xiii), (xiv), (xvi), (xvii); a(l) = 2 
(xiii), (xvi), (xvii); e = p(l ) 
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(P{x, l), f) is a fork : (ii), (xix); cr{g ) = 3 

Since h is coloured either 2 or 3, we get a fork any way. 

Thus, if the vertex /, chosen to be at least three edges away from x, y, z and 
u in G - v, happens to be a leaf, then we can obtain a fork at x that does not involve 
y, z or u. So, assume that all leaves of T are at most two edges away from x, y, z or 
u in G — v. Since the combined valency of x, y, z and u in G - v is only nine, there 
are at most nine vertices one edge away, and at most 18 vertices two edges away from 

x, y, z or u. Thus, there are at most 27 vertices less than three edges away from x, 

y, z or u. As can be readily verified, a maximum of 18 among these can be leaves in 
T. In any tree the number of vertices with two or more children cannot exceed the 
number of leaves. As T has at most 18 leaves, it has at most 18 vertices with exactly 
two children; note that T is a binary tree. In other words, all but at most 36 vertices 
in T have exactly one child. 

Let / be a vertex such that all vertices at most two edges away from / have 
exactly one child in T. If T has more than 36(1 + 3 + 2x3) = 360 vertices, then we 
indeed have a choice for l. As mentioned earlier, T has at least 361 vertices. 

Since / has exactly one child in T, either p(f) = l or p(g ) = l, but not both. 
Case BU f is the child of l 

Here, l ^ p{d)- (See Figure 5 . 4 (a)). We have the following sequence of statements: 


(i) / 0 P(x,< 7 ) 

{1,9} Claim 2 

(ii) f<£P(x,g) 

: (0, P(f) = l 

(iii) e $ P(x,g) 

: (i), l is the only child of e, Claim 3 

(iv) ( P(x,g ), l) is a fork 

: from (i)-(iii), Claim 7 


Case B2: g is the child of l 

Here l ^ p(f). Let h and a be the neighbours of g , other than Z, in G. Since P(x,g) 
is not a fork, &(h) ^ cr(a). As g is coloured 3, without loss of generality, let or(h) — 1 
and cr(a ) = 2 . Since g is a neighbour of / in G, g has exactly one child in T. 

First assume that a is the child of g. Then h is not a child of g. (See Figure 5.4(b)). 
We have the following statements: 

(i) g £ P(x,h) 

(ii) a£P(x,h) 


: g jz p(h ), {g,h} 6 E , Claim 2 
: (i), p(a) = g 
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(a) (b) (c) 

Figure 5.4: (a) Case Bl: (P(x,g)J) is a fork (b) Case B2, a is the child of g: 
(P(x,h),g) is a fork (c) Case B2, h is the child of g, f & P(r,fl): (P(x>a),g,l) is a 


fork 



(iii) 

ItPfah) 

: (i), g is the only child of Z, Claim 3 

(iv) 

l and a are adjacent to g 

: by definition 

(v) 

<r(l) — <r(o) = 2 


(vi) 

( P{x,h),g ) is a fork 

: from (i)-(v) 


Now, consider the case where h is the child of g, and so, a is not a child of g. 

If / ^ P{x, a) then (see Figure 5.4(c)) we have the following sequence of statements: 

(i) g£P(x,a) -.g? p(a) . {s> a} 6 E , Claim 2 

(ii) l P(x, a) : (i)> 9 is the only child of l , Claim 3 

(iii) e £ P(x, a) : (ii), l is the only child of e, Claim 3 

(iv) e and / are adjacent to l : by definition 

(v) fr(e) = cr(/) = 1 : by the choice of <r 

(vi) (P(x ) a),gJ) is a fork :from(i)-(v) 
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So, the only case remaining is where / € P(x, a) and h is the child of g. 



(a) (b) 

Figure 5.5: (a) Case B2, h is the child of g, f =p(p(a)) (b) Case B2, h is the child 

of g, f = p(a): P(x,f ) is a fork 

Since l is adajcent to / in G[Vx], \ d(x, l ) - d(x, f) |< 1. 

• Assume that d(®, l) = d(x> f) - 1. Let f be the parent of /; so, d(®, l) = d(x, f). 
Obtain a tree T\ from T by redefining p(f) as l. Note that T\ is also a BFS-tree 
of G[V t ] . Since f is two edges away from l in G[v T ], f is its only child in T; that 
is, /' is a leaf in T\. Thus, T\ has one leaf more than T. But T, by definition, 
ha.g the largest number of leaves of all BFS-trees rooted at x of G[Vx}- We have 
obtained a contradiction. 

• Assume that d(x,l) = d(x,f) + 1. Since e be the parent of l, d(x,e) = d{x,f). 
Obtain a tree T 2 from T by redefining p(l) as /. T 2 is also a BFS-tree of G\Vr], 
and has one leaf more than T. We derive a contradiction, once again. 

Hence it must be that d(x,l) = d(x,f). 

Since l = p{g) and / <= P(x,a ), d(x,a) > d(xj) = d(x,l) = d(x,g) - 1. 

That is, d(x,a) > d(xj) + 1 = d(x,g). But g is adajcent to a in G[V T j; hence 
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d(x,a) < d(x,g ) + 1 — d(x,f) + 2. So, either / is the parent of a, or / is the grand 
parent of a , in T. 

• Assume that / = p(p(a)); let a' be the parent of a. Clearly, d(x,g) = d(x,a'). 
Obtain a tree T$ from T by redefining p(a) as g. Since a' is two edges away from 
/ in <7[vt], a is its only child in T; that is, a' is a leaf in T 3 . Thus, T z has one leaf 
more than T — a contradiction. (See Figure 5.5(a)). 

Hence it must be that / = p(a). 

We have the following sequence of statements: 


(i) a is adjacent to / 

(ii) l is adjacent to / 

(iii) l and a are not in P(x, f) 

(iv) a(l) = <r(a) = 2 

(v) P(x,/) is a fork 


/ = P(«) 
by definition 

/ = p(a); d(x,l ) = d(x,/) 
from (i)-(iv) 


(See Figure 5.5(b)). But, by Assumption 1, this is a contradiction. 
That proves Proposition V{u, v, x, y , 2 ). 



76 


Chapter 6 

Colouring of Interval Graphs 


A graph G = (V,E) is called an interval graph, if for some set 5 of intervals 
of a linearly ordered set, there is a bijection / : V — + Sr so that two vertices u and 
v are adjacent in G iff f(u) and f(v) overlap. Every interval graph has an interval 
representation in which the endpoints are all distinct [36]. Hence, G can be represented 
by a set of endpoints £ of size 2 | 5 |, where, for each I € there are unique elements 
1(1) and r(I) in £ corresponding respectively to the left and right endpoints of I. We 
define an “inverse” function 2 which gives the corresponding interval for 
each member of £. That is, if e is either 1(1) or r(I) and I € 9 then, 2(e) = I. We 
shall often find it convenient to represent an interval graph as directed, so that for two 
overlapping intervals I\ and h, (h,h) € E iff 1(E) < 1(h). 

Interval graphs are perfect graphs [36]; that is, for an interval graph G = 
(F, E), and for every induced subgraph G' of G, the chromatic number of G' is equal 
to the clique number of G'. 

In this chapter, we consider some parallel algorithmic issues of minimally 
colouring interval graphs. We use the short form “IGC” to denote the problem of 
minimally colouring an interval graph with a known interval representation. 

If in an interval representation of an interval graph, all endpoints are not 
distinct, then an equivalent one of distinct endpoints can be constructed by a constant 
number of invocations to the All Nearest Smaller Values problem [94], provided the 
endpoints are sorted. We briefly outline this algorithm. Let an interval graph G = 
(V,E), be represented by two arrays L and R containing the left and right endpoints 
of the intervals of G given each in a non-decreasing order, either sorted or padded- 
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sorted. The elements of L and R can be merged into an array A in 0(1) time with n 1+e 
processors on a COMMON CRCW PRAM [96, 92]; if endpoints are in padded-sorted 
form, then the merged array will also be in padded-sorted form. Suppose there are more 
than one occurrences in A of an endpoint s. Say, the left most one occurs at location 
i and the right most one at location j. Let r (resp., t) be the nearest smaller (resp., 
larger) on the left (resp., right) of s in A. For, 0 < k < (j — i), change the endpoint 
of location (i -f k) to $ + k * 2 (j~i+i) ^ ^ a right endpoint, and to + k * 2 (f~i T + i) 
if it is a left endpoint. Thus, constructing the new set of endpoints, which are clearly 
distinct, needs only 0(1) time with n 2 processors, or alternatively O(loglogn) time 
with n/loglogn processors on a COMMON CRCW PRAM [11], and hence is in NC 1 . 
So, in the sequel, without loss of generality we assume all 2n endpoints to be distinct. 

6.1 IGC on Bounded Fan-in Circuits 

The complexity of IGC on bounded fan-in circuits is considered in this section. 
It is shown that IGC may not be in NC 1 even when the left and right endpoints are 
separately sorted. 

A linked list is a directed graph in which both the out-degree and the in-degree 
of a vertex can be at most one. In a monotonic linked list v can be an out-neighbour 
of u only if u < v. We will denote by “MLCC”, the problem of labelling each node in 
a collection of disjoint monotonic linked lists, by the first element of the corresponding 
list; that is, MLCC is the problem of finding the connected components in a collection 
of disjoint monotonic linked lists. 

6.1.1 A Reduction from MLCC to IGC 

Consider an instance of MLCC of size n with k components, given in an array 
A[1 .. .re]. From this instance we construct an equivalent instance of IGC as follows. 

Initialise an array R[1...3n] with zeros. Replicate the lists in A in the subar- 
ray B[n+ 1...2n]; that is, A[t] is copied into B[n + i], and if A[i] points to A\j], then 
a pointer is set from B[n -1- i] to B [n -f j] , for 1 < i , j < n. For n + 1 < i < 2n, if B[i] ’s 
out-degree is zero, add an edge from R[i] to B[i + n] and if J?[i]’s in-degree is zero, 
add an edge from B[i - n] to B[i]. Now, B[1 . . . 3n] is a collection of monotonic linked 
lists, where each of the k non-zero elements of B[l...n] (resp., of B[2n + 1...3n]) 
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is the starting point (resp., the ending point) of a list. Using prefix summation, in 
NC 1 , compact B[ 1 . . . 3n] into B[n - k + 1 . . . 2n + k] without disturbing the order of 
its elements. 

Construct a set $ of intervals, so that, for n - k + 1 < i,j < 2 n + k. if B[i\ 
points to B\j] then, [i,j — |] € $. Note that both the left and right endpoints of these 
intervals can now be represented in sorted form in two separate arrays, say L and R. 

Consider the interval graph G defined by 5. Identifying each component of 
B with a unique colour, we get a valid colouring of G\ that is, x{G) < k. Now, 
suppose that there is an optimal colouring of G that for two consecutive edges (u,v) 
and ( v , w) of B, gives different colours to their corresponding intervals I u = [u.v - |] 
and I v = [v,w - |] in G. Since all endpoints are distinct, and there is no endpoint 
between v — ^ and v. both I u and I v share k — 1 mutually adjacent neighbours, and 
we get x(G) = k + 1 > k, a contradiction. That is, an optimal colouring of G will use 
one colour per component of B . In other words, x{G) — k and any optimal colouring 
of G will identify the connected components of B, and hence of A. Thus we have the 
following lemma. 

Lemma 6.1 An instance of MLCC of size n with k components can be reduced to an 
instance of IGC of size 0 ( n) and chromatic number k, in NCf ■ 


6.1.2 A Reduction from IGC to MLCC 


A sequential algorithm for IGC is easy to visualise. Let Q be a queue of size 
X(GQ, which is initialised with all the available x(G) colours. Consider the endpoints of 
the intervals one by one in non-decreasing order. For each left endpoint encountered, 
remove a colour from the front of Q and colour the corresponding interval with it, 
whereas, for each right endpoint, release the colour of the corresponding interval onto 
the back of Q. When the last of the left endpoints has been considered, the graph 

would have been coloured. 


We attempt a direct parallelisation of this algorithm. Let L = {h , . •• ,U and 
R = {n, . . . , r „) respectively be the sets of the left and right endpoints of the intervals 
given in non-decreasing order. For 1 < * < », let U be the rank of r f in L. That is, 
l < < i t < ri < l t . +l . So, when r f releases the colour of the interval I(r,-) onto 

the back of Q, U left endpoints and i right endpoints would have been encountered and 
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the length of Q would be x(C) ~ U + i . Hence, the colour released by r,- will be taken 
up by the (x(G) — t{ + i)-th left endpoint from now on; that is, by the left endpoint 
^U+x(G)-U+i = l x (G)+i ■ I n other words, both I(r,) and X(/ x (< 3 ) + ;) are to get the same 
colour and no interval with a left endpoint between their respective left endpoints will 
get that colour. 

Define a set of pointers p : L -*• L as follows: let p(/(X(r;))) be l x (G)+i, f° r 
1 < i < rz. It is clear that p defines a collection of monotonic linked lists over L. Once 
we find the connected components in this collection we have got an optimal colouring 
of the graph. Thus we have the following lemma: 

Lemma 6.2 IGC can be NC 1 -reduced to MLCC. 

Theorem 6.1 IGC is equivalent to MLCC with respect to NC 1 -reductions. 

6.1.3 A Reduction from List-Colouring to MLCC 

A general linked list L can be visualised as being constituted of alternating 
stretches of forward and backward pointers in an array, and hence can be decomposed 
into two instances of MLCC, one consisting only of forward stretches and the other only 
of backward stretches. Let L 1 and L 11 be the MLCC instances formed by the forward 
and backward stretches of X, respectively. An MLCC algorithm can be used to 2-colour 
L' as follows: 

Find the co nn ected components in L and identify each component with 
the first vertex of the component. For each vertex u in L , let L v be the 
list obtained by removing, where one exists, the edge going out of v in L'. 

Find the connected components in L' v , and identify each component with 
the first vertex of the component. Count the vertices that are not in the 
same component in both V and L' v . Depending on whether the count is 
odd or even, let v have colour 1 or 2. 

Similarly, 2-colour L" also. Now, L is 4-coloured. A 4-colouring can be reduced to a 
3-colouring in NC 1 [32], 

From Theorem 6.1, thus, we get the following theorem: 

Theorem 6.2 The problem of Z-colouring a linked list can be reduced to IGC in NC 1 . 

Linial has shown that for a linked list of length n to be 3-coloured, each of 
its vertices has to scan a path of length ft(log*n) [72]- Besides, the problem of finding 



Chapter 6. Colouring of Interval Graphs 


80 


a path between two vertices apart by a distance of g(n ) in a graph is not in AC 0 
when g(n) —* oo with n; also, this problem is unlikely to be in NC 1 [98]. That is, an 
o(log nlog’ n) depth bounded fan-in circuit that 3-colours a linked list probably doesn’t 
exist. Thus, in view of Theorem 6.2, may be, optimally colouring interval graphs is not 
in NC 1 , after all. 

6.1.4 An Upper Bound 

Consider an instance L of MLCC with k components in it. Divide the input 
array into n/2k segments of size 2k each. Remove the inter-segment links from L. 
Now, each segment forms an instance of MLCC. Use O(logfc) iterations of recursive 
doubling on each segment to find the connected components in it. Since addresses in 
each segment are from a range of size 2k each iteration can be implemented in O(log k) 
time. Identify each component with the first vertex (called a root) in it. Compact the 
roots in each segment. Now each segment has a size of at most k. For each root v, let 
the nearest root, from its subsequent segments, that belongs to the same component as 
v be its out-neighbour; a new instance of MLCC of size at most n/2 is formed. Thus, 
in 0(lognlog 2 k) time all connected components can be found. Informing each vertex 
of the component it belongs to, will take only an additional 0(log n log k) time. 

Hence, when the number of components is restricted to 0(1), MLCC, is in 

NC 1 . 

From Theorem 6.2, thus, we get the following theorem, which complements 
the observation that 3-colouring a linked list, in general, is unlikely to be in NC 1 . 

Theorem 6.3 A linked list that has at most a constant number of stretches, can be 
2-coloured in NC 1 . 

6.2 IGC on PRAM Models 

First, we consider a problem closely related to IGC, namely, that of finding 
the chromatic number of an interval graph with a known interval representation. 

6.2.1 Finding the Chromatic Number 

Chen [16] proves that finding the chromatic number of an interval graph 
requires ft(logn/loglogn) time, through a reduction from Parity. The proof (with 
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slight modifications) runs as follows. Let A[l...n] be an instance of Parity. Con- 
struct (n + 1) intervals 7[0],...,I[n], by setting /[0] = [0,2n -f- 1] and for 1 < j < n, 
I[j ] = [2n + j + \,2n + j + 1] if A\j] = 0 and I[j\ = [j, 2n - j + 1] if A\j] = 1. Find 
the chromatic number x of the graph corresponding to these intervals. The parity of 
A[1 ... n] is even when x is odd and odd when it is even. 

However, if we use a different input representation, that is, if we assume that 
both the left and right endpoints are given to us in separate sorted sets, then it turns 
out that the chromatic number can be found in constant time. 

Let an interval graph G = ( [V,E ) be represented by L = {k,. ..,l n } and R = 
{t*i , . . . , r„ } the left and right endpoints of the intervals of G given in a non-decreasing 
order. The sequences L and R can be cross ranked in 0(1) time with n 1+£ processors 
[96, 92] or in O(loglogn) time with n/loglogn processors on a CREW PRAM [55]. 
For 1 < i < n, let s; be the rank of in R. That is, T\ < ... < r Si < l < r 5i+i . 
Then, i — Si (= U say), is the size of the clique formed precisely by those intervals that 
contain U- In other words, is the size of the clique that !(/,•) forms along with its 
in-neighbours. Observe that any maximal clique C, of an interval graph should contain 
a vertex v, such that the set of all in-neighbours of v is precisely the set C - {t>}. Thus, 
the chromatic number of G can be obtained by taking the maximum of t{ over all i. 
This again can be done in 0(1) time with n 1 "*"' processors [96, 92] or in O(loglogxi) 
time with n/loglogn processors [55, 92] on a COMMON CRCW PRAM. Thus, we 
have: 

Lemma 6.3 The clique number of an interval graph can be found in 0(1) time using 
n i+< processors, 0 < € < 1, on a COMMON CRCW PRAM, provided, the left and right 
endpoints of the intervals are given in sorted order, separately. 

Corollary 6.1 The clique number of an interval graph can be found in O(loglogn) 
time using n/loglogn processors on a COMMON CRCW PRAM, provided, the left 
and right endpoints of the intervals are given in sorted order, separately. 

Even when the endpoints are only padded-sorted, for a graph with a suffi- 
ciently small chromatic number the general lower bound of 0(log n/loglogn) can be 
surpassed. Recall that every vertex v together with its in-neighbours forms a (not nec- 
essarily maximal) clique. (See the arguments preceding Lemma 6.3.) Thus, for each 
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vertex v € F, using O(n) processors, we copy the endpoints of all of v's in-neighbours 
into an array N v [l . . .2n], without changing their relative positions. We further add 
the endpoints of v to the array N v , at their appropriate places. Now, N v contains 
2x v < 2x(G) non-zero entries in non-decreasing order. Compact this array, into 2x v lo- 
cations in 0 (log x{G)/ log log n ) time using n processors [83], and we would have found 
out x v . The chromatic number of the given graph can now be found by taking the 
maximum of x v over all v € V; this can be done in 0(1) time with n 1+£ processors, 

0 < e < 1. The overall time taken is only 0(logx(G)/ log log n ) with n 1+£ processors. 

This upper bound can be proved to be tight by giving a matching lower bound 
as well. We do this by showing a reduction from ordered compaction to the chromatic 
number problem. The ordered compaction problem [83] is defined as follows: given an 
array A of size n, in which only k elements are non-zero, place the k non-zero elements 
of A in an array B of size k, so that the elements appear in the same order in B as 
they do in A. Ragde [83] shows that, not only can ordered compaction problems of size 
n with k non-zero elements be solved in 0(log k/ log log n.) time with n processors on 
a COMMON CRCW PRAM, but a matching lower bound, with a polynomial number 
of processors, also exists. 

The reduction is as follows. Let A[1 . . . n] be an instance of ordered compaction 
that contains k non-zero elements. Create n arrays, A,- for 1 < i < n, each of size n 
by setting Ai\j] = A\j) when j < i and A x [j] = 0 otherwise. For, 1 < i < n, construct 
an interval graph Gi on n intervals Jj[l], — ,/*[«], where, for 1 < j < n, if Ai\j] / 0 
then Ii\j] = [j,2n -j + 1] and if A t [j] = 0 then /,'[;'] = [2n + j - |,2n + j]. Observe 
that intervals /i[r] and Ii[s] overlap iff A,[r] and A,[s] are both non-zero. Find the 
chromatic number x. of G t . Move A\i] to A [*.•]• The first k locations of A now contain 
the non- zero elements of the input in the same order. That is, a more appropriate lower 
bound for finding the chromatic number of an interval graph G is Cl(logx(G)/ loglogn) 
rather than just ft(log n/ loglogn). 

6.2.2 A Lower Bound for IGC 

From the discussion on the chromatic number problem, presented in the previ- 
ous section, a lower bound of fi(log X (G)/ loglogn) on time for IGC, with a polynomial 
number of processors, is obvious. We prove a tighter bound: even when the left and 
right endpoints of the intervals are sorted separately, an Q(logn/ loglogn) time lower 
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bound holds for IGC. 

Lemma 6.4 An instance of MLCC of size n with k — 0(1) components can he reduced 
to an instance of IGC of size 0(n ) and chromatic number k, in 0(1) time on CRCW 
PRAM using a linear number of processors. 

Proof: The proof is similar to that of Lemma 6.1, except in that instead of prefix 
summation, we use ordered compaction to compact B[l . . . 3n] into B[n-k+ 1 . . . 2n+k]. 
Since k = 0(1), the compaction takes only 0(1) time [83]. □ 

Lemma 6.5 An instance of MLCC of size n, even when the number of components is 
two, requires ft (log n/ log log ra) time on a CRCW PRAM, with a polynomial number 
of processors. 

Proof: We show that Parity can be reduced to an instance of MLCC with two com- 
ponents. Consider an instance of Parity of size n, given in an array X[l...n]. Let 
Y[1 . . . 2n + 2] be an array of vertices. Establish pointers on these vertices as follows: 

if X[i] = 0 p(Y[2i - 1]) = Y[2i + 1] 
p(Y[2i}) = Y[2i + 2] 
otherwise p(Y[2i — 1]) = Y[2i + 2] 
p(Y[2i}) = Y[2i + 1]. 

Now, there are exactly two monotonic linked lists in Y. Find the connected components. 
If Y[2n + 2] and 7[2] are in the same component the parity is even; otherwise the parity 

is odd. D 

Thus, we have the following theorem: 

Theorem 6.4 IGC has a lower bound of ft(log nj log log n) on time, on a CRCW 
PRAM, with a polynomial number of processors, even when the left and right endpoints 
of the intervals are sorted separately. 

6.2.3 Upper Bounds for IGC 

Every step of the NC 1 -reduction from IGC to MLCC, of Lemma 6.2, can be 
executed in 0(1) time on a COMMON CRCW PRAM, using a polynomial number 
of processors. Hence, any algorithm for MLCC indicates an algorithm for IGC with 
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the same asymptotic resource bounds, on a COMMON CRCW PRAM. Since MLCC 
is a special case of the list ranking problem, it follows that both IGC and MLCC have 
<9 (log n) time optimal algorithms [6]. An 0(log n) time, n processors, EREW PRAM 
algorithm for IGC has been found before [89]. 

The list ranking problem is FL-complete for NC 1 -reductions, and hence is 
unlikely to be in NC 1 [90]. Besides, it is known that every problem in NC 1 can be 
solved in 0(logn/log logn) time with a polynomial number of processors on a CRCW 
PRAM [53]. Hence the following questions arise (i) whether MLCC is FL-complete for 
NC x -reductions, and (u) whether MLCC can be solved in 0 (logn/ log logn) time on a 
CRCW PRAM with a polynomial number of processors. All these are as yet open. 

However, in this section, we investigate solutions for MLCC, by imposing 
restrictions on the number of components (equivalently, on the chromatic number, 
if we consider IGC). We show that, when the number of components is limited to 
<9 ((log ra) 1-£ ) for some e € (0, 1) an instance of MLCC can be solved in 0(log n/ loglog n) 

time with a polynomial number of processors. 

The 0(logn/ loglog n) time algorithm proceeds by reducing the problem size 
by a factor of 0((logn) £ ) in every step. The following lemmas are made use of by the 
algorithm. 

Lemma 6.6 [22, 91] A prefix sums computation over an array of n bits can be done in 
0(1) time using n2 n processors on a COMMON CRCW PRAM, with a preprocessing 
time of 0 (log n ] log log n) . 

Proof: A single step of table-lookup is used. The input is simultaneously compared 
with every possible input, and the table entry corresponding to the match gives the 
result. Since there are 2" possible inputs, n2" processors are required. The table can 
be constructed in 0 (logn/ log logn) time using 2" simultaneous copies of the standard 
prefix sums algorithm [22]. 

T.. 8.7 An instance of MLCC of sice n can be soloed in 0(1) time using 0(n :2») 

processors on a COMMON CRCW PRAM with a pre-processing time of 0( logn). 

Proof: Suppose an array s*[l...»] contains the input. Of all the 2* possible two- 
colourings of A, we choose a o : {A[l],.. .. A[n]} - {1,2} so that, for every vertex o 
with in-degree zero <r(t>) = 1. .There is a unique such two-colouring. Basically, a vertex 
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gets coloured 1 or 2 depending on whether its rank is odd or even counting from the 
beginning of its corresponding list. Hence, a can be found in 0(1) time using ra2 n 
processors on a COMMON CRCW PRAM. 

Make n copies Ai [1 . . . n], . . . , A„[l . . . n] of the input A. For the i-th copy, if 
A,[i] has in-degree zero, do the following. Make A,[0] the in-neighbour of A,-[i] and find 
a two-colouring <7 t - : {A,[0], . . . , A,-[n]} -+ {1,2} so that, every vertex with in-degree 
zero gets coloured 1; thus, only the vertices of the linked list that start at vertex i will 
get their colour changed. For 1 < j < n, if cr,-(A,-[. 7 ']) ^ <r(A[j]) then let <?(A[;']) = A[i). 
Clearly g(A[i]) = A[i], and we call A[i] a root. Now, q defines the connected components 
of A; that is, vertices u and w are in same component iff q(u) = q(w). 

The number of processors used is n(n + l)2 n+1 = 0(n 2 2 n ). P 

Now, the algorithm is easy to visualise and is embodied in the following the- 
orem: 


Theorem 6.5 An instance of MLCC of size n in which the number of components is 
limited to C(n ) can be soloed in time with 0(nt 2 f ) processors, for any 

t > C(n), on a COMMON CRCW PRAM. 

Proof: Let A[1 . . . n] contain the input instance. For some s < t, divide A into segments 
of size t each. Assign 0(t 2 2*) processors to each segment and solve the i-sized instance 
of MLCC defined by each segment in parallel in 0(1) time using Lemma 6.7. There can 
be at most C{n ) roots in each segment. Compact the roots of each segment in 0(1) time 
using Lemma 6.6. For each root t>, let the nearest root, among its subsequent segments, 
that belongs to the same component as * be its out-neighbour. The compacted roots 
along with the new edges define an instance of MLCC of size at most nC{n)/t. That 
is in 0(1) time the problem size is reduced by a factor of t/C{n). So, the connected 
components can be found in 0(logn/log(t/C(n))) = CKiogt-iogqn)) t * me ‘ tota * 
number of processors used is O^t^n/t) = 0(nt2 t ). 

Corollary 6.2 An instance of MLCC of size n in which the number of components is 
limited to C(n) can be solved in O(j^^j) time with 0(n2‘) processors, for any 
s > C(n), on a COMMON CRCW PRAM. 


Proof: Set t = s/2 in the Theorem. 
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Corollary 6.3 An instance of MLCC of size n, in which the number of components 
is limited to 0((logn) l-£ ) for some e € (0,1) can be solved in 0 (log n/ log log n) time 
with a polynomial number of processors on a COMMON CROW PRAM. 


Proof: Let C(n ) = 0((log n) 1 £ ) be the number of components in the input instance; 
use Theorem 2 with t = log n. □ 
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Chapter 7 

Edge-Colouring of Graphs 

The problem of (A+l)-edge-colouring an arbitrary graph is not yet known to 
be in NC; Karloff and Shmoys [60] give an algorithm that is in NC only when A is at 
most polylogarithmic in n. This algorithm runs in 

0( A 5 log n( A + log n + min{log 3 n, A 2 log A(log* n + A 2 )})) 

time on an EREW PRAM with n + m processors. Liang, Shen and Hu [71, 70] give 
an 0(A 4-5 log 3 Alogn + A 4 log 4 n) time, n 3 log A -f nA 3 processors CRCW PRAM 
algorithm for (A+l)-edge-colouring a general graph. 

Fiirer and Raghavachari [27] give edge-colouring algorithms that are fast, but 
wasteful in colours. They give an 0 (log n log A) time, n + m processors, CREW PRAM 
algorithm for cA-edge-colouring a general graph, c > 1, and an 0(log* n) time, n + m 
processors, CREW PRAM algorithm for A 2 -edge-colouring a general graph. 

In this chapter, we improve these results: the following algorithms for edge- 
colouring a general graph are presented: 

• an algorithm that finds a (A + d)-edge-colouring, 1 < d < A, in 0((logd 4- 
(A/d) 4 )log 2 n ) time, using n + m processors 

• an algorithm that finds a A 1+£ -edge-colouring, 0 < e < 1, in 0(log Alog(log* n)) 
time, using nA 1+e processors 

Also, we show that from Theorem 5.2 it would follow that a cubic graph can be 4-edge- 
coloured in O(logn) time with 0(n) operations on an EREW PRAM. 
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7.1 (A + c?)-Edge- Colouring, d > 1 

An 0((logd-|-(A/d) 4 )log 2 n) time, n+m processors EREW PRAM algorithm 
for (A + d) -edge- colouring an arbitrary graph, 1 < d < A, is presented in this section. 
(The basic idea of this algorithm was discussed, in the context of (A-f-l)-edge-colouring 
of a bounded degree graph, in the author’s M. Tech, thesis [85]; see also [86].) 

The algorithm is based on the proof of Vizing’s theorem by Misra and Gries 
[76]. Besides being faster than the 

0 ( A 5 log n( A + log n + min{log 3 n, A 2 log A(log* n + A 2 )} )) 

time algorithm of [60], our algorithm has a simpler analysis. 

7.1.1 A Proof of Vizing’s Theorem 

Here we briefly sketch a proof of Vizing’s theorem. It is essentially similar 
to the one by Misra and Gries [76], but some of the details are altered to suit for 
our purposes better. (For another proof see [26, 13].) The proof uses induction on 
the number of edges. Basis: every graph with A edges is (A + l)-edge-colourable. 
Hypothesis: every graph with m - 1 edges is (A + l)-edge-colourable. 

Let G = (1 V,E ) be a graph with m edges. For e = {v 0 ,vi} € E, obtain a 
(A + l)-edge-colouring of G - e; this is possible, by the hypothesis. A colour c is said 
to be free at v G V if none of the edges incident at v is coloured c. Clearly, at least one 
colour is free at each vertex of G. Assume that no and v\ do not have a common free 

colour, otherwise, e can be given that colour. 

Let a (c 0 ,Cfc)-fan F(v 0 ) = %,c*), for colours c 0 ,...,c fc and 

distinct vi,...,vk € N(v 0 ), be a data structure that satisfies the following conditions: 

• the colour c; is missing at n,- for 0 < i <k. 

• the edge (y 0 , n;) is coloured c,_i for 1 <i <k. 

• either Ck is missing at no or Cfc = Cj for some j such that 0 < j <k. 

Vertices v 0 = h(F ) and v k = t(F) are respectively called the head and tail of the fan 

F = F(v o). 

Since, for any colour, there is at most one edge of that colour incident with 
any vertex, a colour other than c fc can occur only once in the fan, whereas c k can occur 



Chapter 7. Edge-Colouring of Graphs 


89 


at most twice. If cj, is missing at vq (that is c;. occurs only once in the fan or c* = co; 
in this case the fan is called local ), then for 1 < i <j Ar, the edge {uo,u t } can be given 
the colour c,-, thus completing the proof. (This is called the fan operation). So, assume 
that the fan is not local; c k occurs twice in the fan, and Ck f Cq\ C’ K = cj for some j 
such that 0 < j < k. Let Vj = m(F) be called the middle of F. 

In an edge coloured graph any bichromatic component must be a path or a 
cycle. (A c 0 -Ck component is a maximal connected subgraph induced by edges coloured 
co or c fc .) 

The head and tail of F are the only vertices in F at which colours Co or Ck are 
free, and hence are the only vertices at which a co-c/. path can end. Since, in addition, 
F is a non-local fan, a co-c/. path indeed ends at each of them. These three Co-Cfc paths 
need not all be distinct. But, at least two of them are distinct, because, at least two 
paths are needed to account for three endpoints. 

Let H(F), M(F) and T(F) denote the co-c* paths ending at h(F), m(F) and 
t(F) respectively. Define e(F), the elect of F, as follows: if E(F) = M{F), then let 
e(F) = t(F), otherwise, let e(F) = m(F). Also, let E(F) = M(F) or E{F) = T(F) 
depending on whether e(F) = m(F) or e(F) = t(F). 

Interchange colours Co and Ck in either H(F) or E(F ), but not both. (This is 
the chain operation). 

If e(F) = t(F) and the chain operation was performed on E(F), then 
(vo, Co, vi , Ci . . . , Vk, Co) is now a valid fan, and it is local. 

If e{F) = t{F) and the chain operation was performed on E{F) then 
(vo, Ck, »i , ci . . . , Vk, Ck) is now a valid fan, and it is local. 

If e(F) = m(F) and the chain operation was performed on E(F) then 
(vo, co, i>i , ci . . . , Vj, co) is now a valid fan, and it is local. 

If e (F) = m(F) and the chain operation was performed on -ff(F) then 
(v 0 , ct, vi , Ci . . . , vj, Ck) is now a valid fan, and it is local. 

A fan operation will now complete the proof. 

Observation 7.1 Of the two distinct colours free at v\, only one has been used in the 
proof. Hence the proof would be valid even if v\ had only one free colour. 
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7.1.2 Creating Fans in Parallel 

Consider a graph G = (V, E) that is partially (A-f-l)-edge-coloured so that 
any two uncoloured edges of G are at least three edges apart. In this graph, if fan 
operations are performed at all uncoloured edges simultaneously, no two of them will 
interfere with each other. 

We now show how to construct fans at the uncoloured edges of G. We assume 
that G is given in adjacency list representation; thus, each vertex has a list N(v o) of 
all its neighbours. 

Let e = {i>o,Vi}, be an uncoloured edge; arbitrarily choose one endpoint, say 
vo, of e as the head of the fan F(v o) to be constructed. 

Step 1: Form a directed graph H'(v o) on {«o} U N(v o) as follows: Let (vo, Vi) be an 
edge in H'(v 0 ). For all u € N(v 0 ), let a(u) be a colour that is missing at u; if there is 
more than one colour missing at u, then we choose an arbitrary one of them as a(u). 
For w € N(v 0 ), if the edge {u 0 ,t«} in G is coloured a(u) then let ( u,w ) be an edge in 
H'{v 0 ). Clearly, the out-degree of every vertex in H'(v 0 ) is at most one. 

R.emark: Assume that each chosen head v of a fan has its adjacency list sorted on 
colours, and hence has its incident edges in an array Z(v) indexed by colours. With 
processors assigned one per neighbour of uo, each u £ N(v 0 ), can find out its out- 
neighbour in H\v 0 ) in 0(1) time by performing a concurrent read from Z(v 0 )[<t(u)]. 
On an EREW PRAM, this concurrent read can be simulated in O(logA) time [55]. The 
edge-lists, along with twin pointers, of H'(v 0 ), and hence of the underlying undirected 
graph JT(uo), can be formed in an additional 0(log A) time on an EREW PRAM. 

Step 2: Each component of H'(v 0 ) is a pseudo tree; there is at most one directed 
cycle in each component. Let C denote the component that contains v 0 in H'{v 0 ), and 
let C" denote the cycle in C. To find a fan at v Q , we need only to find the set of vertices 
reachable from u 0 in C. Find C and C using the Euler Traversal Technique of [95]; 
C - C" is a rooted forest. Use the parallel tree contraction technique [1] to find the set 
of vertices reachable from in C - C'-, this set together with C form the set of all 

vertices reachable from vq ’i nC. 

Note thet B'(vo) and the nnderljdng nndirected graph B[v o) have been 
formed in adjacency list representation, and that twin pointers are available; that », 
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the adjacency list entries [x, y] and [y, x] corresponding to an edge {x, y} in H(v 0 ) point 
to each other. Suppose next([x,^j) is the entry following [a: , y] in the adjacency list of 
x. Then the traversal T defined by the pointers “tour- next” is constructed as follows: 

for each adjacency list entry of H(v o) pardo 
tour-next ( [x , y] ) = next ([y,x]) 

Each component of H(vq) that has a cycle, contributes two lists to the traversal 
T. One (say L ') of the two lists contributed by C corresponds to C' — this list does 
not traverse vq. The other list (say L") contributed by C corresponds to C — C' and 
hence traverses vq. Using list ranking, identify L" among the lists of T; this is easy, 
because L" contains v 0 . Since L' is the only list in T that has vertices (of H(v 0 )) in 
common with L", L' can be now identified in 0(1) time on a CRCW PRAM, and hence 
in 0(log A) time on an EREW PRAM with A processors [55]. 

Now that C -C' has been identified, label v 0 with 1 and every other leaf of 
C -C' with 0. For each node of 0 - C', find the maximum among all its descendant 
leaves; all nodes for which the result is 1, along with C', belong to the subgraph C of 
C that is induced by vertices reachable from u 0 through a directed path. By [1] this 
takes O(log A) time on an EREW PRAM. 

Step 3: Clearly, the subgraph C" induced by the set of vertices reachable from v 0 
in C is either a simple directed path, or a simple directed path augmented by a vertex 
disjoint directed cycle. Also, C" gives a fan with v 0 as its head. When C" is a simple 
directed path, the fan is local too; the last vertex of the path is the tail of the fan. If 
the fan is not local, then the vertex with in-degree two in C" (there is exactly one such 
vertex) is pointed to by the middle and tail of the fan. 

Thus, we have the following lemma: 

Lemma 7.1 For a partially (A + l )-edge-coloured graph, where between any two tin- 
coloured edges of G there is no path of length less than three, a fan can be constructed 
at each uncoloured edge in parallel in 0(1 log A) time with 0(n + m) processors on an 

EREW PRAM. 

7.1.3 The Algorithm 

Consider a graph G = (V,jB) that is partially edge-coloured using the palette 
+ 1}, so that between any two uncoloured edges of G there is no path of 
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length less than three. Let E' be the set of uncoloured edges in G. 

Let the class of any fan F to be constructed in G be the ordered pair ( c,d ), 
c < d, if F is either a (c,d)-fan or (d,c)-fan. So, there would be 0( A 2 ) classes of fans 
m G at any time. The set of fans belonging to class (c,d), will be denoted by C cd . 

The following procedure assigns a colour from {1, . . . , A + 1} to to each edge 

of^: 


Procedure Solve(F') 

Loop through the following steps: 

Step 1: Construct fans at every uncoloured edge in G. Perform a fan operation on 
every local fan constructed. 

Step 2: Of the fan-classes of G, find the class with the maximum cardinality; let that 
class be (c,d). 

Step 3: Form a conflict graph H as follows: corresponding to each c-d path of G, let 
there be a vertex in Ft , and let two vertices x and y of 'H be adjacent iff the c-d paths 
corresponding to x and y end at the same fan in G. A c-d path P in G is said to end 
in a fan F, if P is either H(F ) or F(F). 

Remark: Each vertex of H can have (at most) two neighbours; a c-d path may end in 
a fan at both its endpoints. So, H consists of chains and cycles. Note that each fan in 
C c d corresponds to an edge in H. 

Step 4: From every odd cycle of Ti remove one vertex. Rank the remaining graph, 
and for vertices of odd rank swap colours c and d in the corresponding c-d paths of G . 
Remark: Clearly, colours are swapped in at least two thirds of the c-d paths of G. 
Hence, at least two thirds of what were (c, d) or (d. c)-fans have now become (c, c) or 
( d , d)- fans, which are clearly local. (See the proof of Vizing’s theorem in Section 7.1.1.) 
For a fan F € C cd , if there are only two distinct c-d paths among H(F), M(F) and 
T(F), then colours c and d have been swapped in at most one of those two paths. On 
the other hand, if H(F), M(F ) and T(F) axe all distinct, and hence, E(F) f T(F), 
then T( F) also may have had its colours swapped, along with at most one of H(F) and 
M(F); but, this does affect the new fan obtained, because, it does not contain t(F). 
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Step 5: Perform a fan operation on every local fan in G. 

Remark: At least 2/3 * \C c i\ fans are now solved. But, swapping colours c and d in 
some c-d paths of G may have destroyed some fans not in C c d- 


Lemma 7.2 The procedure call “Solve(E')” assigns a colour from {1,...,A + 1} to 
each edge of E' in 0( A 2 log 2 n) time with n + m processors on an EREW PRAM. 

Proof: By Lemma 7.1, Step 1 of the loop takes 0(log A) time. The rest of the steps 
are dominated by an instance each of sorting 0(n) items and ranking a list of 0(n ) 
length. Hence, clearly, each iteration of the loop can be executed in O(logn) time time 
with n + m processors on an EREW PRAM. Each iteration reduces the total number 
of uncoloured edges by a fraction of 0(1/A 2 ). In other words, each iteration reduces 
the number of uncoloured edges to 0((A 2 — 1)/A 2 ) of what had been there before. So, 
the number of iterations required to colour all edges is 0(logn/(log(A 2 /(A 2 - 1)) = 
0(A 2 log n). Hence the lemma. D 

Now we present a procedure that edge-colours a general Graph: 


Procedure Edge-Colour(G, r, A(G)) 


Input: A graph G\r is a parameter, r < logA(G). 

Output: A ¥(r, A(G))-edge-colouring of G\ f(r, A) > A + 1. 

Step 0: If A(G) = 2, then G is collection of lists and cycles. Obtain a 3-edge-colouring 
of G. exit. 


Step 1: Decompose the edge set E of G into two disjoint sets E 0 and E 1 so that 

(a) in Go = (V,E 0 ), for some s 6 V, every vertex other than s has a degree of at most 

fAM] ; an d s has a degree of at most + 1, and 

(b) in Gi = (V,Ei) every vertex has a degree of at most f^'l- 
Remark: The details of this step are given later in this chapter. 


Step 2: 


For an edge e = {s,t} in Go, let G' 0 - Go e. 


Now both Gq and Gi are 


f^M] -degree graphs. 

Execute the following two procedure calls in parallel. 
Call Edge-Colour(Go,r - 1, f^D- 
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Call Edge-Colour^, r - 1, T^]). 

Remark: Now both G'q and G\ are $(r — 1, [— ip^"|)-edge-coloured. 

Step 3: Extend the edge-colouring of G' 0 to Go (e is the only uncoloured edge of Go) 
as follows: create a fan for e with t as its head; perform fan and chain operations, as 
necessary, and give e a colour in {1, . . . ,f (r-1, Note that by Observation 7.1 

this is feasible. Now both Go and Gi are ^(r - 1, [^^])-edge-coloured. 

Add $(r — 1, 1^1) to every colour used in G\. 

Remark: Now G is 2'$ , (r — 1, )-edge-coloured. For 1 < x < 2 $(r — 1, 
let j E x be the set of edges now coloured x in G. 

Step 4: If r > 1, then skip this step; otherwise, proceed. For A(G) + 1 < x < 
2 $(r - 1, f^l), uncolour the edges of E x . 

For each colour x in {A + 2,. . . ,2$(r - 1, [y])} do: 

Form a conflict graph C = (E x ,£) where {e, /} G £ for some e,f G E x iff 
there is a path of length less than 3 between e and / in G. The maximum 
vertex degree a of C is 0(A 2 ); note that there can be at most one edge of 
colour x at any vertex. Find an (a + 1) vertex-colouring a of C. Let E 3 X 
denote the set of edges e in E x such that a(e) = i. For 1 < j < a + 1, in 
turn call Solve(J^). 


Now we describe Step 1 in detail: 


Substep 1.1: If G is not Eulerian, then add a vertex u that is adjacent to every odd 
degree vertex in G. Let G' be the augmented graph thus obtained. If G is Eulerian let 

G' = G. 

Remark: Since the number of odd degree vertices in any graph is even, G' is Eulerian. 

Substep 1.2: Find an Euler tour in G'. If G ^ G' let this tour start at z = u, 
otherwise, let it start at some z G V. Label the edges of the tour 0 and 1 alternately; 
the first edge is labelled 0. Let Go and G x be the subgraphs of G induced by the edges 

labelled 0 and 1 respectively. ^ 

Remark: If G' has an even number of edges, then clearly, Go and Gi are I" 2 1 
graphs. But if G' has an odd number of edges, then both the first and the last edges 
of the Euler tour obtained would be labelled 0, and hence, z would have incident to it 
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two more label-0 edges than label-1 edges. If G £ G' and hence z £ V , then this would 
not matter, because neither G 0 nor Gi would contain 2 , and G 0 and G\ would still be 
[^-degree graphs. 

If G = G and G' has an odd number of edges, then the vertex z has degree 
2 + 1 in Go and j — 1 in G\. Every vertex other than z has degree in both 

Gq and G±. 


Step 4 is executed only when r < 1. Clearly, the number of colours used by the 
procedure call “Edge-Colour(G, 0, A(G))” is $(0,A(G)) = A(G) + 1. 

So, the number of colours used by the procedure call “Edge-Colour(G, r, A (G))” 
is 



Thus, for d > 1, the procedure call “Edge-Colour(G, [log d\ - 1,A)” gives a (A 4- d)- 
edge-colouring of G. 

Let T(r, A) be the time taken by the procedure call “Edge-Colour(G, r, A(G))” . 
A 3-edge-colouring of a 2-degree graph can be found (Step 0) in O(log(log* n)) time 
with 7i processors on an EREW PRAM [48]. An Euler tour in an Eulerian graph can be 
found (Step 1) in O(log n) time on a CRCW PRAM [8], and hence, in O(log 2 n) time on 
an EREW PRAM, with 0(n + m ) processors. Step 2 takes T(r, [A] ) time. Step 3 takes 
<3 (log n) time with n + m processors on an EREW PRAM. An (a + l)-colouring of an 
a-degree graph can be found in 0(aloga(log* n + logo) time, with n processors, on an 
EREW PRAM [32]. Hence, by Lemma 7.2, Step 4 takes 0(A 2 log A(log* n + log A) + 
a A 2 log 2 n ) = 0(A 4 log 2 n) time with n processors on an EREW PRAM. Step 4 is ex- 
ecuted only when r < 1. Thus, after r becomes zero, Step 4 dominates the time taken, 
and hence, T( 0, A) = 0(A 4 log 2 n). Therefore, for r > 1, for an appropriate constant 


T(r,A) = r(r-l,|y]) +clog 2 n 
T^r-2, |"j"P -|- 2c log 2 n 


< 

< •• 
< T 


(°# 1 ) 


4- rclog 2 n 



Chapter 7. Edge-Colouring of Graphs 


96 


= 0((A/2 r ) 4 log 2 n + rlog 2 n) 

Setting r — [logdj — 1, thus we have the following theorem: 

Theorem 7.1 A graph can be (A + d)-edge-coloured in O((log d+(A/d) 4 ) log 2 n ) time 
with 0(n + m) processors on an EREW PRAM. 

7.2 A 1+f -Edge-Colouring, 0 < e < 1 

An O(log A log(log* n)) time, nA 1+£ processor EREW PRAM algorithm for 
^ -edge-colouring an arbitrary graph is now presented. We employ a technique 
similar to the one that Gabow and Kariv [28], and Lev, Pippenger and Valiant [69] use 
to 2 ^ ogA l -edge-colour a bipartite graph. 

Consider a graph G = (V, E) with V = { 1 , . . . , n } . Assume that G is given in 
adjacency list representation. Let N(v) be the adjacency list of v e V, and [v, w] the 
entry in N(v) corresponding to an edge {v, w} e E. 

Visualise each edge of G as having two halves, one half for each end vertex. 
Suppose each vertex of G has a palette {0, . . . , A - 1} of colours. Let each vertex assign, 
from its palette, a unique colour to every half-edge incident with it. This can be done 
as follows. For each v € V, in parallel, rank N(v), starting from zero. Let r v (w) be the 
rank of [v, «;] in N(v ); r v (w ) < A — 1. For an edge {u, w } £ E , let ao(v, w) = r v (w ) be 
the colour assigned to {n, w) by v. 

Thus, each edge receives two colours, one for each half; these two colours need 
not be the same. We now proceed to modify the present colouring so that both halves 
of an edge become coloured the same; but this will cost us some extra colours. 

For an n x A array Ao[l...n,0...A - 1], let A 0 [t’,yo(v,iy)] = [v,w]. Let b 
and B, b < B < A, be two parameters to be chosen later. Suppose the colours in 
the palette are represented as base-6 integers. Such a representation will have 
digits. 

The algorithm now proceeds in phases. In each phase, at a high level, 

for each edge e of G, the least significant digits of the two colours on edge e, which 
are two possibly different base-6 numbers, are replaced by a base-5 number, the same 
in both cases. Also, at the end of each phase, we perform a "shift-right-with-rotate 
on every colour of G, so that each phase deals with a different digit of the original 
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colouring. After phase, the graph would be validly edge-coloured, but each of 

^he r log 6 1 digits of each colour will be a base - B number. 

Now we describe the phases more formally. With respect to the phases, we 
maintain the following invariants: at the beginning of the i-th phase, 1 < i < Q 

is (B/b)'- 1 A-edge-coloured; the colours are from {0, . . . , (f?/f>) i-1 A - 1}. For a colour 
c, let Ri—i(c) be a -digit representation of c, where the i — 1 most significant digits 
are all in base B and all lesser significant digits are in base b. For an edge { v, in} £ E . 
let cr;-i(v, w) denote the colour assigned to {u,tn} by v, at the beginning of the i-th 
phase. Colours 0 V-i(v, w) and <7j_ i(tn,u) agree on the (*—1) most significant digits. The 
adjacency list of G is given in an n x (R/&)‘ -1 A array A,_i[l. ..n,0. ..(B/6) i_1 A-l]; 
A t -_x[u,x] contain [u, in] iff x = w). Here, we do not assume any specific order 

of linking for entries in an adjacency list. Note that the invariants are satisfied at the 
beginning of the first phase. 

The i-th phase, 1 < i < , is now described: 


Phase i 

Step 1: Divide the set of (J9/&)' -1 A presently used colours into classes of size b each; 
a colour c is in the j-th class if [c/6J = j. Two colours ci and c 2 are in the same 
class iff iZ,_i(ci) and i?,_i(c 2 ) agree on all but the least significant digit. There are 
Qi = (B'~ l /b l )A classes. 

Form a set Vi of vertices as follows: for each v € V, and 0 < j < Qi- 1, let u,-j 
be in V]; that is, Vi has Qi copies of each member of V . 

Construct a graph (?,- on V as follows: for each adjacency list entry [u, u>] of 
G, if w) belongs to the j-th class and <7i-i(w,v) belongs to the k- th class, then 

place an edge between v t -j and w^. 

Remark: There is a one to one correspondence between the edges of G and Gi- Also, 
two vertices and w^k are adjacent in Gi only if <7i-\(v,w) belongs to the j-th class. 
Thus, Gi is a 6-degree graph. Since, all the edges of class j (an edge is said to be in 
the same class as its colour) incident at v € V, appear in a subarray of size b of the 
u-th row of Ai- 1, the adjacency lists of Gi can be constructed in O(log&) time on an 
EREW PRAM, provided a processor is stationed at each cell of A;_i . Note that for 
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each adjacency list entry [u,uj] of G. the class number of cr,_ 1 (u’,u) can be determined 
in 0(1) time. 

Step 2: Edge-colour Gi using colours from {0,. . . ,5 — 1}. We assume that there 
is a T(5,6 ) time algorithm for 5-edge-colouring a 6-degree graph. Suppose, the edge 
{vij,Wi t k} gets colour l. Let cr;(v,w) = (l,j) and cr t (w,v) = (l,k). 

Remark: Clearly, i2,_i(er,(v,u;)) and R,_i(cr,(to, u)) agree on the i-most significant 
digits. The total number of colours present in G now is (5/6)* A. 

Step 3: Copy the adjacency lists of G into .4,[1 . . . n.O. . . (B/b)' A — 1], an nx A(5/6)‘ 
array: if Oi(v,w) = x, then let A,[v,x] contain [v,w]. 

The algorithm described above finds a $(5,6) = 5 ^° sA/,log ^ -edge-colouring 
in 

0 ( ! gy (3W)+log ‘ ) ) 

time with n$(5,6) processors on an EREW PRAM; here T(5,6) is the time required 
to 5-edge-colour a 6-degree graph. 

Let Gl = (■£> £) be the line graph of G: each edge of G corresponds to a 
vertex of G L , and two vertices of Gl are adjacent iff the corresponding edges of G are 
adjacent in G. Any fc- vert ex- colouring of Gl is a fc-edge-colounng of G. Since each 
edge of G is adjacent to at most A .(G) - 1 edges at either end, G L has a maximum 
vertex degree of 2A(G) - 2. Thus, a (2A(G) - l)-edge-colouring of G can be found in 
0(A log A(log(log* n) + A)) time with n + m processors [32] on an EREW PRAM. 

Hence, a $(26,6) = (26) ^ og log ^ -edge-colouring of an arbitrary graph can 
be found in 0(61og A(log(log* n) + 6) time with with n$(26,6) processors on an EREW 

PRAM. Note that $(26,6) = (26)^/'°^ is at most A^W, for 6 < 2> g(A ' 1) - 1 - 
Let 6 = 2 2 A; 0 < e < 1; then, £ = 2/ log 6. Thus, we have the following 

theorem: 

Theorem 7.2 A graph can te A'*- edge-coloured in 0( l»gAU*W»)) time with 
nA 1+€ processors on an EREW PRAM; 0 < € < 1. 
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7.3 4-Edge-Colouring a 3-degree Graph 

Let G = (V, E) be a 3-degree graph. Each edge of G is adjacent to at most 4 
edges, at most two at either of its end-vertices. Thus, the line graph Gl = (E,£) has a 
maximum vertex degree of 4. Moreover, Gl cannot have a 5-clique; this is because, five 
mutually adjacent vertices in Gl must correspond to five mutually adjacent edges in 
G, which must then all be incident with the same vertex of G. Hence, Gl is a 4-degree 
Brooks’ graph. Obtain a 4- vertex colouring of Gl , and this would correspond to a 
4-edge colouring of G. 

By Corollary 5.2, we get the following theorem. 

Theorem 7.3 A 3-degree graph can be 4-edge-coloured in O(logn) time with pro- 
cessors on an EREW PRAM. 

The above theorem improves the result of [60] where a 3-degree graph is 4- 
edge-coloured in O(logn) time with n processors on a CRCW PRAM. 
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Chapter 8 

Conclusions 


Several parallel algorithmic issues related mainly to sorting, merging and 
graph colouring were considered in this thesis. 

In chapter 2, we addressed the question, whether k - way merging is exactly as 
hard as sorting k numbers and merging 2 arrays of n elements with the same processor 
advantage in each case, i.e., 

&_Way_Merge(b,p) = 0(Sort(&, ^) + Merge(b,p)) (8.1) 

We presented an algorithm schema for k - way merging. For comparison- based versions 
of sorting, merging and k- way merging, using the schema we showed that Eq. 8.1 is 
satisfied, on all of EREW, CREW and CRCW PRAM models. For integer versions, on 
a CRCW PRAM, when we use the best known integer sorting and merging algorithms 
in our schema, and the keys are integers of (logn) fi ^ bits, again, Eq. 8.1 is satisfied. 
However, for improvements to be anticipated in sorting and merging algorithms, as well 
as, for smaller integer keys, our schema may not satisfy Eq. 8.1. But, we conjecture 
that the equation will nevertheless hold. Theorem 2.12 can be seen as an evidence in 
favour of this conjecture. 

A merge sorting algorithm that runs in to .^(b. n/ioga)) ste ps, with nQ 
processors, 2 < a < n, on the parallel comparison model was presented m Chapter 3. 
This algorithm comes close to achieving the lower bound on sorting on the 

parallel comparison model; moreover it is very different from Cole’s merge sort. It 
remains to be seen if this algorithm can be improved into matching the lower bound. 
On a CREW PRAM, a corresponding question would be whether O(logn) time could 
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be achieved, thereby giving a simple sorting algorithm that does not use pipelining. 

A still harder problem would be designing a simple O(logn) time optimal comparator 
network for sorting. 

An G((loglogn)log*(log* n)) time optimal algorithm on the TOLERANT 
CRCW PRAM, for 3-colouring general rooted forests, was presented in Chapter 4. 
Here optimality was achieved at a loss of speed: a rooted forest can be 3-coloured in 
0 (log (log* n)) time on an n processor CREW PRAM. Moreover, our algorithm depends 
on the approximate compaction problem that has a parallel complexity of 0( log log n). 
It remains an open question whether an o(loglogn) time optimal algorithm for 3- 
colouring general rooted forests can be found. 

In Chapter 4, we also show that a TOLERANT PRAM of size N with a linear 
address space, can be slowed down by any factor A = ft(loglogJV), with no asymptotic 
increase in space or cost. Unlike other CRCW PRAM models, TOLERANT is not 
known to be self-simulating in general [39]; and probably it is not — but this remains to 
be proved. 

An O(logn) time optimal EREW PRAM algorithm for Brooks’ colouring of 
bounded degree graphs was presented in Chapter 5. In addition, a combinatorial result 
was obtained as evidence to suggest that an asymptotically faster algorithm may exist. 
But, how to use this combinatorial result in an algorithm is not clear yet. The lower 
bound problem on Brooks’ colouring is also open; the trivial U(log(log* n)) time lower 
bound on the distributed model is the only one currently known. 

In Chapter 6, we argued that the problem of colouring an interval graph with a 
known interval representation (IGC) may not be in NC 1 . In this context, the following 
questions arise (i) whether IGC is FL-complete for Ntf-reductions, (») whether IGC 
can be solved in 0(log n/ log log n) time on a CRCW PRAM with a polynomial number 
of processors. Both of these are as yet open. However, we gave an O(logn/loglogn) 
time, polynomial processors, algorithm for IGC, for the special case of the chromatic 

number being 0((logn) 1-4 ), 0 < e < 1- 

Two parallel algorithms for wasteful edge colouring were presented in Chap- 
ter 7. Finding a minimal edge colouring is computationally difficult even sequentially; 
The problem of (A + l)-edge-colouring a general graph, on the other hand, has a poly- 
nomial time sequential solution, but, does not render itself to easy paraHelisation; it 
has been an open question for long whether this problem is in NC or not. 
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